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Recent Advances in the Physiology and 
Genetics of Amino Acid-Producing Bacteria 

Mike S. M. Jetten 1 and Anthony J. Sinskey 2 

department of Microbiology and Enzymology, Kluyver Laboratory for Biotechnology, Delft University 
of Technology, Julianalaan 67, 2628 BC Delft, The Netherlands; department of Biology, Room 
68-370, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139 

ABSTRACT: Corynebocterium gtutamicum and its close relatives, C. flavum and C. lactofermentum, have been 
used for over 3 decades in the industrial production of amino acids by fermentation. Since 1984, several research 
groups have started programs to develop metabolic engineering principles for amino acid-producing Corynebac- 
terium strains. Initially, the programs concentrated on the isolation of genes encoding (deregulated) biosynthetic 
enzymes and the development of general molecular biology tools such as cloning vectors and DNA transfer 
methods. With most of the genes and tools now available, recombinant DNA technology can be applied in strain 
improvement. To accomplish these improvements, it is critical and advantageous to understand the mechanisms 
of gene expression and regulation as well as the biochemistry and physiology of the species being engineered. This 
review explores the advances made in the understanding and application of amino acid-producing bacteria in the 
early 1990s. 

KEY WORDS: Corynebocterium glutamicum, lysine, glutamate, pyruvate kinase, aspartokinase, metabolic engineering. 



I. INTRODUCTION 

Amino acids have a wide variety of applica- 
tions. Monosodium glutamate (MSG) and glycine 
are used as flavor enhancers in food, while some 
essential amino acids such as lysine, threonine, 
tryptophan, and methionine are supplemented in 
animal feed (Flodin, 1993). Phenylalanine and 
aspartate are used in the production of artificial 
sweeteners, and many other amino acids serve as 
therapeutic agents in nutritional and metabolic 
disorders or as raw material for the synthesis of 
cosmetics, toothpaste, shampoo, and detergents 
(Martin et al., 1987). 

Amino acids are produced via chemical syn- 
thesis, hydrolysis of natural proteins, and bacte- 
rial fermentations (Kinoshita, 1985). Only a lim- 
ited number of bacterial species are used to obtain 
production strains for the industrial manufacture 
of amino acids (Aida et al, 1986). The detailed 
genetic information and well-developed host-vec- 
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tor systems of Escherichia coli made this organ- 
ism very suitable to construct amino acid-excret- 
ing strains via recombinant DNA technology 
(Beppu, 1986). In this way, E. coli strains were 
constructed that produce considerable amounts of 
threonine, tryptophan, or valine (Aibaet al., 1980, 
1982; Beppu,. 1986; Debadov, 1983; Miwa et al., 
1983). Strains secreting proline, threonine, or iso- 
leucine were derived from Serratia marcescens 
by using classic mutagenesis and screening pro- 
cedures in combination with transduction of de- 
regulated genes (Aida et al., 1986; Nakamori et 
al., 1986). Several Bacillus species have been 
used to breed strains that produce aromatic amino 
acids (Shiio, 1986). 

However, the bacteria most widely used for 
the industrial production of amino acids belong to 
the Corynebacterium-Brevibacterium group of 
coryneform bacteria (Archer et al., 1989; Batt et 
al., 1985a; Eikmanns et al., 1993; Follettie et al., 
1991; Jetten et al., 1993; Kinoshita, 1985; Liebl, 
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1991; Wohlleben etal., 1993). Recent DNA-DNA 
hybridization studies show that there are only 
very minor differences between the B. flavum, 
B.. lactofermentum, and C. glutamicum species 
(Liebl et al, 1991; Eikmanns et al., 1991). These 
differences do not justify the use of different ge- 
nus names, and it has been proposed to consider 
them as one C. glutamicum species (Liebl et al., 

1991) . Throughout this paper, we therefore will 
use the name Corynebacterium. Corynebacterium 
species are Gram-positive, irregular, non- 
sporulating, strictly aerobic, food-grade organ- 
isms that have long occupied a central role in the 
manufacture of amino acids, such as glutamate, 
lysine, threonine, tryptophan, and phenylalanine 
(Kinoshita, 1985). The world-wide production in 
1992 of the flavor enhancer MSG and the food- 
supplement L-lysine/HCl by C. glutamicum was 
estimated at approximately 800,000 and 400,000 
metric tons, respectively (Randell and Andreas, 

1992) . 

In this review, recent findings that contrib- 
ute to the understanding of the relevant physi- 
ology, biochemistry, and molecular biology of 
C. glutamicum and closely related organisms, 
with particular emphasis on the biosynthetic 
pathway for the aromatic and aspartate-derived 
amino acids, are reviewed. Where possible, 
reports about other organisms and amino acids 
are referenced. 



IU BASIC CLONING TECHNIQUES FOR 
CORYNEBACTERIUM GLUTAMICUM 

Although the production of amino acid by 
C. glutamicum and close relatives has been known 
for more than 30 years, the first reports on the 
metabolic engineering of C. glutamicum only 
appeared in 1984 (Miwa et al., 1984; Santamaria 
et al., 1985). Since that time, several research 
groups have independently initiated research pro- 
grams focusing on the development of metabolic 
engineering tools for Corynebacterium species 
(Jetten et al., 1993; Martin et al., 1987; Wohlleben 
et al., 1992). Essential prerequisites for the meta- 
bolic engineering of Corynebacterium are the 
availability of vectors derived from endogenous 
plasmids and efficient DNA transfer systems. The 
availability of these tools allows the cloning, ex- 



pression, disruption, or replacement of genes. By 
using these procedures, the metabolic processes 
that are of biotechnological importance can be 
optimized and used for the construction of more 
efficient production strains, as has been shown 
for E. coli and S. marcescens (Beppu, 1986). 



A. Corynebacterium Plasmids and 
Cloning Vectors 

Numerous Corynebacterium strains have been 
screened for the presence of endogenous plas- 
mids (Martin et al., 1987; Sonnen et al., 1991; 
Wohlleben et al., 1992). Small, cryptic plasmids 
were independently detected in various Coryne- 
bacterium strains. Some of the plasmids are very 
similar and to avoid confusion have been grouped 
into four different families (Sonnen et al., 1991). 
An overview is presented in Table 1. 

The complete nucleotide sequences of repre- 
sentatives of three plasmid families (pBLl , pSRl , 
and pNG2) are available at the moment (Archer 
and Sinskey, 1993; Filipula et al., 1986; Serwold- 
Davis et al., 1990; Yamaguchi et al., 1986). pBLl 
(4.45 kb, isolated from C. lactofermentum) and 
P HYM1519/pSRl (3.1 kb, isolated from 
C. glutamicum) families are most popular among 
the research groups involved in the metabolic 
engineering of Corynebacterium strains (Archer 
and Sinskey, 1993; Eikmanns et al., 1991c; Jetten 
et al., 1994a; Miwa et al., 1985; Sandoval et al., 
1985). Both plasmids have been used to construct 
a multitude of cloning vectors, which replicate in 
Corynebacterium strains and/or E. coli. Although 
a few endogenous antibiotic resistance genes 
(streptomycin, tetracycline, and trimethoprim) 
have been found in Corynebacterium (Batt et al., 
1985b; Katsumataetal., 1 984; Takagi etal., 1986; 
Takedaet al., 1990), several heterologous, more 
convenient antibiotic resistance markers (erythro- 
mycin, chloramphenicol, and kanamycin) are ex- 
pressed efficiently in Corynebacterium and there- 
fore commonly used (Cardenas et al„ 1991 ; Chan et 
al., 1991; Guerrero et al., 1992; Martin et al., 1987; 
Santamaria etal., 1987; Serwold-Davis etal., 1987; 
Yoshihama et al., 1985). 

A comprehensive study on the host range of 
the various plasmids is still lacking, but of all the 
Corynebacterium plasmids tested, pNG2 and its 
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relatives have the widest host range (Serwold- 
Davis et al., 1987). Recently, however, a broad 
host-range derivative of pSRl has been described 
(Archer and Sinskey, 1993) that is able to repli- 
cate in both E. coli and Corynebacterium species. 



B. Mobilizable Vectors 

Recently, a special class of integration vec- 
tors has been developed and is now very fre- 
quently used for gene disruption and gene re- 
placement experiments in Corynebacterium 
(Schafer et al. t 1990; Schwarzer and Puehler, 
1 99 1 ). These so-called suicide vectors lack a func- 
tional Corynebacterium replicon, and can there- 
fore only be maintained when integration occurs 
via homologous recombination into the chromo- 
some. To employ these vectors, very efficient 
DNA transfer via transconjugation or electro- 
poration is required. 

C. DNA Transfer Methods 

At the moment, several methods for the intro- 
duction of cloning vectors and exogenous DNA 
into Corynebacterium are available. These include 
protoplast transformation, transduction, electro- 
poration, and transconjugation. 

I. Protoplast Transformation 
and Transduction 

Protoplast transformations developed for 
amino acid-producing Corynebacterium were 
based on the polyethyleneglycol (PEG)-medi- 
ated standard procedures used for other Gram- 
positive bacteria (Katsumata et al., 1985; 
Santamariaet al., 1985; Yoshihamaetal., 1985). 
The efficiency of this procedure is rather low 
(10 5 transformants per milligram of DNA), and 
the regeneration of protoplasts is often problem- 
atic. Several reports have claimed that the proto- 
plasts generated with the PEG method can also 
be used in protoplast-fusion experiments 
(Karasawa et al., 1986). The improvements ob- 
tained with this method are rather dubious, and 
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the method appears to be very inefficient. So far 
only a few reports have appeared that described 
phage-based transduction of Corynebacterium 
strains (Momose et al., 1976; Ozaki et al., 1984; 
Patek et al., 1988; Sachez et ah, 1986). How- 
ever, the discovery and characterization of both 
temperate nonlytic phage phiGAl and virulent 
corynephages CL31 and Cog may be a good 
starting point to develop new tools for the use of 
transduction in C. glutamic urn (Sonnen, 1992; 
Sonnet et al., 1990a,b; Trautwetter et al., 
1987a,b). 



2. Electroporation 

Instead of protoplast transformation, electro- 
poration and transconjugation are now used as 
more convenient and much more efficient meth- 
ods for transformation. Different research groups 
have each developed their own procedure for 
electroporation (Bonamy etal., 1990; Bonnassie 
et al., 1990; Dunican and Shivnam, 1990; 
Follettie et al., 1993; Hyanes and Britz, 1989, 
1990). Normally, cells used for electroporation 
are grown in rich medium for 3 to 5 h, harvested 
at the early logarithmic phase (OD: 600 nm be- 
low 2), washed extensively with low-salt buffers 
containing 5 to 20% glycerol, and concentrated 
to approximately 10 10 cells per milliliter before 
use. To increase the efficiency of electroporation, 
several compounds (penicillin G, ampicillin, 
isonicotinic acid, or Tween 80) can be added to 
the growth medium to make the cells more per- 
meable for exogenous DNA (Bonassie et al., 
1990; Follettie et al., 1993; Haynes and Britz, 
1990; Wolf et al., 1990). Generally, 100-u.l 
aliquots of the cells are equilibrated with 0.1 to 
1 |Xg of DNA and subjected to a field strength of 
12.5 kV/cm for about 3 to 8 msec. Thereafter, 
the cells are immediately mixed with 0.9 ml of 
an isotonic medium, incubated for 1 to 1.5 h to 
express the antibiotic marker, and plated on a 
sorbitol/sucrose/brain heart infusion medium. 
Transformants generally appear within 48 h. The 
maximum number of transformants obtained with 
this method varies between 4 and 5 x 10 7 per 
microgram of DNA (Haynes and Britz, 1989; 
Follettie et al., 1993). 
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3. Transconjugation 

Another very efficient method for the intro- 
duction of DNA into Corynebacterium is the 
transconjugation of mobilizable suicide plasmids 
from E. coli to various Corynebacterium strains 
(Schafer et al. f 1990, 1994). This method is fre- 
quently used for gene disruption and gene re- 
placement in C. glutamicum (Jager et al., 1992; 
Labarreetal., 1993; Reyes etai., 1991; Schafer et 
al., 1990; Schwarzer and Puehler, 1991). The 
system depends on the mobilization function of 
the broad host-range plasmids derived from RP4 
to transfer plasmids from specialized £. coli strains 
(e.g., SI 7-1) to Gram-positive Corynebacterium 
strains. Special strains with defined deletion mu- 
tations, which lack all vector and antibiotic resis- 
tance marker sequences, can be constructed by 
combining the techniques of gene disruption and 
replacement. Strains constructed in this manner 
are very likely to be considered safe with respect 
to environmental release or industrial use. Sev- 
eral examples of the application of these tech- 
niques, summarized in Table 2, are reviewed in 
this article. 

Numerous observations made during trans- 
formation and transconjugation experiments indi- 
cate the presence of restriction or modification 
systems in Corynebacterium strains (Follettie et 
al., 1991; Liebl, 1991; Wohlleben et al., 1992). 
The problems caused by these systems for clon- 
ing experiments can be overcome by the use of 
restriction-deficient host strains, by heating the 
cells for 9 min at 50°C or exposing the cells to 
environmental stress such as pH, ethanol, or SDS 
(Bonassie et al., 1990; Liebl, 1991 ; Follettie et al., 
1991; Schafer et aL, 1990, 1994). 

D. Gene Expression 

A large number of both homologous and het- 
erologous genes have been successfully expressed 
in Corynebacterium, although comprehensive 
knowledge about translation and transcription 
initiation signals in this organism is still lacking. 
To date there are no indications that translation 
initiation signals from other bacteria do not func- 
tion in Corynebacterium. The information avail- 



able on endogenous transcription initiation is only 
slowly emerging. So far only about ten transcrip- 
tion start points have been mapped (Follettie et 
al., 1991). Visual inspection indicated that there 
are major differences between these promoters 
and those from, for example, E. coli and Bacillus 
subtilis. The most conserved sequence in the pro- 
moter region was a TGTGT motif approximately 
20 base pairs upstream of the transcription initia- 
tion site (Follettie et al., 1991). A TATAAT se- 
quence analogous to other bacterial promoters 
around the -10 position was only moderately con- 
served. Recently, McCormick et al. (1993) pre- 
sented a study in which the TGTGT motif of the 
thrC promoter was targeted by site-directed mu- 
tagenesis. In this way, 14 single base-pair sub- 
stitutions in the TGTGT region were obtained 
and analyzed by comparing the differences of the 
P-galactosidase activity of these thrC-lacZ con- 
structs. Of the 14 mutations, only three had a 
significant effect on promoter strength; a change 
of the TGTGT to TGCTG resulted in a two-fold 
increase in activity, whereas TGTCT or TGTTG 
resulted in a 50% decrease in promoter strength. 
These results are encouraging, but DNA 
footprinting and characterization of more 
C. glutamicum promoters will be necessary to 
obtain a full understanding of transcription initia- 
tion in this organism. 

The efficiency and regulatory properties of 
very strong E. coli promoters such as were 
found to exhibit similar properties in Corynebac- 
terium (Tsuchiya and Morinaga, 1988). Using the 
laclq/Puc system, several inducible, high-expres- 
sion vectors were constructed (Eikmanns et al., 
1991c; Liebl et al., 1992; Gubler, unpublished 
results). In contrast to highly efficient promoters, 
strong transcription termination signals of E. coli 
do not seem to function in Corynebacterium 
(Wohlleben et al., 1993). 



E. Isolation of Corynebacterium 
Glutamicum Genes 

The availability of the aforementioned clon- 
ing systems in Corynebacterium has greatly fa- 
cilitated the cloning of genes involved in both 
intermediate metabolism and several amino acid 
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TABLE 2 • ~ , , ■ 

Genes Disrupted or Integrated Via Transconjugation in C. glutamicum 



Target gene 



Disrupted 
or integrated 



as* 


Intparated 




Di^ruoted 




Intearated 




Disrupted 


cspA 


Disrupted 


dapE. 


Disrupted 


ddh 


Disrupted 


gdh 


Disrupted 


gltA 


Disrupted 


horrP'-thtC 


Integrated 


HvA 


Disrupted 




Disrupted 


lysA 


Disrupted 


ppc 


Disrupted 


pyk 


Disrupted 



Phenotype 

Increased internal lysine 
No lysine production 
Increased lysine production 
No lysine production 
No S-layer formation 
Decreased lysine production 
Decreased lysine production 
No apparent effect 

No apparent effect 

Increased threonine production 

Isoleucine auxotrophy 

No threonine degradation 

Lysine auxotrophy 

No apparent effect 

Decreased lysine production 



Ref. 
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Labarre et al., 1993 
Bormann et al., 1993 
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Reinscheid et al., 1994 
Boles et al., 1993 
Jetten, unpublished 
Schwarzer and Puehler, 1990 
Gubler et al., 1994a 
Peterschwinde et al., 1993 
Gubler, 1994b 



biosynthetic pathways. Table 3 summarizes the 
genes cloned from C. glutamicum. Three strate- 
gies have been followed for cloning Corynebac- 
terium genes: heterologous complementation, 
homologous complementation, and screening of 
libraries with degenerate oligonucleotides. 

Complementation of E. coli auxotrophs by a 
gene bank of C. glutamicum made in E. coli using 
either plasmid or cosmid vectors has been used 
successfully in the cloning of biosynthetic genes 
for tryptophan, lysine, isoleucine, and threonine, 
and of genes encoding the glyoxylate bypass* and 
PEP carboxylase (Table 3). A second approach 
has been the direct cloning of genes in 
C. glutamicum by complementation of amino acid 
auxotrophs of this organism. In this way, several 
genes of the phenylalanine and threonine path- 
way have been obtained (Table 3). 

A third strategy, which has been used only 
sporadically, is the use of nucleotide probes de- 
rived from either N-terminal or internal peptide 
sequences of purified proteins. Recently, Gubler 
et al. (1994) used a combination of the afore-, 
described methods to clone the pyruvate kinase 
gene of C. glutamicum. On the basis of homology 
between published pyruvate amino acids se- 
quences, two degenerate oligonucleotides were 
designed and used to synthesize a 126-bp internal 



fragment of the C. glutamicum pyruvate kinase 
gene (pyk) via a polymerase chain reaction. This 
internal fragment was cloned in an RP4-derived 
suicide vector and used to generate defined Pyk - 
mutants of C. glutamicum. The pyk gene thereaf- 
ter was obtained via complementation of the 
mutant on minimal medium containing ribose as 
carbon source. 

Genes obtained via one .of the aforemen- 
tioned methods can be used to improve produc- 
tion strains by raising the activity of rate-limit- 
ing enzymes or removing the feedback regulation 
of critical enzymes. Use of such genes also al- 
lows for specific probes to be utilized in order to 
determine rate-limiting steps in the biosynthesis 
of a product of interest. 

The biosynthetic pathway of several amino 
acids in C. glutamicum is used to illustrate the 
aforementioned technology. 



III. PRECURSOR SUPPLY 

A. Introduction 

Although great progress has been made in 
the understanding of the molecular organization 
of biosynthetic genes (Jetten et al., 1993; 
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Eikmanns et al., 1993), relatively little is known 
about the enzymes involved in the central car- 
bon metabolism (Figure 1) of C. glutamicum 



(Jetten et al., 1993, 1994a, b; Krulwich and 
Pellicone, 1979). With the limited information 
available, several attempts have been made to 
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FIGURE 1. Central metabolism of C. glutamicum, OAA, ox- 
aloacetate; PEP, phosphoenol pyruvate; PTS, PEP: sugar 
phosphotransferase system; gota, glutamate: OAA transami- 
nase; icd, isocitrate dehydrogenase; mdh, malate dehydronase; 
ode, OAA decarboxylase; pdh, pyruvate dehydrogenase com- 
plex; ppek, PEP carboxylase; pps, PEP synthetase; pyc, 
pyruvate carboxylase; abbreviations for other genes in Table 3. 
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estimate critical branch points and carbon fluxes 
during amino acid production by C. glutamicwn 
(Kiss and Stephanopoulos, 1 992; Stephanopoulos 
and Valino, 1991; Valino and Stephanopoulos, 
1993). This flux analysis suggested that the regu- 
lation of carbon flow at the phosphoenol pyruvate 
(PEP) branch point is the limiting factor in lysine 
production. At this point, three different enzyme 
systems compete for the available PEP, pyruvate 
kinase, PEP: sugar phosphotransferase system 
(PTS), and PEP carboxylase, as outlined in Fig- 
ure 1. To optimize, for example, lysine produc- 
tion, it seemed beneficial to amplify the PEP 
carboxylase {ppc) gene in a strain with a dis- 
rupted pyruvate kinase pyk gene. In such a strain, 
nearly all the carbon could be diverted by PEP 
carboxylase toward the synthesis of oxaloac- 
etate (OA A) and aspartate, and little or no car- 
bon would flow to the TCA cycle. In addition to 
diversion of the carbon flow to OAA, it would 
be necessary to inhibit the action of enzymes 
such as citrate synthase or OAA decarboxylase, 
which convert OAA into undesirable products. 
In order to direct even more carbon toward the 
synthesis of OAA, it also seemed beneficial to 
explore the presence of other enzyme systems 
such as the glyoxylate bypass or pyruvate car- 
boxylase in C. glutamicum. Recent studies con- 
cerning several of the aforementioned aspects of 
the intermediate metabolism of C. glutamicum 
are summarized below. 

B." Pyruvate Kinase 

1. Protein Purification 

Pyruvate kinase converts PEP with ADP, 
irreversibly to pyruvate and ATP, thereby di- 
recting the available carbon toward the TCA 
cycle and energy generation. Using a combina- 
tion of anion exchange and affinity chromatog- 
raphy, pyruvate kinase from C glutamicum was 
recently purified to homogeneity (Jetten et al., 
1994c). The purified enzyme (M r : 232 kDa) is 
composed of four subunits with a M r of 59 kDa 
and has a very high k^ (2520 s -1 ) compared to 
other bacterial enzymes. The enzyme was not 
influenced by fructose- 1 ,6-diphosphate and only 



used Mn 2+ or Co 2+ as cations. In addition to ADP 
(K m = 0.07 mA/), the purified enzyme also used 
GDP (K m = 0.25 mM) and IDP (K^ = 0. 12 nW) 
as substrates. The purified pyruvate kinase of 
C. glutamicum demonstrated complex sigmoidal 
kinetics with PEP as substrate. In the presence of 
the negative allosteric effector ATP, the concen- 
tration at half- maximal rate (S 05 ) increased from 
1.2 to 2.8 mM and cooperation as expressed by 
the Hill coefficient increased from 2.0 to 3.2. 
AMP promoted the opposite effects: the S 05 was 
decreased to 0.4 mM and the enzyme exhibited 
almost no cooperation. These complex kinetics 
show that C. glutamicum is able to regulate pyru- 
vate kinase activity to a certain extent. At high 
ATP levels, the enzyme is inhibited, thereby 
permitting the carbon flow to be directed toward 
biosynthesis, whereas under energy limitation, 
pyruvate kinase is maximally active and directs 
much of the carbon toward the TCA cycle and 
energy generation. 

2. Cloning Procedure 

As mentioned earlier, Gubler et al. (1994b) 
used a very elegant combination of several mo- 
lecular genetic techniques to clone the pyruvate 
kinase gene of C. glutamicum. On the basis of 
several published deduced amino acids se- 
quences, two degenerate oligonucleotides were 
designed and used to synthesize a 126-bp inter- 
nal fragment of the C. glutamicum pyruvate ki- 
nase gene (pyk) via a polymerase chain reaction. 
This internal fragment was cloned in the 
mobilizable pSUP301 vector and used to gener- 
ate defined Pyk* mutants of C. glutamicum via 
transconjugation. Phenotypical analysis of the 
Pyk" mutants revealed that they were no longer 
able to grow on minimal medium containing the 
non-PTS sugar ribose as a sole carbon and en- 
ergy source. The pyk gene could therefore be 
cloned on a 2.8-kb BamHl-Pstl fragment via 
complementation of the mutant on minimal ri- 
bose medium. Sequence analysis of the 2.8-kb 
fragment revealed an open reading frame of 1425 
nucleotides encoding a 475-amino acid residue 
polypeptide with a subunit M r of 59,912 Da. 
Sequence comparison with other bacterial pyru- 
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vate kinases showed an overall identity of 39 
to 44%, with some very conserved regions. The 
information obtained from sequence analysis 
was used to construct more efficient mobilizable 
plasmids for the disruption of the pyk gene via 
marker exchange mutagenesis in lysine pro- 
duction strains. 



3. Performance of Defined 
Pyk~ Mutants 

Both Pyk + and Pyk" strains of C. lactofer- 
mentum ATCC 21799 were used in shake-flask 
fermentations to determine the effect of the pyk 
mutation on lysine production. The Pyk + parent 
strain produced 26 g/1 of lysine together with 3 g/1 
of alanine, 5 g/1 of valine, and 6 g/1 of acetate in 
72 h from 90 g/1 of glucose. The Pyk" mutants 
only produced 15 g/1 of lysine together with 3 g/1 
of valine and 3 g/1 of alanine from 90 g/1 of 
glucose in 72 h.'ln addition to 13 g/1 of acetate, 
the Pyk" strains also excreted 3 g/1 of dihydroxy- 
acetone and 4 g/1 of glyceraldehyde, which might 
be indicative of an accumulation of glycolytic 
intermediates upstream of PEP. This surprising 
finding is in contrast to earlier reports by Mori 
and Shiio (1987), who claimed that the Pyk" 
mutation was beneficial for lysine production in 
C.flavum. This and the following examples will 
clearly show that there is yet much to be learned 
about the regulation of carbon flow in this group 
of organisms. 

C. Phosphoenol Pyruvate: Sugar 
Phosphotransferase (PTS) System 

To date, little information exists about the 
functionality of the PTS system in C. glutamicum. 
Mori and Shiio (1987a, b) published two papers 
on the presence of a PTS in C. glutamicum, which 
was confirmed by Malin and Bourd (1991). Re- 
cently, the cloning of a PTS component was 
reported (Yoon et al., 1993). The unclear 
phenotypical behavior of previously described 
Pyk" mutants indicates that the PTS system might 
not be the only way in which C. glutamicum is 
able to take up glucose. 



D. Phosphoenol Pyruvate Carboxylase 

PEP carboxylase converts PEP with C0 2 to 
oxaloacetate in order to replenish the TCA cycle 
with OAA, which is constantly removed from 
the cycle for the production of amino acids de- 
rived from aspartate. Because of this anaplerotic 
function, PEP carboxylase received considerable 
attention in the literature. The enzyme was puri- 
fied and characterized by Mori and Shiio (1985), 
who showed that the PEP carboxylase of 
C.flavum was stimulated by acetyl-coenzyme A 
and inhibited by aspartate and a-ketoglutarate. 
The ppc gene encoding PEP carboxylase was 
cloned and sequenced later on (Eikmanns et al., 
1989; O'Regan et al., 1989). Amplification of 
the ppc gene in lysine-producing C. glutamicum, 
however, did not result in higher yields. (Sano et 
al., 1987; Cremer et al., 1991). To elucidate the 
role of PEP carboxylase together pyruvate ki- 
nase in more detail, defined ppc and ppc/pyfr 
double mutants were constructed via marker 
exchange mutagenesis (Gubler et al., 1 994; Jetten, 
unpublished results). Independently ppc mutants 
were obtained by two research groups using 
slightly different approaches (Gubler et a!., 
1994a; Peterschwinde et al., 1993). In both cases, 
these ppc~ mutants did not show reduced growth 
rates and lysine production from glucose was 
indistinguishable from the parent strains (Gubler 
et al., 1994a; Peterschwinde et al., 1993). This 
surprising finding is in contrast to the general 
opinion that ppc is the exclusive anaplerotic 
pathway in C glutamicum (Mori and Shiio, 1985; 
Shiio et al., 1993), and indicates that other PEP- 
converting enzymes or even different anaplerotic 
reactions may play an important role in the regu- 
lation of carbon flux in C. glutamicum. 



E. Phosphoenol Pyruvate 
Carboxykinase 

One alternative way for the conversion of 
PEP was identified in both C. glutamicum ppc 
mutants and characterized in more detail (Gubler 
et al., 1994a; Peterschwinde et al., 1993; Jetten 
and Sinskey, 1993). The enzyme PEP carboxy- 
kinase was partially purified using anion exchange 
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chromatography to obtain PEP carboxykinase 
preparations free from interfering OAA-decar- 
boxylating activities. The PEP carboxykinase 
activity in these preparations was completely de- 
pendent on the presence of IDP (K^ 0.5 mAf) 
with a V max of 0.82 u.mol/min/mg in the direction 
of OAA formation. ADP could not substitute for 
IDP, which is in contrast to the situation observed 
in other bacterial PEP carboxy kinases. The find- 
ing that ATP inhibited 60% of the enzyme activ- 
ity at 0.1 mM in the OAA-forming reaction indi- 
cated that the enzyme probably only functions in 
the gluconeogenetic direction and might be regu- 
lated by the adenosine nucleotide pools. Once the 
pck gene encoding for PEP carboxykinase be- 
comes available and defined pckr and pp<r!pclc 
mutants can be constructed, more concrete evi- 
dence about the role of PEP carboxykinase in 
C. glutamicum can be obtained. 



F. Does C. glutamicum Have 
Pyruvate Carboxylase? 

Other alternative pathways for the genera- 
tion of OAA besides PEP carboxykinase could 
be provided by pyruvate carboxylase, which 
converts pyruvate with C0 2 and ATP to OAA. 
To our knowledge, only one report exists in 
which the authors have claimed to be able to 
detect the presence of pyruvate carboxylase in 
C. lactofermentum (Tosaka et ah, 1979). Sev- 
eral groups have tried to confirm this finding, 
but were unable to do so (Dominguez et al., 
1993; Jetten et al., 1994a; Neubacket al., 1993). 
More research, especially with defined double or 
triple mutants (ppc/pykr/pckr) of C. glutamicum 
or C. lactofermentum might be necessary to either 
rule out the existence of a pyruvate carboxylase 
or confirm its presence. Preliminary studies with 
ppcipykr double mutants of strain 21799 
showed that the mutants were no longer able to 
grow on the non-PTS sugars gluconate and ri- 
bose unless provided with pyruvate (Jetten, un- 
published results). The production of lysine from 
pyruvate or acetate in these mutants did not 
seem to be impaired compared with the parent 
strain, which could indicate the presence of 
pyruvate carboxylase. 



G. Glyoxylat Bypass 

The glyoxylate bypass, in which isocitrate 
with one molecule of acetyl-coenzyme A is c ™ 
verted to malate and succinate via glyoxylate by 
the consecutive action of isocitrate lyase and 
malate synthase, could also function as an 
anaplerotic route in C glutamicum. Induction and 
repression of the genes encoding ICL and MS 
(aceA and aceB) have not been studied in 
C. glutamicum (Ozaki and Shiio, 1969). When 
the organism is cultured on C2 (acetate) or C3 
(pyruvate/lactate) compounds, the bypass seems 
to be fully induced, but regulation mechanisms of 
the isocitrate branch-point such as the isocitrate 
dehydrogenase kinase/phosphatase (e.g., in E. coli) 
have not been described (Ozaki and Shiio, 1969). 
Some reports indicate that the aceA-aceB operon 
is not completely repressed in the presence of 
excess glucose, as has been observed in E. coli 
(Gubler et al., 1994a; Ozaki and Shiio, 1969; 
Schrumpf et al., 1992; Reinscheid et al., 1994b). 
This could be due to either a different regulation 
mechanism or the induction of the bypass by 
acetate, which accumulates in the medium as a 
fermentation product. Recently, the genes encod- 
ing aceA and aceB have been cloned from 
C. glutamicum (Reinscheid et al., 1993; Lee et 
al., 1994). Preliminary results showed that the 
genes are organized differently in C glutamicum 
than in E. coli and that the aceB gene of 
C. glutamicum is considerably larger than in 
E. coli. The construction of defined aceA' mu- 
tants in combination with other mutations {ppcf 
pykr) will provide more information about carbon 
flow via the glyoxylate bypass and its contribu- 
tion to OAA synthesis in C. glutamicum. 



H. Oxaloacetate Decarboxylase 

One other important aspect of central carbon 
metabolism in C glutamicum is the often reported 
high OAA decarboxylating activities in crude 
extracts of this organism. Using various chro- 
matographic procedures, the OAA -decarboxylat- 
ing enzyme was purified 300-fold from a 
C. glutamicum pyruvate kinase mutant (Jetten and 
Sinskey, 1995). The molecular mass determined 
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by gel filtration was 118 kDa, while SDS-PAGE 
only showed one subunit of 32 kDa, indicating an 
a 4 subunit structure for the native enzyme. The 
enzyme exclusively decarboxylated OAA and did 
not use any other a-keto acids as substrate. Mn 2+ 
was required for full activity but could be re- 
placed by Mg 2+ , Nt 2+ , Co 2+ , and Ca 2+ . The enzyme 
has a low affinity for OAA (K^ = 2.1 mM) and 
Mn 2+ (K„ = 1.2 mM) and of 158 nmol/min/ 
mg. The enzyme was not inhibited by avidin, 
indicating that biotin is not involved in the decar- 
boxylation. The purified enzyme did not show 
any pyruvate kinase, pyruvate carboxylase, malate 
dehydrogenase, PEP carboxykinase, or malic en- 
zyme activities, demonstrating that it belongs to 
the divalent cation-dependent OAA decarboxy- 
lases. In contrast to several other OAA decar- 
boxylases, the enzyme from C. glutamicum was 
not inhibited by acetyl-CoA, but greatly affected 
by Cu 2+ , ADP, GDP, and HSCoA. This inhibition 
pattern is not easily explained. One could specu- 
late that OAA-decarboxylating activity serves to 
prevent the build up of OAA in the cell. In the 
absence of a sufficient acetyl-coenzyme A sup- 
ply, a high OAA concentration may lead to an 
inhibition of the TCA cycle by shifting the OAA/ 
malate equilibrium. The decarboxylation of OAA 
to pyruvate, together with the oxidative conver- 
sion of the formed pyruvate to acetyl-SCoA, would 
provide enough precursors for citrate synthase to 
continue the operation of the TCA cycle. To what 
extent the OAA decarboxylase contributes to in 
vivo OAA decarboxylation is not known (Jetten 
et al., 1993; Jetten and Sinskey, 1995). 



I. Miscellaneous 

Recently, several more genes encoding en- 
zymes involved in central carbon metabolism 
have been isolated and characterized. Eikmanns 
(1992) located a large cluster of glycolytic genes 
upstream of the ppc gene by sequence analysis. 
The cluster consists of gap (glyceraldehyde-3- 
phosphate dehydrogenase), pgk (3-phosphoglyc- 
erate kinase), and tpi (triose-phosphate isomerase). 
Transcriptional analysis of the gene cluster by 
Northern blots identified four different mRNAs 
corresponding to gap, gap-pgk-tpi, pgk-tpi, and 



pgk-tpi-ppc. Primer extension analysis located 
the transcriptional initiation sites in front of gap 
and pgk. The DNA sequence of these promoter 
structures was aligned with previously deter- 
mined transcriptional start sites of C. glutamicum. 
No homology was found with typical E. coli or 

B. subtilis promoters, but a tentative consensus 
sequence at -35 was identified as TTGACA. 
Mutational analysis by McCormick et al. (1993) 
(see also, Section II.E), however, showed a 
TGTGC motif as the most conserved among 

C. glutamicum promoters. This discrepancy 
clearly points out that additional information 
based on DNA footprinting is necessary before 
more definite conclusions concerning structure- 
function relationships can be made. 

J. Conclusion on Precursor Supply 

Although there have been many papers de- 
scribing mutants defective in enzymes of central 
metabolism (Shiio et al., 1991), our understand- 
ing of the central metabolism in C. glutamicum is 
still very limited. Now that new, more powerful 
genetic tools and more genes involved in the cen- 
tral metabolism have become available, a new 
phase in the study of the intermediate metabolism 
can start. The construction of defined double or 
triple mutants will be necessary to provide strains 
with proper genetic backgrounds. These strains 
can be used to identify which enzyme systems 
contribute to what extent in the biosynthesis of 
various compounds. Similar to recent studies with 
E. coli, the effect of the overexpression of differ- 
ent enzymes {ppc, pck, and pyk) on growth rate 
and metabolite production in C. glutamicum mu- 
tants can be evaluated and used to improve amino 
acid production strains (Chao and Liao, 1993; 
Chao et al., 1993). 



IV. GLUTAMATE 

A. Enzymes Responsible for 
Glutamate Synthesis 

Under certain culture conditions such as bi- 
otin limitation or addition of surfactants, 
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C. glutamicum is able to excrete large amounts of 
L-glutamate into the medium (Clement and 
Laneelle, 1986; Kinoshita, 1985; Shiio et al., 
1963). The generally accepted explanation for 
this high-level accumulation is that the exit of 
glutamate leads to decreases in the internal pool 
of this amino acid, which in turn relieves feed- 
back control mechanisms (Liebl et ai., 1991). 
However, Frings et al. (1993) showed that 
glutamate and glutamine serve as osmoprotectants 
in Corynebacterium strains exposed to high salt 
concentrations. It has been further postulated that 
the low activity and high instability of the a- 
ketoglutarate dehydrogenase from C. glutamicum 
favors the conversion of a-ketoglutarate by 
glutamate dehydrogenase (GDH), especially in 
the presence of high NHJ concentrations (Shiio 
and Ujigawa, 1980). Glutamate can be formed by 
either the previously mentioned GDH or the 
coupled action of glutamine synthase (GS) and 
glutamate synthase (GOGAT) under low NHJ 
concentrations. All three enzymes have been de- 
tected in and purified from C. glutamicum strains 
(Takeda et al., 1990), but only recently has the 
gene encoding GDH been cloned from both 
C. glutamicum and C. melassecola (Bormann et 
al., 1991; Takeda et al., 1990). Sequence analysis 
of the C. glutamicum GDH reading frame showed 
that the gdh gene consists of 1344 bp encoding a 
448-amino acid residue polypeptide. The deduced 
amino acid sequence showed between 47 and 
54% identity with amino acid sequences of other 
organisms (Bormann et al., 1991). Physiological 
studies with overexpressed gdh revealed no indi- 
cations of significant regulation of gdh expres- 
sion. Northern blot analysis indicated a 1.7-kb 
monocistronic gdh transcript, while RNA poly- 
merase footprinting showed a transcription initia- 
tion site 185 nucleotides upstream from the ATG 
start codon. 

Glutamate synthesis in C. glutamicum is gen- 
erally believed to be dependent on the activity of 
GDH because (1) the GS/GOGAT system is re- 
pressed under the high NHJ concentrations of the 
fermentation media employed, (2) the GDH ac- 
tivity is orders of magnitude higher than the GS/ 
GOGAT system, and (3) GDH mutants of 
C.flavum were described as glutamate auxotro- 
phs with negligible glutamate excretion. In order 



to verify this hypothesis, Bormann et al. (1993) 
constructed defined GDH deficient mutants via 
marker exchange mutagenesis and strains with 
tenfold overexpressed GDH activity. Comparison 
of these two strains with the wild type revealed 
that GDH is dispensable for glutamate synthesis 
required for growth and that the level of GDH has 
no influence on glutamate accumulation. The au- 
thors related this surprising result to the compen- 
satory capacity of the GS/GOGAT system and a 
limitation of the glutamate secretion system. 

B. Glutamate Excretion 

The mechanism of glutamate transport into 
the medium has been subject of a long contro- 
versy in the literature. Three different mecha- 
nisms have been proposed to describe glutamate 
efflux. The "leak" model is based on the observa- 
tion that glutamate efflux can be induced by vari- 
ous treatments such as biotin limitation, addition 
of sublethal amounts of penicillin, or addition of 
surfactants, all of which correlate with alterations 
in the cell membrane (Clement and Laneelle, 1986; 
Duperray et al., 1992). The second hypothesis is 
based on a general model for many efflux pro- 
cesses in bacteria, namely, inversion of the up- 
take process induced by changes in chemical po- 
tential or regulation, or by uncoupling due to 
changes in the membrane structure (Clement et 
al., 1984). Kramer and co-workers (Gutmann et 
al., 1992; Hoischen and Kramer, 1989; Hoischen 
and Kramer, 1990; Plakunovetal., 1992; Kramer, 
1 994), however, have proposed a model based on 
the presence of a specific glutamate carrier be- 
cause several of the observations did not agree 
with either the "leak" or the "inversion" model: 
(1) C. glutamicum cells were able to export 
glutamate against an existing gradient, which 
makes it unlikely that the internal/external 
glutamate ratio can be the driving force for 
glutamate efflux (Gutmann et al., 1992), (2) 
glutamate export proceeded even in the absence 
of a membrane potential (Gutmann et al., 1992), 
(3) the transport of glutamate was not coupled 
with the movement of H + , K + , or Cl~ ions, as 
would be expected in the "leak" model (Gutmann 
et al., 1992), (4) the carrier was specific for 
glutamate; no transport of amino acids other than 
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glutamate was observed (Hoischen and Kramer, 
1989; Gutmann et al, 1992), (5) the observed 
difference between the excretion and uptake rates 
(Kramer and Lambert, 1991 ; Kramer et ah, 1991) 
ruled out the inversion of the uptake system, (6) 
no changes were observed in the fatty acid or 
phospholipid contents when cells were switched 
between nonproducing and producing states 
(Hoischen and Kramer, 1990), and (7) there was 
no specific membrane component that regulated 
glutamate efflux (Hoischen and Kramer, 1990). 
In addition, Neuback et al. (1993) could not de- 
tect any differences in the anisotropy of mem- 
branes from producing or nonproducing cells. 
Because no change in anisotropy of membranes 
was observed during any period in the fermenta- 
tion, Neuback et al. (1993) concluded that 
glutamate secretion was not coupled to differ- 
ences in membrane fluidity, making the efflux of 
glutamate via a "leak" very unlikely. 

Gutmann et al. (1992) determined that there 
was a positive correlation between the secretion 
rate and the intracellular ATP pool, indicating that 
ATP or a high-energy compound could be involved 
in the activity of the secretion system. Isolation and 
characterization of the specific glutamate carrier of 
C. glutamicum and cloning and characterization of 
the gene(s) encoding the glutamate transport sys- 
tem will be necessary to provide further evidence 
for the active transport model. 

V. LYSINE 

A. Aspartokinase 

Lysine biosynthesis in amino acid-producing 
Corynebacterium strains proceeds via the 
diaminopimelate (DAP) pathway (Figure 2). 
Aspartokinase (ASK) and aspartate semialdehyde 
dehydrogenase (ASD) catalyze the initial reac- 
tions of this pathway (Follettie et al., 1993; 
Kalinowski et al., 1991; Tosaka and Takinami, 
1978). Regulation of the multibranched biosyn- 
thetic pathway of the aspartate family of amino 
acids in C. glutamicum is relatively simple. In 
contrast to E. coli and Bacillus sp., which have 
three distinct isoenzymes, only one ASK has been 
observed in C. glutamicum so far (Follettie et al., 
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1993; Kalinowski et al., 1991). However, studies 
with a defined ask mutant of C. lactofermentum y 
obtained via marker exchange mutagenesis, 
showed that the ask mutant was still able to grow 
on minimal medium, indicating that more than 
one ASK could be present in this organism (Jetten 
et al., 1995). Also, there is no apparent regulation 
of the synthesis of any enzymes of the lysine 
pathway in C. glutamicum. Instead, this pathway 
is thought to be coordinated by relative specific 
activities, and patterns of inhibition and repres- 
sion at key branch points (Follettie et al., 1993; 
Tosaka and Takinami, 1978). The main control of 
carbon flow into this pathway is mediated by 
concerted feedback inhibition of ASK by lysine 
and threonine (Shiio and Miyajima, 1969). Wild- 
type C. glutamicum does not excrete significant 
amounts of lysine (Schrumpf et al., 1992). Clas- 
sically, lysine-producing C. glutamicum strains 
were obtained by screening for resistance to the 
lysine analog S-(cc-aminoethyl)-D,L-cysteine 
( AEC) and by screening for homoserine and threo- 
nine auxotrophs (Tosaka and Takinami, 1978). 
Most of these resistant strains contain an ASK 
that is resistant to feedback inhibition by lysine 
and threonine (Shiio et al., 1991, 1993). 

In order to elucidate the mechanism of feed- 
back resistance of ASK in C. glutamicum, the 
ask genes from both wild-type and feedback re- 
sistant, lysine-producing C. glutamicum strains 
have been isolated and characterized (Kalinowski 
etal., 1990, 1991; Follettie et al., 1993; Jetten et 
al., 1995). The ask and asd genes were found to 
be organized in an operon and separated by a 
23-base pair intercistronic region. The ask-asd 
operon can be isolated by complementation of 
E. coli strain 5080, which lacks an asd gene 
(Follettie et al., 1993; Jetten et al., 1995). This 
strain therefore cannot grow on LB medium, 
which does not contain the diaminopimelate 
(DAP) needed for cell wall synthesis. The com- 
plete nucleotide sequence of the ask-asd opera- 
tion has been determined for both C glutamicum 
and Cflavum strains (Kalinowski et al., 1991; 
Follettie et al., 1993). Computer analysis of po- 
tential protein-coding regions identifies two long 
open-reading frames of 421 and 344 residue pro- 
teins, respectively. Analysis of the amino acid 
sequence of the first open-reading frame showed 
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significant homology and identity to previously 
characterized aspartokinases in other microor- 
ganisms. Deletion analysis and SI nuclease 
mapping showed that two promoters are required 
for expression of the complete operon. Promoter 
PI is located 38 nucleotides upstream of the 
predicted translation initiation site of the ask 
gene, while a second promoter (P2) responsible 
for 40% of asd expression was found to be lo- 
cated inside the ask gene. The presence of this 
second promoter, a potential ribosome-binding 
site at position 1 105 and a GTG translation ini- 
tiation codon at position 1 116, provided excel- 
lent evidence that the C. glutamicum ask gene 
encoded a second polypeptide. The expression 
of the postulated AKa and AK(3 subunits was 
confirmed by Western blot analysis with 
polyclonal antibodies generated to the C-termi- 
nal 179-amino acid residues of the C. glutamicum 
ask gene product. This organization is similar to 



the askll gene in B. subtilis that encodes the Gl- 
and ^-subunits required for functional asparto- 
kinase enzyme activity. 

In order to determine the molecular basis of 
deregulation of ASK and concomitant AEC resis- 
tance, the nucleotide sequences of wild-type and 
feedback-resistant ask genes were compared 
(Kalinowski et al., 1990, 1991; Follettie et al., 
1993; Jetten et al., 1995). Only four of the ob- 
served mutations in the DNA sequence resulted 
in changes in the amino acid sequence of the 
feedback-resistant ASK. All four mutations are 
located in the coding region for the ask P-subunit. 
Analysis of hybrid ask genes containing one or the 
other of the amino acid differences demonstrated 
that the deregulation of the aspartokinases were 
due to an S301Y, G345D, or R349G substitution 
(Kalinowski et al., 1990; Follettie et al., 1993; 
Jetten et al., 1994). Because all three deregulated 
aspartokinases characterized so far contain muta- 
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tions in the carboxy terminus of the enzyme, a 
regulatory role of the P-subunit has been proposed 
by Follettie et al. (1993). More detailed analysis 
and site-directed mutagenesis are necessary to lo- 
cate the different binding sites of threonine and 
lysine to the aspartokinase and to understand the 
mechanism of the feedback inhibition. 

The effect of deregulated, feedback-resistant 
ASK on lysine production has been studied in 
various high lysine-producing C. glutamicum 
strains (Cremeretal., 1990; Schrumpf etal., 1992; 
Shiio et al., 1991; Jetten et al., 1995). In most 
cases, amplification of the ask gene resulted in 
significant increases in lysine production, indicat- 
ing that in these strains ASK levels are still the 
limiting factor. 



B. Enzymes Involved in Synthesis and 
Conversion of Mesodiaminopimelate 

The pathway of aspartate-derived amino ac- 
ids diverges after the synthesis of aspartate 
semialdehyde, one branch leading to the biosyn- 
thesis of DAP and lysine, the other one to threo- 
nine, methionine, and isoleucine. The condensa- 
tion of ASA with pyruvate to dihydropicolinate 
(DHP) is the committed step in the conversion of 
ASA via DAP to lysine. This reaction is cata- 
lyzed by DHP synthase. The enzyme has a rather 
low affinity for both ASA and pyruvate (Cremer 
et al., 1988), compared to the enzyme competing 
fqr ASA, homoserine dehydrogenase (Archer et 
al., 1991; Follettie et ah, 1988). The gene for 
DHPS, daph, has been cloned by complementa- 
tion of E. coli dap A mutants and was found to be 
closely linked to dapB, the gene encoding for 
the next step in the pathway, DHP reductase 
(Cremer et al., 1988; Pisabarro et al., 1993; Yeh 
et al., 1988). Neither the activity nor the expres- 
sion of the DHPS seems to be affected by lysine 
(Cremer et al., 1988; Yeh et al., 1988). Shiio et 
al. (1987) reported that mutants with defective 
DHPS were reasonable threonine producers, but 
required high levels of meso-DAP for optimal 
growth. Although no meso-DAP transport is 
known to occur in C. glutamicum, the overpro- 
duction of threonine by these mutants would be 
consistent with the notion that the competition 



of DHPS and homoserine dehydrogenase for 
ASA regulates the flow of intermediates into the 
two main branches of this pathway. 

The dapB gene encoding the DHP reductase 
of C. glutamicum has been cloned by comple- 
mentation of E. coli dapB mutants and is closely 
linked to dap A (Cremer et al., 1988; Pisabarro et 
al., 1993; Yeh etal., l988).The^pA-^pB genes 
are clustered on a 3.5-kb Sau3A\-BamH\ frag- 
ment, but separated by an open reading frame of 
750 bp. Pisabarro et al. ( 1 993) speculated that this 
open reading frame may encode for a 29-kDa 
polypeptide involved in the cyclization of 
dihydropicolinate. 

C. glutamicum has a high level of meso-DAP 
dehydrogenase and contains all the enzymes of 
the epimerase pathway for meso-DAP synthesis 
(Schrumpf et al., 1991; Yeh et al., 1988). The 
epimerase pathway in C. glutamicum seems to 
prefer succinylated intermediates. In C.flavum, 
however, higher activities have been reported for 
the acetylated intermediates (Tosaka and 
Takinami, 1978). The gene for meso-DAP dehy- 
drogenase {ddh) has been cloned by its ability to 
restore lysine synthesis in E. coli dapD mutants 
with a defective epimerase pathway (Ishino et al., 
1 987, 1 988; Yeh et al., 1 988). Studies with strains 
containing either a disrupted ddh gene or a dis- 
rupted dapE gene showed that both the epimerase 
and dehydrogenase pathways contribute to lysine 
biosynthesis (Schrumpf et al., 1991; Sonntag et 
al., 1993; Wehrman et al., 1993). Using 13 C la- 
beled glucose, Sonntag et al. (1993) determined 
how the flux was distributed between the differ- 
ent pathways. They showed that the ddh pathway 
contributed about 30% to lysine synthesis in both 
wild-type and lysine production strains. Pseudo 
steady-state experiments also revealed that the 
flux partitioning over the ddh pathway in a pro- 
duction strain decreased from 70% in the initial 
stages of the fermentations to 0% at the end of 
lysine accumulation, probably due to NH} deple- 
tion. Although lysine itself had no effect on the 
activity of the DDH, the enzyme did not contrib- 
ute to lysine synthesis when NHJ was replaced by 
glutamate as an ammonium source, showing that 
this anabolic pathway is directly influenced by an 
extracellular medium component (Sonntag et al., 
1993; Cremer et al., 1988). 
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The final step in the biosynthesis of lysine is 
the pyridoxal phosphate-dependent decarboxy- 
lation of meso-DAP to lysine, catalyzed by DAP 
decarboxylase. The gene encoding DAP decar- 
boxylase (lysA) has considerable homology to 
the gene of E. coli (Yeh et al., 1988). The ex- 
pression of lysA in C glutamicum depends on 
the presence of a 2.35-kb region upstream of the 
lysA coding region (Marcel et al., 1989). On the 
basis of DNA homology, this region was re- 
cently identified as a possible coding sequence 
for arginyl-tRNA synthetase (Oguiza et al., 1993; 
Sharp and Mitchell, 1993). ArgS appears to be 
cotranscribed with the lysA gene in a single 
3000-nucIeotide transcript (Oguiza et al., 1993). 
Oguiza et al. showed that the transcription of the 
argS-lysA operon was repressed by lysine and 
induced by arginine, indicating a peculiar pat- 
tern of interlock between these two pathways in 
C. glutamicum. 

C. Lysine Transport 

Wild-type Corynebacterium strains do not 
secrete significant amounts of lysine into the 
medium (Schrumpf et al., 1992). Effective high 
lysine secretion is only possible when (1) con- 
certed feedback inhibition of lysine and threo- 
nine on aspartokinase is removed and (2) addi- 
tional mutations in the lysine export system 
have occurred. High internal lysine concentra- 
tions and concomitant lysine secretion are 
thought to be the result of these mutations (Broer 
and Kramer, 1991, 1992a, b). The lysine secre- 
tion systems of three different lysine-produc- 
ing C. glutamicum strains were characterized 
and compared to verify this hypothesis (Broer 
et al., 1993). A pseudo wild-type strain (KK- 
25) with a resistant aspartokinase gene intro- 
duced into the chromosome via conjugal trans- 
fer showed only a very low lysine export rate of 
2 nmol/min/mg dry weight, and the strain did 
not secrete lysine when the internal concentra- 
tion was below 1 5 mA/. C. glutamicum produc- 
tion strains DG 52-5 and MH20-22B, obtained 
via classical mutagenesis and screening, had 
comparable high lysine secretion rates (1 1 nmol/ 
min/mg dry weight) and K m values (1 8 ± 5 mM)> 
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but showed subtle differences in dependence 
on ApH and membrane potential. On the basis 
of these observations, it was proposed that at 
least strain MH20-22B carried a mutated form 
of the lysine export permease. The same con- 
clusion was reached in another study, where 
strain MH20-22B was characterized in more 
detail (Schrumpf et al., 1993). The strain was 
obtained after two rounds of mutagenesis and 
screening, the first one based on AEC resis- 
tance, the second on leucine auxotrophy and 
lysine production. The main differences with its 
predecessor are feedback-resistant aspartokinase, 
lack of isopropylmalate dehydratase, and a high 
secretion rate of lysine (10 nmol/min/mg dry 
weight). Application of site-directed mutagen- 
esis to this strain revealed a key role for de- 
regulated aspartokinase and the probable pres- 
ence of a mutated secretion system, as observed 
by Broer and Kramer (1993). 

To study the wild-type lysine secretion sys- 
tem of C. glutamicum in more detail, Erdmann 
et al. (1993) investigated the effect of dipeptide 
uptake on lysine excretion. Uptake of dipeptides 
containing lysine increased the internal lysine 
concentration, thereby triggering the secretion 
process. The secretion of lysine depended on the 
medium composition in which the cells were 
grown. Experiments where protein synthesis was 
inhibited with chloramphenicol indicated thatde 
novo protein synthesis is necessary for the regu- 
lation of lysine secretion. 

Because C. glutamicum mutants can produce 
up to 1 M of the cation lysine in the medium, 
export must not only proceed against a chemical 
gradient but also against the membrane poten- 
tial. As the ATP content of the cells did not 
correlate with the export activity, it was pro- 
posed that lysine is presumably secreted in a 
symport mechanism with two OH", which is 
energetically equivalent to a lysine/2H + antiport 
mechanism. In this way, the positive charge of 
lysine is compensated for and the membrane 
potential can function as the driving force for 
excretion. Similar to elucidation of the glutamate 
transport in C. glutamicum, the availability of 
the gene encoding the lysine permease will 
greatly contribute to our current understanding 
of the lysine export process. 
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VI. THREONINE 

A. Molecular Organization of the 
Threonin Biosynthetic Pathway 

As already mentioned, the relative affinities 
of the competing enzymes homoserine dehydro- 
genase and dihydropicolinate synthase determine 
the distribution of the common substrate aspar- 
tate-|3-semialdehyde (ASA) to the divergent threo- 
nine and lysine biosynthesis pathways (Follettie 
et al., 1993). Carbon flows preferentially toward 
the threonine pathway at the expense of lysine 
biosynthesis as a result of the 1 5-fold-higher sub- 
strate affinity and substrate conversion rate ex- 
hibited by homoserine dehydrogenase over its 
competing enzyme dihydropicolinate synthase 
(Shiio et al., 1989). The activity of homoserine 
dehydrogenase is highly sensitive to allosteric 
inhibition by L-threonine (K { = 0.16 mAf) and 
therefore, in the presence of sufficient threonine, 
ASA predominantly enters the lysine pathway 
(Figure 2). 

The C. glutamicum genes encoding the first 
two enzymes in the threonine pathway, homoserine 
dehydrogenase and homoserine kinase, were iso- 
lated by complementation of E. coli thrB mutant 
(Follettie et al., 1988). The two genes form an 
operon that is expressed from a single promoter 
located immediately upstream of the horn gene 
(Peoples et al., 1988). The hom-thrB operon is 
genetically regulated at the level of transcription 
by a threefold repression by methionine via a 
unique attenuator system (Jetten et al., 1993). 
Two research groups have isolated and character- 
ized genes encoding deregulated feedback-resis- 
tant homoserine dehydrogenase (Archer et al., 
1992; Reinscheid et al., 1991) in order to investi- 
gate the molecular basis for L-threonine allosteric 
inhibition of homoserine dehydrogenase. The 
complete nucleotide sequence of both horr& genes 
was determined, and comparison with the wild- 
type sequence revealed that both mutations re- 
sponsible for the deregulation are located in the 
carboxyl end of the enzyme. Archer et al. (1992) 
determined that the horn* 1 mutation is due to a 
single base deletion in codon 429 within the 
homoserine dehydrogenase gene. The resulting 
frameshift mutation radically alters the structure 



of the carboxy terminus, leading to ten amino 
acid changes and a deletion of the last seven 
residues relative to the wild type. Similar results 
were obtained by Reinscheid et al. (1992), who 
studied the deregulated hom iT gene of the AHV- 
resistant C glutamicum mutant DM 368-3. A 
DNA fragment exchange experiment between 
wild-type and deregulated homoserine dehydro- 
genase genes showed that a 0.23-kb region close 
to the 3' terminus of the hom Ar gene was respon- 
sible for the deregulation. Comparison of se- 
quence data revealed that the mutation in the 
hom dr gene was due to single transition from G 
to A at position 1 133, leading to a replacement 
of glycine-378 by glutamate. 

The final step in the biosynthesis of threonine 
is the conversion of homoserine- phosphate to threo- 
nine. This reaction is catalyzed by the enzyme 
threonine synthase. The thrC gene encoding threo- 
nine synthase was obtained by complementation of 
a C. glutamicum threonine auxotroph (Han et al., 
1989). The thrC gene is expressed as a single 
1 .3-kb transcript that is translated into a polypep- 
tide of 54.5 kDa. Expression of threonine synthase 
appeared to be constitutive in C. glutamicum. The 
predicted amino acid sequence of threonine synthase 
showed homologies of 22 and 15% with the pre- 
dicted thrC gene products of E. coli and B. subtilis, 
respectively. The homology was predominantly 
located in the N-termiiial part of the proteins, and 
it has been suggested that this region may bind the 
pyridoxal-phosphate coenzyme. 

B. Threonine Production by 
Recombinant C. glutamicum Strains 

In order to obtain information on the limit- 
ing steps in carbon flux to threonine and con- 
struct well-defined recombinant threonine-pro- 
ducing strains, the effects of the combined 
amplification of the threonine biosynthetic genes 
on amino acid production in wild-type and lysine 
plus threonine-producing strains has been stud- 
ied (Eikmanns et al., 1991; Reinscheid et al., 
1994; Colon et al., 1993). Amplification of 
threonine genes in wild-type C. glutamicum 
13032 did not result in any threonine secretion, 
reflecting the key role of the still feedback- 
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inhibited aspartokinase of this strain (Eikmanns 
et al., 1991). Introduction of the deregulated 
homoserine dehydrogenase gene hom dl together 
with the wild-type homoserine kinase gene thrB 
on a multicopy plasmid in lysine-producing 
strains resulted in unstable transformant that 
produced only limited amounts of threonine 
(Archer et al., 1992; Eikmanns et al., 1991). In 
order to solve the instability problems, 
Reinscheid et al. (1994) integrated the hom 6s - 
thrB operon into the chromosome of the lysine 
producer MH20-22B via conjugal transfer. In 
this way, they obtained strains with one, two, 
or three copies of the deregulated hom dl -thrB 
operon that were stable for more than 70 gen- 
erations. In the MH20-22B strains with one 
copy of the hom dr -thrB operon, lysine produc- 
tion was reduced from 208 to 129 mM as com- 
pared to the parent strain, with concomitant 
production of 39 mM threonine, 14 mM 
homoserine, 21 mM isoleucine, and 33 mM 
glycine. Strains containing three copies of the 
deregulated operon produced 58 mM lysine, 
65 mM threonine, 21 mM homoserine, 31 mM 
isoleucine, and 30 mM glycine. The authors 
further observed a drastic increase in the inter- 
nal concentration of both threonine (100 mM) 
and homoserine (74 mM). It was therefore con- 
cluded that overall threonine production is most 
probably limited by the efflux of threonine 
(Reinscheid et al., 1994). Similar to the results 
of Reinscheid et al. ( 1 994), Colon et al. ( 1 993) 
observed increased levels of homoserine and 
isoleucine in cells of C. lactofermentum with 
an amplified hom dT -thrB operon. The accumu- 
lation of homoserine is most likely caused by 
the inhibition of the thrB product homoserine 
kinase at higher threonine concentrations 
(Reinscheid et al., 1994). In order to prevent an 
intracellular build up of homoserine, Colon et 
al. (1993) fused the thrB gene to the tac pro- 
moter and regulated hom dT /P tac- thrB expres- 
sion by the addition of IPTG. In this way, com- 
plete conversion of homoserine to threonine 
was observed. Unfortunately, most of the threo- 
nine produced was either converted to isoleu- 
cine or degraded to glycine, resulting in only 
moderately increased threonine production. The 
conversion of threonine to isoleucine was pre- 
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vented by the construction of defined ilvA 
mutants via marker exchange mutagenesis 
(Jetten, unpublished results). The degradation 
of threonine to glycine is not easily solved. At 
least two different enzyme systems have been 
reported to convert threonine to glycine. The 
first system consists of the combined action of 
threonine dehydrogenase and 2-amino-3- 
oxobutyrate-CoA ligase, which converts threo- 
nine via 2-amino-3-oxobutyrate to glycine and 
acetyl-CoA (Bell and Turner, 1976). The second 
possible mechanism is nonspecific conversion 
by the glyA gene product serine hydroxymethyl- 
transferase (Plamann and Stauffer, 1983; Schirch 
etal., 1985). Unfortunately, neither of the genes 
for these systems have yet been characterized 
in C. glutamicum, and hence defined mutants 
obtained via marker exchange mutagenesis will 
not be available soon. The fastest way of ob- 
taining mutants defective in either glyA or tdh 
would be screening for strains unable to grow 
on threonine or glycine as the sole carbon and 
nitrogen source. 

Another approach to develop C. glutamicum 
threonine production strains has been amplifica- 
tion of the thr operon from E. colt in C. flavum 
(Ishida et al., 1993; Patek et al., 1989, 1993). 
Plasmids containing the E. coli thr operon were 
only maintained by applying high trimethoprim 
concentrations in the broth. In this way, Ishida et 
al. (1993) obtained a C. flavum strain that pro- 
duced 64.4 g/1 of threonine from acetate (700 g/1) 
and glucose (100 g/1) with a yield of 19.9% in 
92 h. The strains used by Patek et al. (1993) 
produced 12 g/1 of threonine from 180 g/1 of 
sucrose in 48 h, which is approximately half the 
amount obtained by Reinscheid et al. (1994). 

Many of the problems experienced in the 
construction of C. glutamicum threonine produc- 
tion (redirection of carbon flow, threonine degra- 
dation) have also been encountered in the breed- 
ing of either E. coli (Debadov, 1982) or 
S. marcescens threonine production strains 
(Masuda et al., 1993). Recently, Masuda et al. 
(1993) described the performance of S. marcescens 
strain T-2000 in batch fermentation. Strain T-2000 
was obtained by transformation of strain T- 1 165 
with plasmid pSK301 containing a mutant threo- 
nine operon. Strain T-l 165 itself is defective in 
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two threonine-degrading enzymes and already con- 
tains the deregulated threonine operon on the 
chromosome. Introduction of the low-copy-num- 
ber plasmid pSK301 did not give the expected 
increase in threonine yield. Therefore, the effect 
of different nitrogen sources on threonine produc- 
tion by strain T-200 was investigated. By chang- 
ing the nitrogen source from the normally em- 
ployed urea to ammonium, the authors were able 
to obtain an increase in threonine titers from 55 to 
100 g/1 from 300 g/1 of sucrose, showing that 
maintaining sufficient levels of ammonium in the 
fermentation broth can be very critical. Without 
mentioning any specific bottlenecks, the authors 
even expected greater improvements in threonine 
yield if they could optimize the precursor flow 
toward the threonine pathway. 

VII. 1SOLEUCINE 

A. Molecular Organization of Genes 
Involved in Isoleucine Biosynthesis 

The biosynthesis of isoleucine in C. gluta- 
micum and other organisms starts with the con- 
version of L-threonine to 2-oxobutyrate (Fig- 
ure 2) catalyzed by the *7vA gene product 
L-threonine deaminase (LTD) (Mockel et al., 
1992). The next step in the synthesis is per- 
formed by an acetohydroxy acid synthase 
(AHAS), which naturally uses both 2-oxobutyrate 
and 2-acetolactate as substrate. In Enterobacteri- 
aceae, up to five different AHAS isozymes have 
been found, each with a different complex regu- 
lation pattern (Barak et al., 1987). To date, only 
one enzyme is known in C. glutamicum (Cordes 
et al., 1992; Eggeling et al., 1987). 

Both LTD and AHAS are feedback inhibited 
by isoleucine, while AHAS is also inhibited by 
leucine and valine (Cordes et al., 1992). Further, 
the expression of AHAS seems to be regulated 
by all three branched-chain amino acids 
(Eggeling et al., 1987). Sahm and co-workers 
studied the production of isoleucine from 
2-oxobutyrate in the late 1980s (Eggeling et al., 
1987; Scheer et al., 1 987), and recently described 
the isolation and molecular characterization of 
the genes involved in branched amino acid syn- 
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thesis in C. glutamicum (Cordes et al., 1992; 
Keilhaueret al., 1993; Mockel et al., 1992). The 
genes encoding LTD, AHAS, and isomerore- 
ductase (IRT) were cloned via complementation 
of E. coli mutants defective in one of the isoleu- 
cine biosynthetic genes (Cordes et al., 1992). 
Deletion analysis of the obtained clones and ex- 
pression in a homologous C. glutamicum back- 
ground showed that the /7vB (AHAS) and ilvC 
(IRT) genes were clustered on a 7-kb chromo- 
somal fragment and probably are transcribed in 
the same direction. Southern blot analysis showed 
that the i7vA gene is localized at least 6 kb away 
from the /VvBC genes, indicating that C. gluta- 
micum does not have a clustered ilv operon as 
has been reported for £. coli. 

Analysis of overexpressed LTD showed that 
the enzyme is composed of four identical sub- 
units and exhibits complex sigmoidal kinetics 
(Mockel et al., 1992). The negative allosteric 
effector isoleucine increased both the coopera- 
tivity. as expressed in terms of the Hill coeffi- 
cient (n), from 2,4 to 3.7 and the substrate 
concentration (S 05 ), where half-maximal rates 
are observed, from 21 to 78 mA/ threonine. The 
positive effector valine promoted the opposite 
effects: reduction of S 05 to 12 mA/ threonine 
and no cooperativity (n = 1 ). Sequence analysis 
of the (7vA gene revealed an open reading frame 
encoding a polypeptide of 436 amino acids with 
an M r of 46,599 Da, which is in agreement with 
the molecular mass determined for the subunit 
of the native enzyme. Comparison of the de- 
duced (7vA amino acid sequence with sequences 
of other threonine deaminases showed that the 
C. glutamicum protein has a considerably re- 
duced carboxy terminal part. On the basis of 
this size reduction and structural analogies with 
other threonine-con verting enzymes, the authors 
postulated an important role for the carboxy 
terminal part in the allosteric control of threo- 
nine deaminase (Moeckel et al., 1992). DNA 
sequence analysis and marker exchange mu- 
tagenesis of the /7vBC operon revealed that the 
AHAS protein is encoded by two genes, /7vB 
and i/vN, as has been observed in other bacteria 
(Keilhaueret al., 1993). Northern RNA blot analy- 
sis of the t/vBNC cluster identified three different 
transcripts of 3.9, 2.3, and 1 . 1 kb, corresponding to 
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//vBNC, z/vNC, and ilvC, respectively. The ilvC 
transcript enccxling isomeroreductase was by far 
the most abundant. Deletion analysis of the 
operon revealed the presence of three active 
promoters, which, however, were not further 
characterized. Addition of 2-oxobutyrate to the 
medium resulted in increased amounts of 
//VBNC and z/vNC transcripts, which correlated 
well with the observed increase in AHAS activ- 
ity. The presence of an z/vNC transcript has not 
been observed in other bacteria. It was there- 
fore concluded that separate synthesis of the 
small ilvN subunit of AHAS might be neces- 
sary for maximal AHAS activity and could 
contribute to the allosteric regulation of the 
enzyme, as has been postulated for E. coli 
(Weinstock et al M 1992). 

B. Transport of Isoleucine 

1. Excretion 

In contrast to the situation described for 
glutamate excretion in C. glutamicum, no spe- 
cial culture conditions or difficult treatments 
are necessary to trigger the efflux of isoleucine. 

C. glutamicum accumulates and secretes sub- 
stantial amounts of isoleucine when the precur- 
sor 2-oxobutyrate is added to the medium 
(Ebbinghausen et al., 1989; Eggeling et al., 
1987). As described above, 2-oxobutyrate in- 
duces AHAS, thereby promoting its own con- 
version, and 2-oxobutyrate bypasses the regu- 
lation of LTD. Similar to glutamate transport, 
three mechanisms have been proposed for 
isoleu-cine secretion into the medium: passive 
diffusion, inversion of the uptake system, and 
active secretion by a specific carrier. 

Ebbinghausen et al. (1989) found that 
C. glutamicum cells were able to export isoleu- 
cine against an existing gradient, which makes 
it unlikely that the internal/external isoleucine 
ratio can be the driving force for isoleucine 
efflux. Further, the transport of isoleucine was 
not coupled with the movement of K+ ions, 
showing that the cells are not "leaky". The 
observed difference between isoleucine excre- 
tion and uptake rates ruled out inversion of the 
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uptake system. No amino acids other than iso- 
leucine were excreted, indicating that a specific 
carrier is involved in the transport of isoleu- 
cine. Ebbinghausen et al. ( 1 989) further showed 
that the membrane potential was the driving 
force for isoleucine excretion, making the spe- 
cific carrier model even more likely. 

2. Uptake 

Like other Gram-positive bacteria, C. gluta- 
micum possesses a single transport system for 
the uptake of all three branched-chain amino 
acids (Ebbinghausen et al., 1989b). As C. gluta- 
micum is fully capable of synthesizing enough 
branched-chain amino acids for growth, the role 
of this low-activity transport system in C. gluta- 
micum is not clearly understood. However, 
Ingraham et al. (1983) have suggested that such 
a system could function as a cyclic retention 
mechanism to scavenge amino acids leaking out 
of the cell. To prove the presence of separate 
uptake and excretion system and to study the 
regulation of the uptake system, Boles et al. 
(1993) constructed defined ilvA mutants via 
marker exchange mutagenesis. These mutants 
are dependent on the addition of isoleucine or 
2-oxobutyrate to the medium. It was proposed 
that these defined mutants could be used to screen 
for mutants defective in isoleucine uptake by 
screening for growth in the presence of 
2-oxobutyrate and nongrowth when supplied with 
limiting isoleucine. However, screening of more 
than 20,000 "colonies did not result in mutants 
with the expected phenotype. The failure of the 
selection procedure was discovered when the 
isoleucine uptake rate of /7vA mutants was deter- 
mined in the presence or absence of 2-oxobutyrate 
and branched-chain amino acids. It was found 
that increasing the internal isoleucine concentra- 
tions greatly enhanced the activity of the uptake 
system. The uptake system became fully induced 
when the internal isoleucine concentration rose 
above 4 to 8 mM. Addition of chloramphenicol 
to short-term fermentations with itvA mutants 
prevented induction of the uptake system, indi- 
cating that the observed induction is due to 
de novo synthesis of the carrier protein, and not 
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caused by regulation of the carrier activity. The 
availability of such an induction system will 
greatly contribute to the design of isolation pro- 
cedures for the specific isoleucine carrier. Boles 
et al. (1993) explained the observed regulation 
with the cycle retention model of Igraham et al. 
(1983), which postulates that amino acid uptake 
systems with low uptake rates function as an effi- 
cient way to scavenge amino acids that are lost 
by nonspecific efflux from the cytosol. 

VIII. PHENYLALANINE, TYROSINE, 
AND TRYPTOPHAN 

A. Molecular Organization of 
C. glutamicum Genes Involved in 
Aromatic Amino Acid Biosynthesis 

The biosynthetic pathway of the aromatic 
amino acids phenylalanine, tyrosine, and tryp- 
tophan in Corynebacterium and other organisms 
involves the participation of at least 17 different 
enzymes and shows complex regulatory mecha- 
nisms (Figure 3). The properties of the enzymes 
and some of the regulatory aspects have been 
reviewed elsewhere (Shiio, 1986). Here, we dis- 
cuss the isolation and characterization of 
C. glutamicum genes involved in the biosynthesis 
of aromatic amino acids and their application in 
the construction of production strains via recom- 
binant DNA technology. Biosynthesis of the aro- 
matic amino acids in C. glutamicum proceeds via 
a common pathway to chorismate, from which 
the pathways to phenylalanine, tyrosine, and tryp- 
tophan diverge. The control of the metabolite flow 
through this pathway primarily occurs via end- 
product inhibition of four enzymatic steps: 
3-deoxy-D-arabino-heptulsonate-7-phosphate 
(DAHP) synthase (DS), anthranilate synthase 
(AS), chorismate mutase (CM), and prephenate 
dehydratase (PD) (Figure 3). The first enzyme of 
the common pathway, DAHP synthase, catalyzes 
the condensation of PEP and erythrose-4-phos- 
phate to DAHP. This enzyme is synergistically 
inhibited by a combination of phenylalanine and 
tyrosine. Its expression is somewhat repressed by 
tyrosine (Shiio, 1986). DAHP synthase forms a 
polypeptide complex with chorismate mutase, 
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which catalyzes the conversion of chorismate to 
prephenate. Therefore, it is not surprising that the 
chorismate mutase- DAHP synthase complex is 
also inhibited by a combination of tyrosine and 
phenylalanine. PD and AS, which initiate the 
pathways to phenylalanine and tryptophan, are 
subject to end-product inhibition by phenylala- 
nine and tryptophan, respectively. 

The genes encoding both wild-type feedback- 
sensitive and deregulated Afunctional DAHP 
synthase-chorismate mutase have been cloned by 
complementation of auxotrophic mutants (Follettie 
and Sinskey, 1986; Ikeda and Katsumata, 1992). 
Recently, the wild-type am gene encoding DHAP 
synthase was sequenced and showed strong ho- 
mology with the aro genes of E. coli and Salmo- 
nella typhimurium (Chen et al., 1993). Compari- 
son of the wild-type sequence with the sequence 
of genes encoding deregulated DAHP synthase 
might reveal more details about binding sites or 
the inhibition mechanism. In addition to the de- 
regulated DHAP synthase-chorismate mutase 
genes, Ikeda and Katsumata (1992) also isolated 
the gene encoding deregulated prephenate 
dehydratase. The three deregulated genes were 
assembled on one multicopy plasmid pKFI and 
used to transform strain KYI 0865. Strain 
KYI 0865 normally produces 7 g/1 of tryptophan, 
but when transformed with plasmid pKFI, the 
metabolite flow was successfully redirected to- 
ward the synthesis of phenylalanine, resulting in 
the production of 10 g/1 of phenylalanine. The 
effect was even more pronounced in fed batch 
fermentations, where the parent strain produced 
18 g/1 of tryptophan and the transformant secreted 
28 g/1 of phenylalanine. This example shows that 
metabolite flow in an organism can efficiently be 
redirected by amplification of deregulated branch- 
point enzymes. 

The tryptophan-specific pathway starts with 
conversion of chorismate to anthranilate by the 
trpE gene product anthranilate synthase. The en- 
zyme is subject to feedback inhibition by tryp- 
tophan. Both wild-type and deregulated tryptophan 
genes form C. glutamicum and C. lactofermentum 
have been isolated via complementation of E. coli 
tryptophan auxotrophs (Del Real et al., 1985; 
Heery and Dunican, 1993; Follettie and Sinskey, 
1986; Matsui eta!., 1986, 1987). Similar to E. coli, 
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FIGURE 3. Biosynthetic pathway for aromatic amino acids. DHAP, 3-deoxy- 
D-arabino-heptulose-7-phosphate. 



the trp genes of C. glutamicum are clustered in 
the order trpE, frpD, trpQ, trpB, and trp A (Heery 
and Dunican, 1993). The trpE gene is preceded 
by the trp promoter and attenuator (Matsui et al., 
1987; Sano and Matsui, 1987). The trp attenuator 
consists of a 17-codon open reading frame con- 
taining three tandem trp codons and is located 
immediately downstream of the trp promoter. 
Sequence analysis of the C. lactofermentum trpE 
gene revealed that the feedback resistance of the 
enzyme is caused by a missense mutation in codon 
38 of the trpE gene (Matsui et al., 1987). How- 
ever, Heery and Dunican (1993) did not observe 
any difference in the deduced amino acid se- 
quence of the anthranilate synthase gene of a 



feedback-resistant C. glutamicum strain. Instead, 
they found that this C. glutamicum strain con- 
tained a nonsense mutation in the trp leader se- 
quence that could be responsible for a constitu- 
tive an ti termination response. 

In order to improve tryptophan production 
by strain KY 10894, Katsumata and Ikeda ( 1 993) 
applied metabolic pathway engineering to this 
strain. Successive introduction of the genes en- 
coding deregulated DAHP synthase, deregulated 
anthranilate synthase, and phosphoribosyl-trans- 
ferase on a multicopy plasmid ultimately re- 
sulted in a 54% increase in tryptophan produc- 
tion, to 43 g/l. In contrast to tryptophan-producing 
strains of Bacillus and E. coli, this C. glutamicum 
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did not show reduced tryptophan secretion rates 
during the stationary phase. The authors ascribed 
the absence of this decline to the abundant pres- 
ence of stable DAHP synthase. Because high- 
level tryptophan production by this strain was 
still dependent on the presence of selection pres- 
sure by spectinomycin, the authors will have to 
ensure plasmid stability in the absence of the 
antibiotic or use conjugal transfer to stably inte- 
grate multiple copies of the deregulated genes 
into the chromosome. 



IX. CONCLUSIONS AND PROSPECTS 

At the moment, nearly all modern molecular 
biology techniques are available for C. glutamicum. 
The use of these techniques in the manipulation of 
Corynebacterium strains will lead to the construc- 
tion of new and hopefully better production strains 
for primary metabolites. 

Recent studies with Ppc- and Py/c-defective 
mutants of C. glutamicum have put the regulation 
of carbon flow at the PEP branch point in a new 
perspective. Extensive metabolic flux analysis stud- 
ies with double or even triple mutants will be 
necessary to evaluate the contribution of the indi- 
vidual enzyme systems to the synthesis of OAA 
and to obtain a more complete understanding of the 
regulation of central metabolism in C glutamicum. 
This knowledge will be helpful in diverting precur- 
sor supply toward the biosynthesis of aromatic and 
aspartate-derived amino acids. 

The studies with gdh' mutants of C gluta- 
micum showed that the GS/GOGAT system has 
enough compensatory capacity to synthesize 
glutamate, and that the excretion of glutamate 
might be limited by the secretion process. Kramer 
and co-workers (1993) have ended the contro- 
versy of the glutamate secretion in favor of ac- 
tive transport. Availability of the genes encod- 
ing glutamate and lysine transport will greatly 
contribute to our understanding and ability to 
manipulate amino acid secretion processes in 
C. glutamicum. 

Studies of the ask-asd operon have indicated 
that regulation of aspartokinase may be mediated 
by the P-subunit of ASK. Genetic engineering of 
the coding region of this subunit can provide 
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completely deregulated aspartokinases, which are 
an absolute prerequisite for high lysine and threo- 
nine production strains and for a better under- 
standing of lysine- and threonine-binding sites. 

The integration and amplification of a deregu- 
lated homoserine dehydrogenase gene showed that 
lysine-secreting strains can be manipulated to pro- 
duce threonine. Before threonine hyperproducing 
C. glutamicum strains can be obtained, the threo- 
nine degradation mechanisms will need to be un- 
raveled and inactivated. 

The genes for the biosynthesis of isoleucine 
and the aromatic amino acids have been cloned 
and their expression studied. The application of 
genes encoding feedback-resistant enzymes of the 
aromatic pathway has resulted in improved pro- 
duction strains for phenylalanine and tryptophan. 
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The two main contributions of this article are the solidi- 
fication of Corynebacterium glutamicum biochemistry 
guided by bioreaction network analysis, and the deter- 
mination of basal metabolic flux distributions during 
growth and lysine synthesis. Employed methodology 
makes use of stoichiometrically based mass balances 
to determine flux distributions in the C. glutamicum 
metabolic network. Presented are a brief description 
of the methodology, a thorough literature review of 
glutamic acid bacteria biochemistry, and specific results 
obtained through a combination of fermentation stud- 
ies and analysts-directed intracellular assays. The latter 
include the findings of the lack of activity of glyoxylate 
shunt, and that phosphoenolpyruvate carboxylase {PPC) 
is the only anaplerotic reaction expressed in C. glu- 
tamicum cultivated on glucose minimal media. Network 
simplifications afforded by the above findings facilitated 
the determination of metabolic flux distributions under 
a variety of culture conditions and led to the following 
conclusions. Both the pentose phosphate pathway and 
PPC support significant fluxes during growth and lysine 
overproduction, and that flux partitioning at the glucose- 
6- phosphate branch point does not appear to limit lysine 
synthesis. © 1993 John Wiley & Sons, Inc. 
Keywords: flux analysis * metabolic engineering * Cory- 
nebacterium glutamicum • lysine 

INTRODUCTION 

Standard practice for strain improvement of industrial mi- 
croorganisms involves exposing a product-producing strain 
to a mutagen, then selecting for those mutants which exhibit 
a desired phenotype and testing the derived strains for 
improved productivity or yield. 13 In this type of procedure, 
the alterations of the parent strain's genome cannot be easily 
resolved due to the nonspecific effects of the mutagens 
employed, so that the genetic modifications that lead to 
improved strains may go undetermined or be incorrectly 
attributed to the phenotype selected for. Although genetic 
engineering techniques 3 can be used to remove the ambigui- 
ties associated with nonspecific mutagens, the enzymes that 
should be modified to increase product yield or productivity 
are often unknown, so that the extra effort required by 
genetic engineering techniques may not be rewarded with 
an equivalent improvement in strain enhancement. This is 
especially true when strain improvement requires modifi- 
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cations of the primary metabolic enzymes which are not 
obviously associated with the product of interest. Conse- 
quently, mutation-selection is still often a faster approach 
for strain improvement. As case in point involves lysine 
production by glutamic acid bacteria. 25 

It is well established that strains of glutamic acid 
bacteria (typically Corynebacterium glutamicum, Brevibac- 
terium flavum, and Brevibacterium lactofermentum) can 
be obtained that produce significant amounts of lysine 
(ca. 30% molar yield) by selecting for mutants that 
exhibit S-(2-aminoethyl) L-cysteine resistance (AEC, a 
lysine antimetabolite), which dramatically enhances se- 
lecting strains with feedback-resistant aspartate kinase. 11 * 74 
Deregulation of lysine synthesis, however, does not result 
in the maximum yield (75% molar), 68 so that mutation- 
selection techniques have focused on modifying enzymes 
of the primary metabolism, such as citrate synthase, 59 
the pyruvate dehydrogenase complex 44,75 (PDC), pyruvate 
kinase 61,65 (PK), and phosphoenolpyruvate carboxylase 88 
(PPC). Although these modifications have resulted in a 
strain that produces lysine at a 55% molar yield, 62 it is not 
known in what respect(s) this strain differs from the original 
parent strain. Indeed, genetic engineering approaches have 
yet to make significant improvements in lysine yield 1 1,32,50 
due to the inability of deciphering the results obtained from 
mutation-selection. 

Although techniques such as metabolic control theory, 15 
biochemical systems theory, 67 whole-cell kinetic models, 23 
stable 21 and radio 6 isotope tracers, and linear analysis 45 
have been applied to metabolic networks, significant in- 
formation regarding potential metabolic bottlenecks can be 
obtained from mass balance techniques. 36 * 38,46 - 76,80,81 - 83 The 
latter are readily applied and do not require information 
regarding enzyme kinetics. In this study we utilized these 
mass balance techniques in conjunction with intracellular 
assays to solidify the biochemistry of C. glutamicum ATCC 
21253, and to examine metabolic flux distributions during 
growth and lysine overproduction to establish the basal case 
for subsequent perturbation studies. 

MATERIALS AND METHODS 
Microorganism and Medium 

The organism C. glutamicum ATCC 21253, obtained from 
the American Type Culture Collection (ATCC, Rockville, 
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The extensive use of l3 C enrichments in precursor metabolites for 
flux quantification does not rely on NADPH stoichiometrics and 
can therefore be used to quantify reducing power fluxes. As an 
application of this concept, the NADPH fluxes were quantified 
in an L-lysine producer of Corynebacterium glutamicum grown 
into metabolic and isotopic steady state with [l- l3 C]glucose. In 
this case, where the organism's NADPH-dependent glutamate 
dehydrogenase consumes reducing power, the NADPH flux 
generated is 210% (molar flux relative to glucose uptake rate) with 
its major part (72% of the total) generated via the pentose 
phosphate pathway activity. An isogenic strain in which the 
glutamate dehydrogenase of C. glutamicum was replaced by the 
NADH-dependent glutamate dehydrogenase of Peptostreptococcus 
asaccharolyticus was made and the metabolite fluxes were again 
estimated. The major response to this local perturbation is a drasti- 
cally reduced NADPH generation of only 139%. Most of the 
NADPH (62% of the total) is now generated via the tricarboxylic 
acid cycle activity. This shows the extraordinary flexibility of the 
central metabolism and provides a picture of the global regulatory 
properties of the central metabolism. Furthermore, a detailed 
analysis of the fluxes and exchange fluxes within the anaplerotic 
reactions is given. It is hypothesized that these reactions might also 
serve to balance the total reducing power budget as well as the 
energy budget within the cell. C i w Academic Pr*« 

Key Words : reducing power supply; lysine production; anaplerotic 
reactions; pentose phosphate pathway; flux quantifications; isotope 
studies. 



1. INTRODUCTION 

The fueling reactions within the central metabolism serve 
to provide energy as well as precursor metabolites for the 
cell. This dual function couid be of particular importance 
when metabolite overproduction is attempted, since removal 
of precursors and an altered energy budget require a read- 
justment of the fluxes in this part of the metabolism. 
Metabolic engineering tries to make specific predictions of 
where increased activities could lead lo a product increase. 
Such predictions are corroborated for reactions within 



assembling pathways, where allosterically controlled 
enzymes and branching points are major targets for obtain- 
ing flux increase (Eggeling et aL 1996). Examples are the 
successful construction of amino-acid-producing-bacteria 
such as Corynebacterium glutamicum (Eggeling el £//., 
1997a) or Escherichia co!i (Eggeling and Sahm. 1998). 
Moreover, as a further target for improving metabolite 
production, a specific product exporter has recently been 
identified (Vrljic el a!., 1996). This completes the identifica- 
tion of the cellular control steps for converting precursors of 
the central metabolism into the desired extracellular 
product. However, in contrast to the assembling pathways 
and the transport process, the targets in the central 
metabolism to be engineered are less well defined. This is 
due to ( i ) the many branching points and cyclic and parallel 
reactions, (ii) the not easily accessible tools for quantifying 
the in vivo fluxes, (iii) the fact that enzyme activities in the 
network may be overlooked, and last but not least (iv) the 
very flexible responses of the living cell. Therefore, it is still 
a very ambitious goal to predict targets for fiux increase in 
this part of the metabolism. A recent success in engineering 
the fiux of glucose-derived carbons toward aromatic amino 
acid formation was based on more qualitative physiological 
knowledge {Flores et a/., 1996) and the successful applica- 
tion in producer strains has not yet been documented 
(Barry, 1996). 

Of course, one aspect of metabolic engineering is a quan- 
tification of all fluxes. We have developed a comprehensive 
methodology distinguished by (i) the analysis of fluxes in 
the steady state, (ii) the intensive use of carbon atom quan- 
tification, made possible through [ l- 13 C]glucose and NMR 
spectroscopy, and (iii) metabolite balances (Marx et aU 
1996). This approach is completed by (iv) sophisticated 
mathematical tools also allowing the estimation of reverse 
fluxes and statistical parameters (Wiechert et «/., 1997). In 
particular, the high resolution obtained by such analyses in 
formulating new questions, whose answers could then in the 
future be utilized for metabolic engineering. This is the spe- 
cial benefit of intracellular flux analysis. In our previous 
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investigations, we simply quantified the fluxes in isogenic 
C. glutamicum strains in a cause-effect study, where the cen- 
tral metabolism itself was not engineered, but only the final 
biosynthesis sequence or the export (Marx et aL 1996. 
1997). This served to probe the flexibility of the fueling 
reactions in response to different carbon burdens. 

However, in most flux analyses no access to the reducing 
power budget has been possible to date. In contrast. 
NADP(H) stoichiometrics are generally used as additional 
constraints to make solutions for carbon fluxes possible at 
all. In contrast, our approach, depending heavily on infor- 
mation on the fate of carbon atoms, does not require these 
constraints. Access to the reducing power budget is thus 
possible, since some carbon fluxes, e.g.. through the glucose- 
6-phosphate dehydrogenase or the isocitrate dehydrogenase 
reaction, are inevitably linked to the NADPH flux. There- 
fore, carbon fluxes can be used to derive NADPH stoichio- 
metrics. We have already performed such calculations in 
previous work to a certain extent (Sonntag et 1995). In 
the present paper, we describe our direct intervention in the 
NADPH budget to once again quantify the flexibility within 
the fueling reactions in a cause-effect study. Whereas we 
previously quantified the response of the central metabolism 
due to different carbon flux burdens ( Marx et aL, 1997), this 
time we focus on the response of the central metabolism to 
different reducing power demands. 

MATERIALS AND METHODS 

Genetic Engineering 

The Pepto.sirepJococcus asacclutrolyticus glutamate dehy- 
drogenase gene was cloned according to the published 
sequence (Snedecor et aL 1991) by use of the polymerase 
chain reaction and the primers 5'-ATGGATCCCCGAGT 
GAGAAATCACGGTG-3' and 5'-CTGAGCTCAGCTGA 
TTATATGAGTTT-3'. The 1.4-kb product obtained was 
digested with BamH\ and Sad and ligated with the shuttle 
vector pEKEx2 treated identically (Eikmanns et aL 1991 ). 
Transformation of E. colt DH5a yielded kanamycin-resis- 
tant clones which were confirmed by restriction analyses to 
contain pEKExpgdh of 9.7 kb (Fig. 1 ). With this vector the 
specific glutamate dehydrogenase activity (NADH as sub- 
strate, E. coli as host) was increased from 0.05 to 1.1, and, 
after isopropyl-/y-i>thiogalactoside addition, to 5.4 //mol/ 
min and mg protein. 

For construction of the gene deletion vector pK\9moh- 
sacBJgdh. part of the gdh gene was isolated as a 1.25-kb 
Sphl fragment (Bormann et al.. 1992) and ligated with 
pUC18. A BglU and EcoRV site within gdh was used to 
remove a 61 1 -bp internal fragment to result in pUCl&Jgdfi. 
The Sph \ fragment, now only 0.6 kb in size, was reisolated 



Stars el ill. 



Sphl 




FIG. 1. The expression vector pEKEXpgdh to express the N A In- 
dependent PvpiosfnptoiotTU* asuuhurolyiivus glummaie dehydrogenase 
gene jfrf/i in C. gltnuniiaim. Abrevialions; Kuth kanamycine resistance: 
tacf. the constitutive!}- expressed repressor of the hit operon; W(7\ Ttic 
promotor. 

and ligated with the vector pK\9mohsacB, which had been 
previously cleaved with Sphl and treated with alkaline 
phosphatase. The resulting vector pKl9mohsacBAgdh was 
introduced into the mobilizing donor strain E. coli SI 7-1 
(Simon et al.. 1983), and matings were performed with C. 
glutamicum MH20-22B ( Schrumpftv 1992) according to 
Schafer et al (1990). The single transconjuganl obtained 
was cultivated overnight in liquid LBG for proliferation, 
and then dilutions were spread on LBG containing 10% 
sucrose (Schafer et al.. 1994) to select for clones where the 
second recombination had taken place. One hundred Suc r 
clones were assayed for their kanamycin phenotype. and 
eventually 15 Suc r Kan'' clones were subjected to glutamate 
dehydrogenase activity determinations. 

The enzymes necessary for in vitro recombinations were 
obtained from Boehringer ( Mannheim, Germany ) and used 
as recommended. Plasmids were isolated as described 
( Schwarzer and Puhler. 1990). E. coli was transformed after 
pretreatment with CaCU/RbCU and C. glutamicum by 
electroporation (Liebl et al., 1989). 

Continuous Culture 

The strains (homologous mutant strain C. glutamicum 
MH 20-22 hjgdh pBKl.9gdh-l and heterologous mutant 
strain C. glutamicum MH2Q-22Bdgdh pEKExpgdh) were 
grown in continuous culture as described elsewhere ( Marx 
et aL 1996). The medium was identical to that previously 
used with 43 mM glucose and 1.4 mM L-leucine. except that 
the ammonium sulfate concentration was 600 mM. and 
isopropyl-//-i>-thiogalacioside ( 1PTG ) (0.5 mM ) and methi- 
onine S-sulfoximine ( 14.8 mM ) were present. After inocula- 
tion itnd consumption of glucose, the medium together with 
[ l-°C]glucose was fed continuously at a rale of 0.082 h' 1 . 
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table: i 

Extracellular Fluxes in Two L-Lysine- Producing Strains of 
C.gluiamicum Differing Only in the Possession of an NAD or 
NADPH-Dependent Glutamate Dehydrogenase 

Gluiumaie dehydrogenase 



Ca rbon" X<lh\ N AD I .(,- f //;| N A D P I 



Glucose 100 100 

Biomass 33 28 

Product IK 32 

CO ; 45 .W 

Sum % W 

) p S (mo!p nn>ls 'l O.IK 0.32 

i/5 1 nmol s • min " 1 ■ mg^i, ) 23.0 26.9 



Wi. Abbreviations: i>dhi NAD) are the data for the strain with 
the NAD-dependent glutamate dehydrogenase of Peptosfrepioanais 
asaivharolytktis. yrf/if NADP) those for the strain with the NAD-depend- 
em enzyme. )' xs ( )' PS i is the yield of conversion of the substrate glucose 
• Si to biomass (XI. and product (P). respectively, f/ s is the specific sub- 
strate uptake rate. 

J Expressed as percentage of consumed glucose carbon. 

C0 2 evolution as well as the dissolved oxygen concentra- 
tion of the culture, was recorded on-line. In the steady state, 
glucose was entirely consumed. Samples of the cultures were 
collected to estimate lysine, growth. NADH- and NADPH- 
dependent glutamate dehydrogenase activities, and gluta- 
mine-oxoglutarate aminotransferase activity. In addition, 
in the steady state the 13 C0 2 content of the effluent gas was 
quantitated by nondispersive infrared spectroscopy. From 
the analyses, the carbon balances of the cultures were 
derived which are given in Table 1 as percentages of flux 
relative to the glucose input rate. 

Enzy me Da en i j ina turns 

Cells were harvested from the chemostat or exponentially 
growing cultures. They were washed with 0.9% NaCl. 
resuspended in the respective assay buffer, and disrupted by 
bonification. The cell debris was removed by centrifuga- 
tion. and the resulting crude extract was used directly 
for glutamate dehydrogenase activity determinations. For 
determination of the glutamine-oxoglutarate aminotrans- 
ferase activity (GOG AT), crude extract was gel-filtrated 
iPD-10 columns. Pharmacia! to remove ammonium salts, 
using as eluent buffer 0.1 M potassium phosphate (pH 7.5), 
1.44 mM cysteine. 0.5 mM EDTA. 

The glutamate dehydrogenase activity was assayed by 
following the oxoglutarate-dependem decrease in the extinc- 
tion at 340 n m. The assay contained, in I ml. 100/imol 
Tris HCI buffer (pH 8.0). 0.25 //mol NADPH <or NADH. 
respectively |. 20 //mo) NH 4 CI. and 10// mol 2-oxoglutarate. 
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The GOGAT activity was assayed in a system containing, in 
I ml. 100 //mol Tris/HCI buffer (pH 7.6). 0.1 //mol dithio- 
threitol. 0.25 //mol NADPH, 10//mol 2-oxoglutarate, and 
10 //mol glutamine. The protein concentration in the 
extracts was determined by the biuret method using bovine 
serum albumin as standard. 

Determination of Fractional Enrichments in Carbon Atoms 

High-resolution 'H NMR spectra were obtained on an 
AMX-400 WB spectrometer { Bruker. Karlsruhe. Germany) 
operating at 400. 1 3 M Hz and equipped with a multichannel 
interface and a 5-mm inverse GRASP probe head. For the 
determination of 13 C enrichments in the protonated carbon 
atoms, as well as in the nonprotonated carbon atoms. 
l H NMR difference spectroscopy with and without GARP- 
1 composite pulse ,3 C decoupling was used. The 13 C enrich- 
ments of protonated carbons were calculated as the ratio of 
the n C satellite signal area in the difference spectrum to the 
total integral of each proton in the GARP-1 decoupled 
spectrum. For the determination of ,3 C enrichments in the 
nonprotonated carbons, spin-echo NMR spectroscopy was 
applied (Wendisch et «/., 1997). Typical settings of the 
inversion pulse sequence were 15-s relaxation delay. 2-s 
H 2 0 presaturation pulse; 9//s 90° transmitter pulse; and 
256-//S sine-shaped spoiler delay, with a selective 18-//S 1 80° 
transmitter pulse, and a 180° decoupler inversion pulse 
adjusted sample-specific. GARP-1 decoupling was optionally 
used. 

The 12 C0 2 concentration in the exhaust gas was measured 
on-line by non-dispersive infrared spectroscopy lURAS 
10E, Hartmann & Braun, Frankfurt, Germany). The 13 C0 2 
concentration was quantified ofT-line by infrared spectros- 
copy with appropriate standards and settings {COI, Fischer 
Analysen GmbH. Leipzig. Germany). 

Flux Determinations 

The precursor demand for cell mass formation is as 
described (Marx et aL 1996), as well as the biochemical 
model used ( Marx et aL 1997). The mathematical treatment 
has been described in detail by Wiechert et ai ( 1997) using 
matrix notations of the carbon-labeling balance equations. 
The program NMRFlux was used (Wiechert and de Graaf. 
1997) to convert the biochemical equations and the equa- 
tions of carbon atom transitions into the corresponding flux 
equations. Similarly, the experimentally determined fluxes 
and enrichments were incorporated. Using a sophisticated 
parameter-fitting algorithm ( Wiechert et aL 1997). which is 
based on the Levenbcrg-Marquard algorithm, a nonlinear 
least-squares analysis with mutiplc initializations was 
performed to determine the global minimum. 
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RESULTS 

Et Jgineering of Stra ins 

The rationale of the NADPH flexibility study in 
C glutumkum was to use an NADH-consuming glutamate 
dehydrogenase instead of the original NADPH-dependent 
enzyme of the organism (Oshima ei a!., 1964; Bormann et 
aL 1992). Such an NADH-consuming enzyme is present in 
P. asaccharolyticus, where it serves catabolic purposes 
(Snedecor et a/., 1991). Of course, when dictated by the 
equilibrium constant, the enzyme also catalyzes the forma- 
tion of glutamate. Based on the published sequence and 
using the polymerase chain reaction, we amplified the 
P. asaccharolyticus glutamate dehydrogenase gene from 
genomic DNA. A Bam HI site was introduced 39 bp in front 
of the gene and a Sad site 163 bp downstream of it. This 
served to clone the amplified fragment into the correspond- 
ing sites of the expression vector pEKEx2 ( Eikmanns et aL 
1991} to yield pEKExpgdh. With this vector, NAD- 
dependent and IPTG-inducible glutamate dehydrogenase 
activity was obtained in C. glutamicum (see below). 

To delete the original NADPH-dependent enzyme from 
the C glutamicum chromosome, the mobilizable but non- 
replicative vector pK\9mobsacBdgdh was constructed (see 
Materials and Methods). It carries the homologous 
C glutamicum gene, from which an internal fragment had 
been deleted. It was introduced into the lysine producer 
C. glutamicum MH20-22B by intergeneric mating. This 
lysine producer was chosen because (i) the NADPH 
demand for lysine production is high, (ii) the strain is very 
well characterized (Schrumpf et aL 1992 ), and (iii ) we have 
already established several detailed flux quantifications for 
this high-level producer (Marx et aL 1996, 1997). After 
mating, one transconjugant (resistant for kanamycin) was 
obtained with the vector integrated in the chromosomal gclh 
due to homologous recombination. This clone was sub- 
jected to a second round of positive selection ( resistance for 
sucrose) to result in clones which had lost the vector. Fifteen 
of the numerous clones obtained were assayed for glutamate 
dehydrogenase activity, and three of them were identified as 
having no enzyme activity due to the appropriate location 
of the second recombination. The loss of the chromosomal 
fragment in the constructed strain C glutamicum MH20- 
22BJgdh was confirmed in a Southern blot (not shown). 
This strain was transformed with plasmid pEKExpgdh. 
encoding the heterologous P. asaccharolyticus enzyme. As 
the control, the deletion strain was transformed with 
pEK1.9s<///-l. encoding the homologous C. glutamicum 
enzyme ( Bormann et aL 1992). In the following the strains 
will be referred to as the heterologous mutant and the 
homologous mutant, respectively. 
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Gene Expressions 

To verify expression of the glutamate dehydrogenase 
genes, enzyme activity determinations were performed. In 
the strain MH20-22BJ^//; with deleted dehydrogenase, no 
enzyme activity was detected (not shown). In the heterolo- 
gous mutant with the P. asaccharolyticus gene, the NADH- 
dependent enzyme was synthesized (without IPTG addi- 
tion) to yield a high specific dehydrogenase activity of 
1.80//mol/min/mg protein. In addition to the dehydro- 
genase activity residual glutamine-oxoglutarate amino- 
transferase activity was detectable in this heterologous 
mutant (activity of 0.033 //mol/min/mg protein); it was 
almost zero in the homologous mutant. Since this activity, 
together with the glutamine synthetase, could contribute 
to NADPH-dependent glutamate synthesis, the inhibitor 
methionine S-sulfoximine was used ( Weisbrod and Meister, 
1973), In vivo this inhibitor specifically prevents glutamate 
synthesis via the transferase/glutamine synthetase system, 
since growth of C. glutamicum with inactivated gclh is 
abolished at a concentration of 10 mM, whereas growth of 
the wild type is not (M. Tesch, personal communication). 
Consequently, both the heterologous mutant and the 
homologous mutant were grown in the continuous cultures 
in the presence of 15 mM methionine S-sulfoximine. As 
verified in extracts of cells taken from the cultures during 
establishment of the steady state, the dehydrogenase 
activities in both cultures were already comparable after 3 
and 4 volume changes. In the steady state, a specific activity 
of 5.08^mol/min/mg protein was determined for the NAD- 
dependent enzyme of the heterologous mutant, and a 
specific activity of 6.97//mol/min/mg protein for the 
NADP-dependent enzyme of the homologous mutant was 
determined. 

NMR Quantifications 

Continuous cultures of the mutants were supplied with 
99% enriched [ l- n C] glucose. In the steady state all 
glucose was consumed. Making use of the fact that in the 
metabolic and isotopic steady state the enrichments in the 
precursors within the metabolic network are identical to 
those in the precursors stored as building blocks in the cell 
material, we hydrolyzed the total cell material and isolated 
the building blocks on a preparative scale to allow enrich- 
ment determinations. In total, the fractional enrichments in 
carbon in nine amino acids, in glucosamine, and in the 
sugar moiety of guanosine were quantified by 'H-detected 
NMR spectroscopy, including spin-echo NMR for non- 
protonated carbons of amino acids (Tables 2 and 3). 

As a technique generally applicable to quantifying the 
enrichment in CO.. a method for estimatinc the enrichment 
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TABLE II 

Measured U C Enrichments in Cell Constituents Isolated from the Heterologous Mutant, gdh( NAD) 



Metabolite , l, C enrichment ( % I in carbon 





C-l 


C-2 


C-3 


C-4 


C-5 


C-6 


Gly 


3.2 ±0.3 


2.2 ±0.3 










Alu 


4.3 ± 1.0 


4.5 + 1.0 


37.0 ±0.8 








Gl\ 


18.6 ±0.6 


30.1 ±0.6 


18.3 ±0.2 


37.5 ±0.4 


3.6 ± 1.0 




Arc 


n.d. 


30.7 ± 1.0 


18.4+ 1.0 


35.1 ± 1.0 


2.9+ 1.0 


16.4 ±0.5 


Thr 


12.6 ±0.8 


n.d. 


n.d. 


16.2 ±1.5 






Asx 


n.d. 


18.1 + 1.0 


30.7 ± 1.5 


n.d. 






Phe" 


n.d. 


n.d. 


34.3 ± 2.5 


n.d. 


20.5 ± 2.5 


1.1 ± 1.0 




n.d. 


4.8 ± 1.0 


n.d. 


37.4 ± 1.5 






Lys 


n.d. 


13.2+0.6 


30.2 ± 1.0 


15.5+1.0 


n.d. 


7.6 + 0.4 


Guanosine 


24.5 ±0.8 


\~d. 


n.d. 


n.d. 


25.1 + 1.0 




Glucosamine 


7S.0±4.0 


n.d. 


n.d. 


n.d. 


n.d. 


n.d. 


CO; 


16.5 ±0.5 













Sow. n.d., not determined. 

" The mean value Cor C-5 is derived from C-5 plus C-9. The enrichment in C-7 and C-8 is 1.1 + 1.0%. 
* The mean value for C-4 is derived from C-4 plus C-4\ 



in C-6 of arginine was developed because the C-6 enrichment 
stems exclusively from C0 2 incorporated via carbamoyl- 
phosphate synthase activity. In a series of experiments, the 
spin-echo delay (TE) was varied to quantify the amplitudes 
of the ,3 C satellites in the difference spectra obtained. A 
cosine dependency of the amplitudes of the difference spec- 
tra on TE is observed due to the evolution of the long-range 
J HCNC during the spin-echo delay (Fig. 2). According 
to the difference signal of 30.7% for the heterologous 



mutant, a spin-echo delay of 142 ms and calibration result 
in a fractional enrichment of 16.5+0.5%. This is in 
excellent agreement with the enrichment in C0 2 as deter- 
mined by nondispersive infrared spectroscopy (Table 2). It 
shows that in the culture the intracellular C0 2 is in equili- 
brium with the extracellular C0 2 . Access to the fructose 
6-phosphate enrichments was made possible by the quanti- 
fications in glucosamine. Together with the available 
enrichments in e4p and p5p, this further improves the flux 



TABLE III 

Measured U C Enrichments in Cell Constituents Isolated from the Homologous Mutant, ^rfA(NADP) 



Metabolite 13 C enrichment ( % ) in carbon 





C-l 


C-2 


C-3 


C-4 


C-5 


C-6 


Gl\ 


2.0 ±0.2 


1.7+0.3 










Ser 


n.d. 


n.d. 


23.0+2.0 








Ala 


2.0 ± 0.5 


2.2 ± 1.0 


25.7 ± 0.7 








Glx 


15.7 ±0.6 


21.1 ±0.3 


ll.6±0.3 


25.8 + 0.7 


2.2 ± 1.0 




Arc 


n.d. 


20.5 ±0.5 


11.6± 1.0 


23.9 + 1 .0 


1.9 ± 1.0 


21.3 + 0.4 


Thr 


8.0 ± 0.5 


12.7 ± 1.0 


20.4 ± 1.7 


15.3 ± 1.0 






Phe" 


n.d. 


1.1 ± 1.0 


25.5 i 1.0 


n.d. 


12.8 + 2.0 


1.1 ±1.0 


Yal* 


n.d. 


2.6 ± 1.0 


2.6+ 1.0 


26.0 + 2.5* 






Lys 


n.d. 


8.2 ± 0.4 


22.2 ± 1.(1 


14.8 ± 1.0 


n.d. 


3.7 ± 0.4 


Guanosine 


7.5 ±0.4 


2.9 ± 1.3 


1.1 + 1.0 


1.! ± 1.0 


16.4 ±0.6 




Glucosamine 


62.3 ± 2.5 


n.d. 


n.d. 


n.d. 


n.d. 


n.d. 


CO; 


22.4 + 0.5 













"The mean value for C-5 is derived from C-5 plus C-9. The enrichment in C-7 and C-S is 1.1 + 1.0%. 
h The mean value lor C-4 is derived from C-4 plus C-4'. 
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TABLE IV 

"C Enrichments in Precursor Metabolites of the Heterologous Mutant as Derived from the Experimentally Determined Enrichments 
(Standard Characters) and as Predicted by the Solution of the Mathematical Model for Flux Quantification (in Italics) 



Metabolite '-'C enrichment ( % ) in carbon 





C-l 


C-2 


C-3 


C-4 


C-5 


C-6 


ibp 


78.0 + 4.0 


n.d. 


n.d. 


n.d. 


n.d. 


n.d. 


f6p 


MO.O 


1.4 


4.6 




/..* 


9.6 


p5p 


24.5 ±0.8 


n.d. 


n.d. 


n.d. 


25.1 ±1.0 




p5p 


24.4 


23 


2.4 


1.6 






e4p 


1.1 ±1.0 


1.1 ±1.0 


1.1 ± 1.0 


20.5 ±2.5 






e4p 


2.9 


2.2 


/J 








gap 


3.2 ±0.3 


2.2 ± 0.3 


n.d. 








gap 


3.1 


1.9 


J7.9 








pyr 


4.3+1.0 


4.6 ± 1.0 


37.2 ±1.5 








pyr 


4.0 


3.6 ■ 


37.2 








akg 


18.6 ±0.6 


30.1 ±0.6 


18.3 ±0.2 


37.5 ±0.4 


3.6 ± 1.0 




akg 


163 


30.6 




J7J 


3.7 




oaa 


12.6 ±0.8 


18.1 ± 1.0 


30.7 ± 1.5 


16.2 + 1.5 






IKJU 


11.7 




30.6 


16.3 






lys 


n.d. 


13.2 + 0.6 


30.2 ± 1.0 


15.5 ± 1.0 


n.d. 


7.6 + 0 


Irs 


9.5 


14.0 


32.5 


/<U 


35.3 


7.* 


CO, 


16.4 ±0.5 












CO : 


16.4 













24 




0 40 80 120 160 200 

Spin-Echo relaxation delay TE (ms) 

FIG. 2. Spin-echo difference signal amplitude of '*C-6 satellites ot H-5 
of arginine as a function of the spin-echo relaxation delay, from arginine 
isolated from the homologous mutant ( ■ I and from the heterologous 
mutant i • l after growih on [ l- n C]-glueose. 



determinations within the pentose phosphate pathway 
(PPP). 

The fractional enrichments in the precursor metabolites 
of the heterologous mutant thus derived are given in 
Table 4. For those positions where several data were 
available (for instance, C-l pyr derived from C-l ala and 
C- 1 val ), a weighted mean value was used. Weighting was 
performed according to the measurement error In a similar 
manner, the data for the precursor metabolites of the 
homologous mutant were derived and are given in Table 5. 

Estimated Carbon Fluxes 

The current model used to estimate the carbon fluxes con- 
tains 14 fluxes removing or supplying metabolites from the 
network, with 11 of them being accumulating fluxes for 
biomass synthesis. Another 16 fluxes connect the 14 intra- 
cellular metabolite pools in the central metabolism ( Marx et 
aL 1997). For each experiment a detailed statistical analysis 
was carried out according to the procedures described by 
Wiechert et at. (1997). In particular, confidence regions 
and correlations were computed for the estimated fluxes 
(Table 6). The sums of squared deviations were 32.9 and 
48.7 for the heterologous and the homologous mutant. 



Intratx'lhilar Fluxes muter DUfcrvm NADPH Demands 



MI-TAIKlLIf l-NCilM-J-KISU I. J5-48 ( 

article no. MT9R0106 



TABLE V 



»C Enrichments in Precursor Metabolites of the Homologous Mutant as Derived from the Experimentally Determined Enrichments 
(Standard Characters) and as Predicted by the Solution of the Mathematical Model for Flux Quantification (in Italics) 



Metabolite "C enrichment |%l in carbon 



p5p 

pip 76 



e4p 
e4 P 

gap 
gap 
pyr 



akg I6J 



C-l CO 



1.1 + 1.0 1.1 + 1.0 1.1 ±1.0 12.8 ± 2.0 

22 I* 14 l4J 

2.0 ±0.2 1 .7 ±0.3 23.0 ±2.0 

2.(1 I- 7 263 

2.0 ±0.5 2.4+1.0 25.7 ±0.7 

2.0 /-f 2 ft. 2 



oaa 
ouu 

Ivs 

lys *■< 

C02 21.3 ±0.4 

CO, 232 



C-3 C-4 C-5 C-6 



ftp 62.3 ±2.5 n.d. 

ff>p '-J 

7 5 +04 2.9±1.3 1.1 ±1.0 I.I ±1.0 16.4 ±0.6 

_\/ /.A M 



akg 15.7±0.6 21.1 ±0.3 ll.6±0.3 25.8±0.7 2.2±1.0 



21,2 I tJi 26.2 I* 



n.d. n.d. n.d. n.d. 

2.5 1.4 U 



8.0 ±0.5 12.7+ 1.0 20.4 ±1.7 15.3 ±1.0 

7M H.6 21.2 M2 

nd 8.2 ±0.4 22.2+1.0 14.8 ±1.0 n.d. 3-7 ±0.4 

22.5 H-2 24< * 44 



TABLE VI 

Statistical Analysis of Net or Exchange Fluxes in the Two Isogenic 
C glutamicum Strains. 



Flux rate I •/• ) of lower and upper boundary 



Flux 


Reaction 


Homologous mutant 


Heterologous mutant 


V, 


Glucose-6-P DH 


56.8 


95.9 


21.0 


31.8 




Isodtraie lyase 


0.0 


0.0 


0.0 


5.6 




Diaminopimelatc DH 


I1J 


17.3 


6.2 


9.2 


vr 


Giucose-6-P i some rase 


96.9 


X 


0.0 


X 


V™ 


Glyceraldehyde-3-P DH 


79.6 


X 


36.8 


132.5 


v r 


TK(x5p + e4p « g3p + IISp) 


3.1 


14.9 


3.5 


8.0 


vr 


TA 


0.0 


9.4 


6.4 


20.3 




TK(x5p + ri5p — s7p + g3p) 


54.4 


74.7 


30.0 


46.6 


v,7 


Fumarase 


214 


41.2 


16.2 


50.2 




Anapterotic exchange 


3.8 


17.0 


22.6 


36.5 



Son: The 90% confidence interval for the upper and lower boundary is 
given. The values are expressed as perceniages relative to the glucose 
uptake flux (set to 100%). Abbreviations: DH dehydrogenase; TK. trans- 
kciolase: TA. transaldolase. The fluxes are specified in detail in Marx vt al. 
11996). 



respectively, showing good correspondence between 
measured and computed values. 

The carbon fluxes estimated in the homologous mutant 
are shown in Fig. 3. The flux entering the PPP is as high as 
76%. The excess carbon flux not required for the synthesis 
of polyolphosphates for biosynthetic purposes is over the 
nonoxidative transaldolase and transketolase reactions to 
f6p and gap, which thus rejoin the glycolysis. At the level of 
the gap 162 mol% is present out of the theoretical maxi- 
mum of 200 mol%. This equals 81 % of the maximum. Thus 
at the level of the gap, 19% of the carbon was used for 
biomass formation and carbon dioxide generation. In the 
heterologous mutant (Fig. 4), the flux entering the PPP is 
strongly reduced to 26%. This can be traced back to the 
enrichment in C-3 of pyr (compare Tables 4 and 5), since 
less C-l of the substrate [ l- ,3 C]glucose is removed from the 
network in this mutant. However, at the level of the gap, 
again 176mol% is present. Thus the strong change in the 
oxidative part of the PPP has limited consequences for the 
gap availability, and for the (luxes towards pep and pyr. In 
both mutants the elyoxylic acid cycle carries no flux. 
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Glucose 




FIG. 3. In vivo flux distribution in the central metabolism for C glutamicum MH20-22B using its NADP-dependent gluumate dehydrogenase 
(homologous mutant MH20-22Bdy<#r pEK 1 .9gdh- 1 ). Numbers in oval symbols near thick lines give the estimated net fluxes; those near the thin arrows 
give the measured fluxes required for biomass synthesis. Numbers in hexagons give the estimated additional exchange fluxes in reversible reactions. All 
fluxes are given as molar metabolite flux expressed as a percentage of the glucose uptake rate, which was 26.9//moI/g dry. Abbreviations: g6p, glucose 
6-phosphate; H6bp. fructose 1.6-bisphosphate; ffip. fructose 6-bisphosphate; gap, glyceraldehyde 3-phosphatc; pep, phosphowo/pyruvate; pyr, pyruvate: 
oaa. oxaloacetate; fum, fumarate; icit, isocitrate. akg, ketoglutarate; p5p. pentose 5-phosphate; s7p. sedoheptulose 7-phosphate; e4p. erythrose 
4-phosphate: LL-DAP. LL-diaminopimelatc; DL-DAP. DL-diaminopimelatc. 
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FIG. 4. in vito flux distribution in the central metabolism for C. glutamicum MH20-22B using the NAD-<Iependent glutamatc dehydrogenase or 
P. asuccharofyticus (heterologous mutant C ghuannewn MH20-22B4?<//r pEKEXpgdh ).The glucose flux, set to 100% was 23.0//mol/g dry. Description 
and abbreviations as in Fig. 3. 



whereas the flux entering the tricarboxylic acid cycle (TCC) 
is significantly increased to 86 mol%, compared to the 
homologous mutant (58%). This could be due to an 
increased NADH demand in the heterologous mutant. In 
both mutants a high and comparable part of the lysine is 
synthesized via the dehydrogenase branch of the split 



pathway, which contributes 48 % to total lysine synthesis in 
the homologous mutant and 43% in the heterologous 
mutant. 

The exchange fluxes are shown in hexagons in Figs. 4 
and 5, with their 90% confidence interval given in Table 6. 
For the heterologous mutant we estimated a large backflux 
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Glucose 



TransA 




Lysine 



Glutamate 



TransA 1 

Lysine 



Peptostreptococcus 
TransA GluDH 

y asaccharolyticus 

Glutamate 



FIG. 5. The NADPH-generaiing and -consuming duxes with N AD PH -dependent glutamate dehydrogenase (A) and NADH-dependent glutamate 
dehydrogenase ( B}. The numbers in large ovals are the NADPH (luxes in mol % relative to the glucose uptake flux. 



of 29.2 mol% (confidence interval 22.6-36,5) within the 
anaplerotic reactions. In the homologous mutant this value 
is only 10.3% (interval 3.8-1 7.0). This difference is reflected 
in the higher enrichment in C-2 of pyr in the heterologous 
mutant, which is due to the higher incorporation of oxalo- 
acetate carbons via backflux into pyruvate. On the other 
hand, the enrichment in C-2 of oxaloacetate (oaa) decreases 
with high exchange flux (Fig. 6). In the computer-aided 
analysis all data are considered and contribute to the 
estimate of this exchange flux v~ . Further high-exchange 
fluxes are present within the PPP, with the transketolase 
reaction (p5p + p5p«->s7p + gap> exceeding the net flux 
more than threefold. In the TCC a moderate backflux from 
oxaloacetate to fumarate was estimated. 



77ft' Estimated NADPH Fluxes 

For the NADPH-generating fluxes only three reactions 
are relevant. These are the glucose-6-phosphate dehydro- 
genase. giuconate-6-phosphate dehydrogenase and iso- 
citrate dehydrogenase reactions, where the carbon flux is 
directly equivalent to molar NAPH generation (Fig. 5A). 
These generating fluxes add up in the homologous mutant 
to a molar flux of 210.5%. The malic enzyme was not con- 
sidered because only a maximal value can be derived from 
the current flux analyses (see Discussion). To calculate the 



consuming fluxes we counted the required amounts of 
NADPH from known stoichiometrics of the various biosyn- 
thetic reactions and the cell composition ( Marx et a/., 1996), 
Approximately 1 5,500 // mol NADPH is necessary for the 
synthesis of 1 g dry weight. Of this amount about 50% is 
consumed by the glutamate dehydrogenase reaction, which 
functions mainly to deliver its amino group via transami- 
nation to cell mass synthesis (6350^mol/g dry), and to 
a lesser extent to provide glutamate and glutamine for 
protein synthesis (800/imol/g dry). Based on these stoichio- 
metrics the NADPH-consuming flux for a-ketoglutarate 
reduction is 31 mol%, and that directly required for reduc- 
tion of ceil material via dehydrogenases is 44 mol%. The 
fluxes for L-lysine synthesis were calculated separately. The 
aspartate DH and dihydrodipicolinate DH carry a flux of 
29.9 mol % (Fig. 3) and the diaminopimelate DH a flux of 
14.3%, adding up to 74.1%. The transaminase flux is 
46 mol % (29.9% aspartate amino transferase, plus 15.6% 
succinylaminotransferase ). It follows that the sum of the 
NADPH-consuming fluxes is 194.9%. which is almost 
identical to the generating fluxes of 210.5%. 

The flux balance for the heterologous mutant is shown in 
Fig. 5B. Here the glutamate dehydrogenase flux does not con- 
tribute to biosynthetic purposes due to the N AD-dependent 
dehydrogenase of P. asaccharolyticus. The NADPH flux of 
dehydrogenases required for biosynthetic purposes is com- 
parable (51.6%). whereas the dehydrogenases operating 
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Bidirectional flux v )6 between 
C r and C 4 -units (%) 

FIG, 6. Degree of "C enrichment in C-2 of oxaloacetaie as a function 
of the exchange flux (exchange between C-4 and C-3 units). The 
dependency is given together with measured values, estimated values, and 
their statistical errors. 

in L-lysine synthesis consume less NADPH, since in this 
particular experiment the total carbon flux toward L-lysine 
was reduced. However, importantly, as the response of the 
reduced NADPH demand due to the NAD-dependent glu- 
tamate DH, the generation of NADPH is drastically 
reduced. The total NADPH generating flux is 138.7 mol%. 
The consuming flux is calculated to be 95.3%, which gives 
an excess of 43.4 %, This flux is exclusively accessible by the 
given isotope analysis and its relevance will be discussed 
below. About one fourth of the reduced NADPH demand in 
the heterologous mutant can be attributed to reduced 
demand for biomass synthesis and three fourths to that for 
L-lysine synthesis. 

DISCUSSION 

Although it is often attempted, it is not possible to predict 
in a necessarily reductionist approach to describing the 



living cell how specific enzyme activities must be preset to 
result in the maximal output of the desired product. In con- 
trast, the real benefit of such an approach is to be found else- 
where (Popper, 1974).This might be in our example the 
quantification of the intracellular flexibility or a deeper 
insight into the exchange fluxes. 

We were able to identify such an exchange flux for a bac- 
terium in vivo within the anaplerotic reactions ( Marx et al y 
1996). The statistical analysis of this exchange flux in 
C.glutamicum shows that v^ (opposed carboxylating and 
decarboxylating reactions) is different in the two glutamate 
DH mutants analyzed. The 90% confidence interval for the 
heterologous mutant is 22.6-36.5% molar exchange flux, 
and the interval for the homologous mutant is 3.8-17.0% 
(Table 6). This difference results from the different levels of 
the l3 C content for the oxaloacetate C-2 position (Tables 4 
and 5). The enrichment in this position responds very sen- 
sitivly to the v^ exchange rate. The corresponding quan- 
titative relation for the homologous mutant is shown in 
Fig. 6. The measured enrichment was 12.7 ± 1.0. Therefore, 
the molar exchange for the individually considered flux of 
v^ is 0-9.3% (relative to glucose uptake). Since the simula- 
tion considers all fluxes and exchanges, the estimated 
exchange between the C-3 units and the C-4 precursors as 
integrated in the total flux system is 10.3%, with a 90% 
confidence interval of 3.8-17.0%. The difference in the 
exchange flux between the two mutants had to be given as 
the sum of the individual fluxes ihterconverting the C-3 
units and the C-4 units, since the partially identical educts 
and products do not yet permit any individual flux quan- 
tification. The situation is even more complex, since only 
recently has it become possible to clone the pyruvate car- 
bo xylase-encoding gene pyc { Peters- Wend isch et al % 1998), 
which is therefore, in addition to phosphoe/w /pyruvate car- 
boxylase (p/>c), the second reaction contributing to the net 
synthesis of oxaloacetate (Fig. 7). As mentioned, the parti- 
tion of the carboxylating flux between the two carboxylase 
reactions identified is not yet known, nor is the physiologi- 
cal role of the backflux at this extensive set of enzymes 
around the anaplerotic reactions. In rats the backflux from 
pyruvate to oxaloacetate was estimated in isotopomer 
studies to be of approximately the same size as that from 
oxaloacetate to phosphoeno/pyruvate (Katz et aL 1993). 
The backflux in C glutamicum may indeed be of fundamen- 
tal physiological significance (Fig. 7). Thus, an energy-con- 
suming futile cycle without contributing to biosynthesis 
( Rognstadt, 1 996} may be formulated via phosphoevw/pyru- 
vate carboxykinase and pyruvate carboxylase. The inhibi- 
tion of both enzymes by ADP corresponds to the idea that 
the ATP- and GTP-consuming cycle is reduced at a low 
energy charge. In glutamate production we found a very 
slight futile cycle flux between C-3 units and C-4 units of 
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FIG. 7. Overview of ihe flux-carrying reactions at the anaplerotic 
node (hat could contribute to the reverse flux of 29.2% in the heterologous 
mutant, including a hypothetical transhydrogenase cycle and an energy- 
consuming cycle. 

only 17.1% (Marx et al, 1997). This could therefore be 
attributed to a low energy charge or increased conservation 
metabolism under the stress of biotin limitation. Another 
cycle which could be constructed consists of the malic enzyme 
and malate dehydrogenase resulting in the transfer of reduc- 
ing equivalents from NAD to NADP. This cycle represents a 
transhydrogenase activity. There is contradictory informa- 
tion on C. glutamicum for a specific transhydrogenase 
(Cocaign-Bousquet and Lindley, 1996; Park et al % 1997; 
Vallino 1991). In general, the sum of the carboxylating 
and decarboxylating reactions apparently serves as a fine 
tuning within the central metabolism to balance the fueling 
reactions with the catabolic reactions (Kornberg, 1966). 
Although the conversions between C-3 units and C-4 units 
seem to be important for the metabolic design of aromatic 
amino acid production in E. coli (Patnaik et al, 1995), it is 
still unclear how they can be selectively used to increase 
amino acid overproduction. It is striking that with high C-4 
unit formation the reverse flux is lowest in the homologous 
mutant of C glutamicum ( 10.3%, see Fig. 3), as is the case 
under glutamate-producing conditions (Marx etal, 1996). It 
is thus conceivable that the net conversion can be promoted 
not only by overexpression of, for example, pyruvate car- 
boxylase, but also by reduction of the reverse fluxes. 

As shown by the stoichiometry in Fig. 5, more NADPH 
is formed than is consumed. However, based on the 90% 
confidence interval, also the greater discrepany calculated 
for the heterologous mutant is still within the error range. 
Thus formation of an NADPH surplus cannot be quantified 
with certainty. Nevertheless, in all flux situations analyzed 
for C. glutamicum to date we have determined an excess of 
NADPH formation and never a deficit (Marx et al, 1997; 
Sonntag et al, 1995). Several additional reactions that are 
difficult to quantify might consume NADPH such as the 



system identified in E. coli detoxifying peroxide radicals 
( Luchi and Weiner, 1 996 ), or an N APDH-oxidizing activity 
which at least has been detected in extracts of C. glutamicum 
(Vallino, 1991 ). Also, a weak reaction of the P. assacharoly- 
ticus glutamate dehydrogenase using NADPH instead of 
NADH might be attributed to the greater discrepancy of the 
heterologous mutant, since the enzyme has a very low 
activity of 0.9% with NADPH as substrate. 

In any case, the cell (Fig. 5) actually responded by 
reduced redox generation through our selective intervention 
of diminished demand for redox equivalents. There are 
several examples of cases where externally added electron 
mediators might affect the intracellular redox metabolism as 
inferred from an altered product spectrum (Neijssel 1977; 
Kwong and Rao, 1992). An internal alteration of the 
NADPH balance was obtained in yeast by the synthesis 
of an NADPH-dependent xylose reductase, which led to 
altered by-product formation (Meinander et aL 1996). Our 
quantitative data show that for the directly altered NADPH 
demand the NADPH generation is reduced by about one- 
third. It thus follows that C. glutamicum is enormously 
flexible not only with respect to the carbon fluxes in the cen- 
tral metabolism (Marx et al, 1997) but also with respect to 
the redox fluxes. This could be a reason for the particular 
suitability of this bacterium for metabolite overproduction. 
Since our data show that the cell adjusts the fluxes to the 
reduced NADPH demand and does not use a potentially 
available increased NADPH supply to increase lysine 
production, we regard this as an indication that lysine 
production is not essentially NADPH-limited. However, 
this question must be addressed in further studies by direct 
intracellular quantification of the involved metabolites. 

In total, we have now generated a unique set of five dif- 
ferent flux situations of isogenic strains specifying the 
response of the metabolism to an altered carbon or reducing 
power flux. Acting only on one enzyme activity, for instance, 
a kinase (Marx et al, 1997) or a dehydrogenase (this work), 
results in an extreme range of fluxes in the central 
metabolism. The oxidative PPP carbon flux varies from 25 
to 76%. As a general trend, it can be seen that at a high 
NADPH demand a high PPP flux results. Furthermore, 
there is an inverse relation of the PPP flux to the TCC flux. 
A global picture of flux regulation in C. glutamicum can thus 
be derived. The TCC flux is first adjusted according to the 
carbon requirements. The isocitrate dehydrogenase activity 
necessarily linked to the latter provides a basal level of 
NADPH flux. In extreme cases this is manifest in glutamate 
production, where a high input into the TCC is present 
due to the high a-ketoglutarate requirement. However, 
there is still high cycling activity due to the a-ketoglutarate 
DH activity present, which was originally thought to be 
absent in C. glutamicum (Kinoshita. 1985). This activity 
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serves to provide sufficient NADH for ATP generation. 
Whereas in the TCC cycle the fluxes are adjusted to the 
specific carbon and NADH demand, the PPP flux is 
adjusted according to the NADPH demand not fulfilled by 
the TCC With high TCC flux (glutamate production) the 
PPP flux is thus reduced, whereas at low TCC flux (lysine 
production) the PPP flux is high to adjust the NADPH 
balance. The specific regulators controlling the entry into 
the PPP are the glucose-6-P DH and 6-P-gluconate DH, 
both controlled in their catalytic acitivity by the NADPH 
concentration via allosteric mechanisms (Sugimoto and 
Shiio. 1987). Due to the flexibility of the central metabolism 
and its design to provide precursors as well as NADPH and 
NADH, it appears that the major targets to be altered to 
obtain maximal productivities are the biosynthesis path- 
ways themselves, together with the export of the metabolite 
to be produced. What might limit maximal product forma- 
tion beyond these reactions cannot be predicted by mecha- 
nistic flux studies, but relates to questions concerning the 
cell in its entirety, for instance, the tolerance of the cell 
toward physical features of the product, or how controlled 
reduction of precursor use for biomass formation can be 
obtained. 

ACKNOWLEDGMENTS 

We lhank E. Bormann. K. Krumbach, and S. Peirovic for help with con- 
structions and analyses, M. Tcsch for sharing results prior to publication, 
and W. Wiechert (University of Siegen). D. Molenaar (University of 
Diisseldorfl, and V. Wendisch (University of California. Berkeley) for 
discussions. 

REFERENCES 

Bern. A. (1996). Improving production of aromatic compounds in 
Escherichia coli by metabolic engineering. TJBTECH 14. 250-256. 

Bormann. £.. Eikmanns. B. J„ and Sahm. H. ( 1992). Molecular analysis of 
the Corynebacterium gluiamicum gdh gene encoding glutamate 
dehydrogenase. Mot. Microbial. 6. 317-326. 

Cocaign-BousqueL M.. Guyonvarch, A., and Lindley. N. D. (19961. 
Growth-rate dependent modulation of carbon flux through central 
metabolism and the kinetic consequences for glucose-limited 
chemostat cultures of Corynebacterium gluiamicum. Appt. Enriron. 
Microbiol. 62, 429-436. 

Eggeling. L.. DeGraaf. A. A., and Sahm, H. ( 1996). Quantifying and direct- 
ing metabolite (lux: Application to amino acid overproduction. Adv. 
Biochem. Eng. 54, 2-30. 

Eggeling. L.. Oberle. S., and Sahm, H. (1997 ). Improved L-lysine yield with 
Corynebacterium gluiamicum: Use ofdapA resulting in increased flux 
combined with growth limitation. Appi Microbiol. Bioiechnol. 49, 
24-30. 

Eggeling. L„ Morbach, S., and Sahm, H. (1997). The fruits of molecular 
physiology: Engineering the L-isoleucine biosynthesis pathway in 
Corynebacterium gluiamicum. J. Bioiechnol. 56, 167-182. 

Eggeling. L.. and Sahm. H. (1998). Amino acid production: Principles of 
metabolic engineering. In "Metabolic Engineering" (Y. S. Sang and 
E. T. Papoutsakis. Eds.), Biotechnology Intelligence Unit. R. G. 
Landes. 



mi-taw >uc i:nc;ini:i.rin(J I. 35-48 ( 1999) 

ARTK'LI; NO. MT9801O6 

Eikmanns. B.. Kleinertz. E.. Liebl. W.. and Sahm. H. (1991). A family of 
Corynebacwrium gluiamicum/ Escherichia coli shuttle vectors for gene 
cloning, controlled gene expression, and promoter probing. Gene 102, 
93-98. 

Flores. N.. Xiao. J.. Berry, A., Bolivar, F.. and Valle. F. ( 1996). Pathway 

engineering for the production ofaromatic compounds in Escherichia 

coli. Nature Biotcchnot. 14, 620-623. 
luchi, S., and Weiner, L. (1996). Cellular and molecular physiology of 

Escherichia coli in the adaption to aerobic environments. J. Biochem. 

120. 1055-1063. 

Katz, J., Wals, P., and Lee, W. P. (1993). Isotopomcr studies orglucone- 
ogenesis and the Krebs cycle with ,J C-labeIed lactate. J. Biol. Client. 
268, 25509-25521. 

Kinoshita, S. (1985). Glutamic acid bacteria. In "Biology of Industrial 
Microorganisms" (A. L. Demain and N. A. Solomon, Eds.), 
pp. 1 15-142, Benjamin/Cummings, Redwood City, CA. 

Koniberg, H. L. (1966). Anaplerotic sequences and their role in meta- 
bolism. In "Essays in Biochemistry" (P. N. Campbell and G. D. 
G revile. Eds), Vol II, pp. 1-31, Academic Press, New York. 

Kwong, S. C. W„ and Rao, G. (1992). Effect of reducing agents in an 
aerobic amino acid fermentation. Bioiechnol. Bioeng. 40. 851-857. 

Liebl, W., Bayerl, A., Stillner. U., and Schleirer. K. H. (1989). High 
efficiency elect roporation of intact Corynebacterium gluiamicum cells 
FEMS Microbiol. Lett. 65, 299-304. 

Marx. A., de Graaf, A. A., Wiechert, W., Eggeling. L.. and Sahm, H. ( 1996). 
Determination of the fluxes in the central metabolism ofCoryncbac- 
terium gluiamicum by NMR spectroscopy combined with metabolite 
balancing. Bioiechnol. Bioeng. 49, 1 1 1-129. 

Marx. A., Striegel. K„ de Graaf, A. A., Sahm, H., and Eggeling, L. ( 1997). 
Response of the central metabolism of Corynebacterium gluiamicum 
to diflerent flux burdens. Bioiechnol. Bioeng. 56, 168-180. 

Meinander. N.. Zacchi, G.. and Hah-Hagerdal. B. (1996). A heterologous 
reductase affects the redox balance of recombinant Saccharomyces 
cerevisiae. Microbiology 142. 165-172. 

Mori, M.. and Shiio, I. (1985). Purification and some properties of 
phosphoenolpyruvate carboxylase from Brevibavterium flavum and 
its aspanate-overproducing mutant. J. Biochem. Tokyo 97. 1119- 
1128. 

Mori, M., and Shiio, 1. (1985). Synergistic inhibition of phosphoenol- 
pyruvate carboxylase by aspartate and 2-oxoglutarate in Brecihuc- 
terium Jlarum. J. Biochem. Tokyo 98, 1621-1630. 

Neese, R. A., and Helleretein, M. K. (1995). Calculations for gluconeoge- 
nesis by mass isotopomer distribution analysis. / Biol. Chem. 276. 
14464-14466. 

Neijssel, O. M. (1977). The efTect of 2.4-dinitrophenol on the growth of 
Klebsiella uerogenes NCTC 41 8 in aerobic chemostat cultures. FEMS 
Lett. 1.47-50. 

Oshima. K., Tanaka. K... and Kinoshita, S. (1964). Studies on L-glutamic 
acid fermentation. XI. Purification and properties of L-glutamic 
dehydrogenase from Micrococcus gluiamicus. Agrk. Biol. Chem. 
Tokyo 28. 714-722. 

Park. S. Sinskey, A. J., and Stephanopoulos, G. (1997). Metabolic and 
physiological studies of Corynebacterium gluiamicum mutants. 
Bioiechnol. Bioeng. 55, 864-879. 

Patnaik, R., Spiizer. R. G., and Liao, J. C. ( 1995 ). Pathway engineering for 
production of aromatics in Escherichia coli: Confirmation of 
stoichiometric analysis by independent modulation of AroG. TktA. 
and Pps activities. Bioiechnol. Bitteng. 46. 361-370. 

Peters- Wendisch. P.. Kreutzer, C. Kalinowski. J., Patek. M., Sahm. H.. 
and Eikmanns. B. J. (1998). Pyruvate carboxylase from Coryncbac- 
lerium gluiamicum : Characterization, expression and inactivation of 
the pyc gene. Microbiol. UK 143. 1095-1 103. 



47 



metabolic i:n(;ini:i:rinc; I. 35-48 ( 1 999 > 
ARTICLE NO. MT980I06 



Marx et al. 



Peters-Wendisch. P. C. Wcndisch. V. F.. Paul. S„ Eikmanns, B. J., 
and Sahm. H. (1997). Pyruvate carboxylase as an anaplerotic 
enzyme in Corynehacterium glutamicum. Microbiol. UK 143. 1095- 
1 103. 

Popper. K. ( 1974). "Studies on the Philosophy of Biology." (F.J, Aydaand 
T. Dobzhansky, Eds.), pp. 259-283, MacMillan, London. 

Rognstad. R. (1996). Futile cycling in carbohydrate metabolism. I. Back- 
ground and current controversies on pyruvate cycling. Biochem. Arch. 
12. 71-83. 

Sauer. U„ Haizimanikatis. V., Hohmann. H„ Manneberg, M„ van Loon, A., 
and Bailey (1996). Physiology and metabolic fluxes of wild-type and 
riboflavin-producing Bacillus suhtilis. Appl. Environ. Microbiol. 62. 
3687-3696. 

Schafer. A., Kalinowski. J„ Simon, R., Seep-Feldhaus, A., and Puhler, A. 
( 1990). High-frequency conjugal plasmid transfer from gram-negative 
Escherichia coli to various gram-positive cor/neform bacteria. J. Bac- 
terial. 172. 1663-1666. 

Schafer, A., Tauch, A. Jager, W., Kalinoski, J., Thierbach, G„ and 
Puhler. A. (1994). Small mobilizable multi-purpose cloning vectors 
derived from the Escherichia coli plasmids pK18 and pK19: Selec- 
tion of defined deletions in the chromosome of Corynehacterium 
glutamicum. Gene 145, 69-73. 

Schmmpf, B., Eggeling, L.. and Sahm, H. ( 1992). Isolation and prominent 
characteristics of an L-lysine hyperproducing strain of Corynehac- 
terium glutamicum. Appl. Microbiol. Biotechnol. 37, 566-571. 

Schwarzer, A., and Puhler, A. ( 1990). Genetic manipulation of the amino 
acid-producing Corynehacterium glutamicum strain ATCC 1 3032 by 
gene disruption and gene replacement. Bio/Technology 9, 84-87. 

Simon. R.. Priefer, U., and Puhler, A. ( 1 983 ). A broad host range mobiliza- 
tion system for in vivo genetic engineering: Transposon mutagenesis 
in gram negative bacteria. Bio/Technology 1, 784-791. 

Snedecor. B.. Chu, H., and Chen, E. (19911. Selection, expression, and 



nucleotide sequencing of the gluiamale dehydrogenase gene of Pep- 
tostrcplovoccus asaccharolyticus. J. Bacteriol 173.6162-6167. 

Sonntag, K„ Schwinde. J., de Graaf. A. A.. Marx, A.. Eikmanns. B. J., 
Wiechert, W.. and Sahm. H. ( 1995). 1 3-C NMR studies of the fluxes 
in the central metabolism of Corynehacterium glutamicum during 
growth and overproduction of amino acids in batch cultures. Appl. 
Microbiol. Biotechnol. 44, 489-495. 

Sugimoto, S.. and Shiio, I. (1987). Regulation of gIucose-6-phosphate 
dehydrogenase in Brerihacierium flavum. Agric. Biol. Chcm. Tokyo 
51, 101-108. 

Vallino, J. J. (1991). "Identification of Branch-Point Restrictions in 
Microbial Metabolism Through Metabolic flux Analysis and Local 
Network Perturbations." Ph.D. thesis. Massachusetts Insitute of 
Technology, Boston. 

Vallino, J. J., and Stephanopoulos. G. ( 1993). Metabolic flux distributions 
in Corynehacterium glutamicum during growth and lysine over- 
production. Biotechnol. Bioeng. 41, 633-646. 

Vrljic, M.. Sahm, H., and Eggeling. L. ( 1996). A new type of transporter 
with a new type of cellular function: L-Lysine export in Corynehac- 
terium glutamicum. Mol. Microbiol 22. 815-826. 

Weisbrod, R. E., and Meister. A. (1973). Studies on glutamine synthetase 
from Escherichia coli. J. Biol. Chem. 248. 3997-4001 

Wendisch, V. F„ de Graaf. A. A., and Sahm, H. ( 1997). Accurate deter- 
mination of 1J C enrichments in nonproionated carbon atoms of 
isotopically enriched amino acids by 1 H nuclear magnetic resonance. 
Anal Biochem. 245, 196-202. 

Wiechert, W., Siefke, C de Graaf, A. A., and Marx. A. (1997). Bidirec- 
tional reaction steps in metabolic engineering. II. Flux estimation and 
statistical analysis. Biotechnol Bioeng. 55, 1 18-135. 

Wiechert. W., and de Graaf, A. A. ( 1997). Bidirectional reaction steps in 
metabolic engineering; Modeling and simulation of carbon isotope 
labeling experiments. Biotechnol Bioeng. 55. 101-117. 



48 



Metabolic Engineering 



Gregory N. Siephanopoulos 

Massachusct: • institute of Technology 
Cambridge, Massachusetts 



EDITORS 
Anthony J. Sinskey 

Massachusetts Institute of Technology 
Cambridge. Massachusetts 



Martin L, Yarmush 

Massachusetts General Hospital 
Boston, Massachusetts 



ASSOCIATE EDITORS 
Barry C. Buckland Hans Westerhoff 

Meivk and Company Free University 

Rahwaw New Jersey Amsterdam, The Netherlands 



Douglas Cameron 

University of Wtuttnsin 
\ ft tdison. Wisconsin 

Arnold L. Demain 

Massat Itusetts Institute of Technology 
Cambridge, Massachusetts 

Peter Dunnill 

University College London 
London. United Kingdom 

Jean Marc Engasser 

tnslitut Sationat ftilytcchnidue 

dt Lorraine 
Vandoeuvre les Nancy Cedes. Fram e 

David Fell 

Oxford Hnntl.es University 
Oxfttrd. United kingdom 

Peter Gray 

University of Ken South Wales 
Sydney. Australia 

J. J. Heijnen 

Technology University of Delft 
Delft. The Netherlands 

Lonnie Ingram 

University of Florida 
Gainesville. Florida 

Thomas Jeffries 

USDA Forest Sen-ice 
Madison. Wim initio 



EDITORIAL BOARD 
Joanne Kelteher 

Grorge Washington University 

Medical Center 
Washington. DC 

Chaitan S. Khosla 

Stanford University 
Stanford. California 

Susan Leschine 

University of Massachusetts ai Amherst 
Amherst. Massachusetts 

James Liao 

University of California 
Lf»v Angeles, California 

Jens Nielsen 

Denmark Technical I 'nirersity 
Lynghy. Denmark 

Bernhard O. Palsson 

University of California. San Diego 
La Jolla. California 

E. Terry Papoutsakis 

Sarthwestcnt University 
F.vanston, Illinois 

Steve Picattagio 

DttPaiu 

Witiniitgton, Delaware 



Matthias Reuss 

Unhersiiai Stuttgart 
Siangan. Germany 

Herman Sahm 

tnslitut Jiir Hint ethnologic dcr KF1 
Jidicit. Genuany 

Ka-Yiu San 

Rice University • 
Houston. Texas 

Kenneth N. Timmis 

GRF-Xationul Research Cent w for 

Biotechnology 
Brunnsctnt eig. Gei nmny 

Ronald Tompkins 

Massachusetts General Hospital 
Huston. A tussat husetts 

Donald Trimbur 

(ienccor International 
Pah Alto. California 

J. Craig Venter 

Institute for Genomic Research 
Rot kvillr. Maryland 

Chris Walsh 

Harvard Medical School 
Hosttnt. Massat husetis 

Vernon R. Young 

Massat husetts Inslitutt ■ of let -tun tfogy 
Cambridge. Massat itusetts 



Cover figure. Central an;<erohic metabolic pathway <»f £'. coli. The iWnicftUilioti end-products include acetoin pnxluccd hy cells containing ihc ALS plasmid. The llu\es through each 
pathwax arc designated t>, through u i: . (For detail*. mX the paper hy Vang et at. ( IW*>. Mctab. Fug. 1. 26- U.) 



The Journal or Biological Chemistry 
Vol. 257, No. 3, Issue of February 10, pp. 1189-1195, 1982 
Printed in U.S. A. 

13 C Nuclear Magnetic Resonance Studies of the Biosynthesis by 
Microbacterium ammoniaphilum of L-Glutamate Selectively Enriched 
with Carbon-13* 

(Received for publication, July 13, 1981) 

T. E. Walker, C. H. Han, V. H. Kollman, R. E. London, and N. A. Matwiyoff$ 

From the Los Alamos National Laboratory, University of California, Los Alamos, New Mexico 87545 



Doc. Ref. AS3 
Appl.No. To be assigned 



13 C NMR of isotopicaUy enriched metabolites has 
been used to study the metabolism of Microbacterium 
ammoniaphilum, a bacterium which excretes large 
quantities of L-glutamic acid into the medium. Biosyn- 
thesis from 90% [l- 13 C]glucose results in relatively high 
specificity of the label, with [2,4- 13 C 2 ]glutamate as the 
major product. The predominant biosynthetic pathway 
for synthesis of glutamate from glucose was deter- 
mined to be the Embdcn Meyerhof glycolytic pathway 
followed by P-enolpyruvate carboxylase and the first 
third of the Krebs cycle. Different metabolic pathways 
are associated with different correlations in the enrich- 
ment of the carbons, reflected in the spectrum as dif- 
ferent "C- l3 C scalar multiplet intensities. Hence, inten- 
sity and 13 C- l3 C multiplet analysis allows quantitation 
of the pathways involved. Although blockage of the 
Krebs cycle at the a-ketoglutarate dehydrogenase step 
is the basis for the accumulation of glutamate, signifi- 
cant Krebs cycle activity was found in glucose grown 
cells, and extensive Krebs cycle activity in cells metab- 
olizing [l- 13 C]acetate. In addition to the observation of 
the expected metabolites, the disaccharide a,a-treha- 
lose and a,j8-glucosylamine were identified from the 13 C 
NMR spectra. 



Several varieties of bacteria have been described which 
utilize glucose or acetate as a sole carbon source and which 
excrete large quantities of amino acids (I). One of these is 
Microbacterium ammoniaphilum, which has an impaired 
Krebs cycle due to a reduced activity of the enzyme(s) nec- 
essary for processing a-ketoglutarate and its metabolites. It 
can be induced to excrete L-glutamate late in the growth cycle 
of cells cultured in media of high nitrogen content at suffi- 
ciently low concentrations of biotin that the cells develop 
altered membranes which allow facile transport of L-gluta- 
mate formed by transamination of a-ketoglutarate (2-4). We 
concluded that these characteristics of M. ammoniaphilum 
made it an excellent potential vehicle for the biosynthesis of 
L-glutamate selectively labeled with 13 C, if it were grown on 

* This work was performed under the auspices of the United States 
Department of Energy. The proton NMR spectra were obtained at 
the Colorado State University Regional NMR Center, funded by the 
National Science Foundation Grant CHE-78-18581. The costs of 
publication of this article were defrayed in part by the payment of 
page charges. This article must therefore be hereby marked "adver- 
tisement" in accordance with 18 U.S.C. Section 1734 solely to indicate 
this fact. 

% Recipient of National Institutes of Health Research Grant RR- 
00962-05 from the Division of Research Resources, Department of 
Health, Education and Welfare. 



specifically enriched glucose or acetate as the predominant 
source of carbon. 

Our interest in these systems developed in the continuing 
program at the Los Alamos National Laboratory on the use 
of microorganisms and cell-free extracts for the large scale 
synthesis of natural products, uniformly or specifically labeled 
with 13 C (5). Current work is focused on those l- amino acids 
for which efficient organic synthesis methods are not now 
available and for which there is a need either in human 
metabolic and nutritional studies or in the investigation of the 
structure and dynamics of proteins enriched with labeled 
amino acids. In the use of biological systems for the selective 
enrichment of natural products, there are two prime practical 
considerations: 1) the optimization of the incorporation of the 
l3 C label from the precursor into the product; and 2) the 
mmimization of the degree to which the label becomes ran- 
domized in the product due to the flux of the precursor and 
its metabolites through various pools. Especially at high levels 
of enrichment of products for mass spectrometric and nuclear 
magnetic resonance studies, randomization of the label leads 
to complications arising from multiple isotope clusters and 13 C 
nuclear spin coupling patterns (6-8). 

In this paper we address the latter consideration in a I3 C 
NMR analysis of the flow of 13 C label through the Krebs cycle, 
the phosphogluconate pathway, and the glyoxylate shunt into 
selectively labeled L-glutamate excreted by M. ammonia- 
philum grown in media containing [l- l3 C]glucose or [1- 13 C] 
acetate as the predominant carbon source. The work illus- 
trates the general utility of selective, in contrast to random, 
13 C labeling of natural products biosynthetically. In addition 
it is of interest because it demonstrates that it is not generally 
valid to assume in biosynthetic labeling experiments that the 
label from 13 C-enriched precursors will occupy the sites of a 
natural product in a statistically independent manner. 

EXPERIMENTAL PROCEDURES 

Materials— I l- 13 C]Glucose enriched to 90 atom % (9) and [1- 13 C] 
acetate (10) enriched to 70 atom % were gifts of V. Kerr and T. 
Sanchez, respectively, both of the Los Alamos National Laboratory. 
Corn steep liquor was a gift of Grain Processing Company of Musca- 
tine, IA. Cultures of M. ammoniaphilum (ATCC 15354) were obtained 
from the American Type Culture Collection, Rockville, MD. 

Nuclear Magnetic Resonance Spectroscopy — Pulsed ,3 C NMR 
spectra were obtained with a Varian XL- 100- 15 spectrometer (25.2 
MHz) interfaced to a Data General SuperNova computer. The spectra 
were recorded at 25 °C with a spectral width of 5000 Hz and 4096 
spectral points. The lock signal was obtained from -5% D 2 0 contained 
in the sample. Spectra of the production media were obtained on the 
supernatant after removing the cells by centrifugation at 12,000 X g 
for 10 m. The supernatant could then be run directly, appropriately 
diluted, or concentrated by reducing the pH and evaporating in 
vacuo. Standard samples were dissolved in 0.05 M potassium phos- 
phate buffer and the pH adjusted to 7.5. Proton NMR spectra for 
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DiO solutions were obtained on a Nicolet spectrometer in the FT 
mode operating at 360 MHz. "C chemical shifts were measured with 
respect to external tetramethyl silane and 'H shifts with respect to 
the methyl resonance of internal trimethylsilylpropanesulfonic acid. 

l3 C NMR Spectroscopy of Cell Suspensions— M. ammoniaphilum 
was grown on a slant of nutrient agar at 30 °C for 3-4 days. One 
loopfull of the microorganism was precultured in 30 ml of a rich 
medium containing 5% corn steep liquor, 5% casein, 0.1% MgSCv 
7H 2 0 and 0.1% KH 2 P0 4 at pH 7.3 for 24 h at 30 °C in a 500-ml 
Erlenmeyer flask in a shaker. The fermentation was carried out at 
30 °C in a baffled 500-ml Erlenmeyer flask containing 30 or 60 ml in 
a production medium (pH 7.3) which contained 5% glucose, 1.2% urea, 

0. 5% KHzPO,, 0.1% MgSO<7H 2 0, 0.001% FeS(V7H 2 0, 0.001% 
MnSCv4H 2 0, biotin (1.5 /ig/liter), and a 2% innoculum of the pre- 
culture medium. 

M. ammoniaphilum cells were grown for 24 h at 300 rpm in an 
orbital shaker, harvested by centrifugation at 12,000 x g for 10 min, 
washed twice with 0.2% KC1, and suspended in fresh media containing 
[l- 13 C]glucose or [l- l3 C]acetate. Aeration was provided by continu- 
ously bubbling a stream of pure oxygen gas through a very thin 
capillary tube which reached to within 1-2 mm of the bottom of the 
spinning 12-mm NMR tube. The total volume in the NMR tube was 
4 ml to prevent bubbles breaking the surface from disturbing the 
magnetic field, and to allow for evaporation of sample due to oxygen 
bubbling. Other details are given in the figure captions. Spectra were 
continually obtained for 2-3 h increments over a 20-h period, at which 
time all of the glucose had been consumed and the glutamate level 
was reasonably stable. 

RESULTS 

,3 C and l H NMR Spectra— The time dependence of the I3 C 
NMR spectra of a suspension of M. ammoniaphilum contain- 
ing [l- 13 C]g!ucose near the stationary phase of growth is 
summarized in Fig. 1. Since the production of glutamate in 
large quantities by this microorganism does not begin until 
that stage is attained (1), the initial growth was accomplished 
with a natural abundance of glucose to rruriimize loss of the 
13 C label. In addition to the expected resonances from 13 C- 
enriched bicarbonate, glutamate, succinate, and lactate, the 
spectra in Fig. 1 exhibit prominent peaks from two unusual 
l3 C-labeled products which are eventually consumed: treha-. 
lose and glucosylamine. The accumulation of trehalose, which 
presumably functions as a storage carbohydrate, has been 
observed previously by 13 C NMR in growing yeast (11, 12), 
bacteria (13), and differentiating amoeba cultures (14). To our 
knowledge, glucosylamine formation in cell culture has not 
been observed and, at present, we do not know whether its 
synthesis is under enzyme control and what role it plays in 
the control of glucose or nitrogen metabolism. 1 The formation 
of lactate and succinate by M. ammoniaphilum depends on 
the oxygen tension which, for the cultures appropriate to Fig. 

1, fell sharply in the later stages of glucose consumption. In 
fully aerated cultures, the lactate and succinate levels are 
sharply reduced (Fig. 2). The time dependence of the l3 C 
NMR spectra of the cell-free production media (see "Experi- 
mental Procedures") obtained under fully aerobic and par- 
tially anaerobic conditions was quantitatively like the spectra 
of the cell suspensions. This observation, together with the 
lack of ,3 C-enriched resonances in the spectra of medium-free 
packed cells obtained at the end of the metabolism experi- 
ments, suggests that the I3 C spectra monitored (Figs. 1 and 2) 
were those of extracellular substrates and metabolites. A 
qualitative appreciation of the flow of labeled carbon in this 
system can be obtained from the summary of the time de- 
pendence of the intensities of the 13 C resonances of the glucose 
substrate and the 13 C-labeled metabolites tabulated in Fig. 3 

1 The appearance of glucosylamine in the 13 C NMR spectra of 
metabolizing cells is not inconsistent with the ease of its formation 
under a variety of conditions (15-17). In addition, cell-free extracts of 
M. ammoniaphilum will catalyze the formation of aj?-[l- 13 C, ls N] 
glucosylamine from [l- ,3 C]glucose and [ ,5 N2]urea. 
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Fig. 1. Time dependence of the proton-decoupled Fourier 
transform 1S C NMR spectra (25.2 MHz) of a suspension of M. 
ammoniaphilum initially containing [1-"C (90 atom %)]glu- 
cose. T ~ time interval for accumulation of the spectra after transfer 
of the cells to fresh medium containing [l- I3 C]glucose (see "Experi- 
mental Procedures"). Enriched ,3 C resonances are: HC0 3 ~ ion; Glc, 
a and J? anomers of [l- ,3 C]glucose; T, C-l of.a,a-[l,l'- 13 C2]trehalose; 
Git, C-2 (56.0 ppm), C-4 (34.6 ppm), and C-3 (28.3 ppm) of glutamate; 
S, C-2 of succinate; L, C-3 of lactate; and GA, jff- and a- C-l resonances 
of glucosylamine. In the T = 0-2 h and 7-9 h spectra, the glucose C- 
1 resonances are truncated and the natural abundance glucose reso- 
nances are apparent. 



for a particular experiment. The yield of glutamate from 
labeled glucose in experiments like these typically is 35-40%. 
The l3 C chemical shifts of the metabolites are tabulated in 
Table I. 

The 360 MHz proton NMR spectrum of a solution of 
[l,2,3,4- 13 C 4 ]glutarnate derived from M. ammoniaphilum 
grown on [l- l3 C]glucose is shown in Fig. 4A. Despite the high 
level of 13 C incorporation demonstrated by the proton spec- 
trum, the C-2 and C-4 resonances in the 13 C spectrum are 
predominantly singlets (Fig. 4B), suggesting that the 13 C label 
is introduced into specific sites (C-2, C-4) with a minimum of 
scrambling to other sites. This site specificity is even greater 
in the early production phase of glutamate. A plot of the 
intensity ratio of C-3 to C-4 versus time is presented in Fig. 5, 
which shows that C-3 is more highly labeled at later times. 
An even higher degree of site-specific incorporation is illus- 
trated in Fig. 6, which is a reproduction of the 13 C NMR 
spectrum of a solution of [l,5- 13 C2]glutamate obtained from 
M. ammoniaphilum grown on [l- I3 C]acetate. 

Major Metabolic Pathways in M. ammoniaphilum— The 
metabolic pathways of M. ammoniaphilum have been estab- 
lished primarily through "C labeling studies carried out by 
Shiio and co-workers (Fig. 7). A glucose carbon source is 
metabolized primarily via the Embden Meyerhof pathway to 
phosphoenolpyruvate and pyruvate (18, 19). Accumulation of 
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Fig. 2. Time dependence of the proton-decoupled "C NMR 
spectra of a fully aerated suspension of M. ammoniaphilum 
grown on [1-"C (90 atom %)]glucose. See Fig. 1 legend. 
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Fig. 3. Time dependence of the intensities of the "C reso- 
nances of the enriched glucose substrate and the metabolites 
identified in Fig. 1. See text also. 
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Fig. 4. Proton NMR and proton-decoupled "C NMR spectra 
of [l,2,3,4-"C 4 ]glutamate. a, proton NMR spectrum (360 MHz) of 
[1 ,2,3,4- 13 C4]glutamate derived from M. ammoniaphilum grown on 
[1- ,3 C (90 atom %)]glucose. For each proton, the center of the multi- 
plet arises from 12 C-H moieties with fine structure caused by H-H 
scalar interactions. The doublets with similar fine structure are due 
to the 13 C-H splitting from moieties containing the ,3 C label. The 13 C 
populations calculated from the ratio of the doublet to the singlet 
intensities are: C-2, 34%; C-4, 39%; and C-3, 14%. b, proton-decoupled 
"C NMR spectrum (25.2 MHz) of [l,2,3,4- l3 C 4 ]glutamate derived 
from M, ammoniaphilum grown on [l- ia C (90 atom %)]glucose. The 
spectrum illustrates the nonrandom distribution of the "C label 
among the C-2, C-3, and C-4 sites. If the label were distributed 
randomly with the 13 C abundances derived from the spectra in B, 
then, for example, the C-l and C-5 signals should consist of approxi- 
mate 1:2:1 multiplets. 

glutamate is due primarily to the reduced activity of a-keto- 
glutarate dehydrogenase (20) which, in effect, shuts down the 
Krebs cycle. As a result, a-ketoglutarate levels are elevated 
and, in the presence of sufficient ammonia, glutamate levels 
will consequently become elevated. Excretion of glutamate 
into the medium is found to be dependent on growth on 
limiting biotin levels (2, 3); the effect of the biotin is related to 
the production of faulty, "leaky" membranes which permit 
this excretion (21). 

As is apparent from Fig. 7, several pathways connecting 
glucose to glutamate are available to the organism. The short- 
est pathway for direct conversion follows the glycolytic deg- 
radation of glucose to pyruvate and phosphoenolpyruvate, 
formation of acetyl CoA and oxaloacetate via decarboxylation 
and carboxylation, respectively, and the first third of the 
Krebs cycle which combines the oxaloacetate and acetyl CoA 
to yield citrate and, subsequently, a-ketoglutarate. On the 
basis of this pathway, the glutamate I3 C-labeling pattern 
resulting from any initial glucose l3 C-labeling pattern is readily 
determined: [l- 13 C]glucose leads to the synthesis of [3- 13 C] 
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Table I 

"C NMR chemical shifts of metabolites 
Samples were dissolved in water or 0.05 m phosphate buffer, pH 
-7 to 7.5, 25 °C. 



GLUCOSE 



Compound 






Chemical shift 






C-l 


C-2 


C-3 


C-4 


C-5 


C-6 


a-D-Glucose 


93.2 


72.5 


73.9 


70.7 


72.5 


61.7 


0-D-Glucose 


97.0 


75.3 


76.9 


70.7 


76.9 


61.8 


a,a-Trehalose° 


94.2 


73.2 


73.7 


70.8 


72.2 


61.7 


a-D-Glucosylamine 6 


84.2 












^-D-Glucosylamine* 


86.2 












L-Glutamate 


175.8 


55.8 


28.3 


34.6 


182.2 




L-Glutamine 


175.0 


55.4 


27.5 


32.0 


178.5 




L-Lactate 


183.4 


69.7 


2L3 








Succinate 


183.3 


35.3 


35.3 


183.3 







° Assignments are from Ref. 14. 

* Shifts for ajS-glucosylamine C-2 through C-6 were not observed. 
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Fig. 5. Time-dependent ratio of C-3 to C-4 intensities in ["C] 
glutamate. Data were taken from the ,3 C NMR spectra of the 
experiment depicted in Fig. 2. 
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Fig. 6. Proton-decoupled "C NMR (25.2 MHz) spectrum of 
[l,5- 1 *Cj]glutamate obtained from M. ammoniaphilum grown 
on [1-"C (70 atom %)]acetate. The cells were grown on natural 
abundance glucose for 24 h and then transferred to the medium 
containing [l- 13 C]acetate and cultured for an additional 96 h. Only C- 
1 and C-5 are labeled. From the multiplet to singlet intensities of C- 
4, the 13 C enrichment at C-5 is calculated to be 70 atom %. From the 
corresponding intensities for the C-2 resonance, the l3 C enrichment 
at C-l is calculated to be 35 atom %. 
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Fig. 7. Major metabolic pathways in Af. ammoniaphilum 
involving glucose, acetate, and glutamate. Glucose labeled at C- 
1 produces [3- ,3 C]pyruvate via the Embden-Meyerhof pathway 
(EMP) and unlabeled pyruvate via the hexose monophosphate shunt 
(HMS). [3- 13 C]Pyruvate enters the tricarboxylic acid {TCA) and 
glyoxylate shunt (GS) cycles as [3- 13 C]oxaloacetate and/or [2- 13 C] 
acetate, and can result in the formation of [2- ,3 C]glutamate, [4- i3 C] 
glutamate, and [2,4- 13 C]glutamate via a-ketoglutarate formed in a 
third of a turn of the TCA cycle. Formation of glutamate after one or 
more turns of the TCA cycle will tend to randomize the label because 
of the formation of the symmetrical intermediates, succinate and 
fumarate. Similar scrambling will occur through the flow of citrate 
through the glyoxylate cycle. 

pyruvate and [2- 13 C] acetate and, consequently, to [2,4- 13 C 2 ]- 
labeled glutamate. Reference to Fig. 4 for the glutamate 
released into the medium indicates a pattern which is sub- 
stantially in agreement with this prediction, supporting this 
as the predominant route of glutamate synthesis from glucose. 

As is apparent from Fig. 7, the pathway summarized above 
is not the only route of glutamate synthesis, and the impor- 
tance of other pathways is suggested by deviations from the 
above pattern, e.g. labeling of C-l and C-3 of glutamate, which 
also occurs to a limited extent. In the following section, we 
consider the general approach of estimating the contributions 
of alternative pathways from the available NMR data. 

Quantitative NMR Analysis of Metabolic Pathways — 
Quantitation of the 13 C population at specific sites in multiply 
labeled molecules requires the accurate evaluation of the 
ratios of singlet-to-multiplet intensities in both 'H and 13 C 
NMR spectra. This is usually a straightforward task, but it 
must be recognized that the satellite resonances can exhibit 
different Ti and/or nuclear Overhauser effect values than the 
center resonance. The differential T! effect in *H NMR spectra 
can be quite significant. For example, a center band ( l2 C-H) 
to satellite ( 13 C-H) l H T! ratio greater than 4 has been 
observed for a l3 C-enriched sugar methine proton (22) . In such 
cases, the pulse delay must be sufficiently long that the center 
band is not overpulsed. A similar, but smaller effect occurs in 
the 13 C NMR spectra of enriched molecules. Although the 
13 C- 13 C dipolar interaction is negligible, relative to the 13 C-H 
dipolar interaction for proton bearing carbon atoms (23), the 
former interaction can be large for quaternary carbons and 
lead to significant differences in the satellite ( l3 C- l3 C) and 
center band ( 13 C- 12 C) Ti values. Again, if the pulse delay is 
not long enough, the intensity of the center band will be 
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underestimated due to overpulsing. Another related potential 
source of error in the analysis of 13 C multiplet intensities of 
quaternary carbon atoms is a reduction in the nuclear Over- 
hauser effect value due to 13 C- 13 C dipolar interactions and 
which leads to a reduction in the satellite ( 13 C- 13 C)/center 
band ( I3 C- 12 C) intensity ratio (7). In the evaluation of the 
relative intensities of the multiplets discussed in the following, 
these potential problems have been avoided. 

The satellite-to-center band intensities, in the 'H NMR 
spectrum (Fig. 4a) of the L-[l,2,3,4- 13 C4]glutamate derived 
from M. ammoniaphilum grown on [l- 13 C]glucose, directly 
give the total I3 C enrichment of C-2 (34 atom %), C-4 {39 atom 
%), and C-3 (14 atom %). Although the relative total integrated 
intensities of the 13 C resonances of C-2, C-3, and C-4 are 
consistent with these enrichments, this high level of enrich- 
ment is not reflected in the 13 C- I3 C multiplets in the 13 C NMR 
spectrum (Fig. 4 b), as noted previously. The apparent dis- 
crepancy arises from the presence of 13 C Isotopomers (e.g. 
I3 C 4 12 C 3 13 C2 and 12 C 4 13 C 3 I3 C 2 ), in which there is a high degree 
of negative correlation in the 13 C labeling that does not affect 
the proton spectrum, but has a major influence on the ap- 
pearance of l3 C multiplets. For example, the C-2 enrichment 
obtained from the H-2 satellites is 34 atom %, whereas the 
apparent enrichment calculated from the C-l satellites (Fig. 
46) is only 25 atom %. This apparent discrepancy can be 
explained by a negatively correlated labeling of C-l and C-2: 
[3- 13 C]oxaloacetate forms [2- 13 C]glutamate via a ^ turn of the 
Krebs cycle, whereas in a single turn of the Krebs cycle, [4- 
13 C]oxaloacetate is formed which can be the precursor for [1- 
13 C]glutamate. Note that these reactions will not introduce 
label into C-l and C-2 in the same molecule. 

Having the absolute 13 C enrichments from the l H NMR 
spectra and the multiplet intensities from the 13 C NMR spec- 
tra, one can use several approaches to interpret the data 
quantitatively, on the basis of the metabolic pathways oper- 
ative. In the present study, we have used the method of 
mixtures with the essential elements summarized below. 

1) Variables are assigned to each alternative path with the 
sum normalized to 1.0. 2) a labeling pattern is associated with 
each path, based on the biochemical pathways expected. 3) a 
table is formulated in which the various possible labeled 
species resulting from each pathway are included. 4) from the 
above table, a second table giving the carbon enrichments and 
multiplet probabilities can be constructed and then compared 
with the data for quantitation. 

The mixture analysis is most simply illustrated for the 
quantitation of the two competing pathways for glucose deg- 
radation. 1) the dominant Embden Meyerhof pathway, and 2) 
the hexose monophosphate shunt. We assign a probability / 
to the former pathway and 1 — / to the latter. If the labeling 
probability of the starting material [l- I3 C]glucose is desig- 
nated on, the labeling probability of pyruvate C-3 or the acetyl- 
CoA produced by pyruvate decarboxylation will be a,/2 due 
to the formation of two trioses via this pathway. The labeling 
of acetyl-CoA via pathway (2) will be zero since the label is 
removed as C0 2 . Although the labelings of pyruvate or acetyl- 
CoA derived from pyruvate are not directly observed, analysis 
of the various biosynthetic pathways indicates that regardless 
of the particular path, glutamate C-4 is always labeled from 
C-2 of acetyl- Co A and so provides a direct measure of this 
pool. If the C-4 enrichment observed from the ^ spectrum 
(Fig. 4a) is ao, we have 



Mai/2) + 0 = ao 



(1) 



Quantitation of the glyoxylate and Krebs cycle activities 
can be approached as above using a mixture analysis, although 
the case is somewhat more complex. Table II summarizes the 
important labeling pathways and the labeling patterns which 
result. As indicated in the table, several compromises have 
been made in the analysis since in some cases different labeling 
paths can result in identical labeling patterns and thus cannot 
be distinguished. In addition, a small quantity of the labeling 
patterns indicated under w and y can also be produced after 

Table II 

Possible labeling patterns for C-2, C-3, and C-4 of [ X3 C] glutamate 
derived from [1 13 C] glucose 



Pathway 
designa- 
tion 


Labeling 
pattern* 


Source of label 


Pathway 


X 


o * * 
C1-C2-C3-C4 


Phosphoenolpyruvate carbox- 
ylase + V& Krebs cycle 


% 
58 


y 


* * 

* • 

C1-C2-C3-C4 


Fumarase activity plus single 
turn of glyoxylate cycle 


26 


z 


* • * * 

Cl-C2"C3-C< 


Multiple turns of Krebs or gly- 
oxylate cycles 


10 


w 


* * * 

Ct-C2-C3-C4 


Single turn 






* * 

C|-C2-C3-C< 


Krebs cycle 


6 



"Asterisk (*) indicates a labeling probability a derived from the 
enriched pyruvate pool. Zero (°) indicates a labeling probability fi 
derived from the enrichment of the bicarbonate pool. 

Table III 

Glutamate labeling and multiplet prob abilities 

Analytical expressions 



Enrichments 

C-l 

C-2 

C-3 

C-4 
Multiplets 



x0 + (z + w/2)-q 
(x + y/2 + w/2)*a 
(y/2 + z + w/2)-a 
(x + y + z + w)-a = a 



za(l-a) 2 + - 



C-l (doublet/singlet) (xjSa + zc?)/(x0(l - a) + zaU - a) + ~ 
C-2(doublet/singlet)" [xa0 + 2za J (l - a)}/(xa(l - fi) + | a + 

;■) 

C-2(quartet/singlet)* za 3 /[\xa(l - £> + |« + za(l - a) 2 + j a j 
C-3(doubIet/singlet> |j a 2 + 2za 2 (l - a) + ^ « 2 ]/[f «U - «) + 

za<l-a) 2 + |a<l-a)J 
C-3(triplet/singlet) z« 3 /^ a(l - a) + za(l - a) 2 + ^ a<l - a) j 
C-4 (doublet/singlet) ^| + z + ^aV^xa + ya^l - ^ + 

H)] 



za(l — a) + wa 



For the present study, a, = 90%, ao = 39%, giving f = 87%, or 
13%, of the glucose is metabolized via the hexose monophos- 
phate shunt. 



0 Represents total doublet intensity, although two separate dou- 
blets are resolvable due to different J a and t/23 values. 
* Quartet since c/12 7* t/23. 
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additional turns of the Krebs and glyoxylate cycles. Neverthe- 
less, the approach summarized represents a reasonable ap- 
proximation for quantitative pathway analysis. Since spectral 
inspection {Fig. 46) indicates that the number of multiply 
enriched molecules (z) is small, these approximations are 
probably fairly good. As can be seen from the table, negative 
correlations of the labeling of C-2 and C-3 can result from the 
Krebs or glyoxylate cycle, or from back reactions which reach 
fumarase. Krebs cycle activity will also result in a negative 
correlation between C-l and C-2 labeling. Analytical expres- 
sions for carbon labeling and relative multiplet probabilities 
can be derived in a straightforward manner, and the results 
are summarized in Table III. Using these expressions and the 
value of ao = 0.39 and the relative multiplet intensities in the 
13 C spectrum (Fig. 4b), values of 58, 26, 10 and 6% are obtained 
for x, y, z, and w, respectively. The relative contribution of 
these pathways should result in 13 C populations at C-2 and C- 
3 of 33 and 10%, respectively, in good agreement with those 
obtained directly from the *H NMR spectrum (Fig. 4a). 

We were unable to grow M. ammoniaphilum on acetate as 
the sole carbon source; however, cells grown on unlabeled 
glucose were washed and incubated with [l- l3 C]acetate. A 13 C 
NMR spectrum of the extracted glutamate is shown in Fig. 6. 
The primary features of this spectrum are the labeling ratio 
of 2:1 for C-5:C-1 and the absence of label for C-2, C-3, and C- 
4. These observations parallel the "C observations (24) on the 
glutamate labeling patterns obtained with another glutamate 
producer, Brevibacterium flavum. 

DISCUSSION 

Activity of the Hexose Monophosphate Shunt—The deter- 
mination that 13% of the glucose metabolized by M. ammon- 
iaphilum occurs via the hexose monophosphate shunt agrees 
well with the 10% value determined by analysis of the pyru- 
vate carboxyl group in metabolically blocked cells (18, 19). 
The result is significantly lower than the 25% figure deter- 
mined by recovering U C0 2 from [1- I4 C]- and [6- 14 C]glucose- 
grown cells. There is a variety of approximations inherent in 
both approaches (25). In particular, the hexose monophos- 
phate pathway may be somewhat underestimated since some 
of the pentose produced via this path will be used for nucleo- 
tide synthesis, rather than for production of pyruvate. 

Activities of the Krebs and Glyoxylate Cycles in the Me- 
tabolism of [l- l3 CJGlucose— The major metabolic routes 
found for the degradation of glucose involve glycolysis to 
pyruvate and phosphoenolpyruvate, formation of acetyl-CoA 
and oxaloacetate via decarboxylation and carboxylation, re- 
spectively, and the first third of the Krebs cycle to yield a- 
ketoglutarate and, subsequently, glutamate. Significant activ- 
ity of the glyoxylate shunt is required to explain the C-3 
labeling and reduced C-2 labeling of glutamate. The observed 
Krebs cycle activity was initially unexpected since this cycle 
is presumably shut down due to the reduction of a-ketoglu- 
tarate dehydrogenase activity. Evidence for the activity of 
this cycle comes primarily from the observed negative corre- 
lation in the labeling of C-l and C-2 of glutamate. 

The time dependence of the I3 C NMR spectra (Fig. 2) and 
the time-dependent increase in the C-3/C-4 intensity ratio 
(Fig. 5) is particularly interesting. Since the C-4 enrichment 
level is constant (ao = 0.39), this time-dependent ratio reflects 
a time-dependent increase in the relative contributions of the 
y, z, and/or w pathways. Although, in principle, it is possible 
to quantitate these as a function of time, there are several 
problems with this approach in practice. 1) the spectra cor- 
respond to cumulative pools of glutamate and therefore reflect 
cumulative average pathways; 2) 'H NMR spectra of purified 
glutamate in parallel with the 13 C data were not obtained; 3) 



the signal/noise in the 2-h spectra was insufficient to permit 
analysis of the C-l multiplet, which is of value for the dem- 
onstration of Krebs cycle activity. In view of these problems, 
a simplified approach was taken in which the glyoxylate and 
Krebs cycle activities were both represented by y since it 
becomes impossible to distinguish between y and w in the 
absence of C-l data. Based on this analysis, the cumulative 
contribution from the y pathway increases from -0 to a steady 
state value of -35%, with z increasing to -12% and x decreas- 
ing to a value of -53%. This increased cycling may reflect the 
decreasing glucose concentration which might release the 
inhibition of enzymes in both the glyoxylate and Krebs cycles, 
i.e. the time-dependent changes may reflect catabolite repres- 
sion. A similar time-dependent result was recently reported in 
a study of the glutamate derived from rat hearts perfused 
with [2- I3 C]acetate (26). 

Another interesting result is that the 13 C spectra show that 
the labeling at C-4 also appears to show a negative correlation 
with C-5. On the surface, there is no immediately obvious 
route for C-5 to be labeled, and one would expect C-5 to be at 
natural abundance regardless of the labeling at C-4. Instead, 
the analysis of the doublet center band intensities at C-5 leads' 
to an enrichment of 32 atom % 13 C at C-4, compared to the 
higher value of 39 atom % obtained by the J H satellites. This 
apparent discrepancy can be explained using the following 
pathway in which oxaloacetate breaks down to acetyl-CoA, 
which leads to an anticorrelated (vis a vis C-4) C-5 enrichment 
of glutamate (Structure 1). 
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The amount by which the acetyl-CoA pool is diluted with 
material derived from the oxaloacetate pool can be calculated 
from Equation 2. 

D/S = + { 'C n2a2 + fay (2) 

(1 " n [ {^T/" + 2a(1 " a) ] + f ? (1 ~ a) 

where D/S is the doublet to singlet ratio for 13 C-5, f is the 
fraction of acetyl-CoA derived directly from glucose, 1 - fthe 
fraction of acetyl-CoA derived from oxaloacetate, «, y, z, and 
w are as defined in Table II, and y is the natural abundance 
13 C level for C-5. Solving this equation for / from the data in 
Table II gives a value of -96%, suggesting that only -4% of 
the acetyl-CoA pool is derived from oxaloacetate. Finally, 
given f one can calculate the enrichment at C-5 from Equation 
3. 

% n C-5=fy+a(l-f)(y/2 + w/2 + z) (3) 

This value cannot be directly measured from the data, but a 
value of 1.8% 13 C can be obtained from Equation 3. Using a 
different approach, one can also calculate the enrichment at 
C-5 from the apparent enrichment at C-4 using the C-5 
multiplets and the actual enrichment at C-4 using the proton 
data. The discrepancy between these two values and the 
known distribution of label between the C-5 doublet and 
singlet from the amount of negatively correlated labeling at 
C-2 and C-3 of oxaloacetate, allow one to calculate an enrich- 
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ment of 1.7% at C-5, in excellent agreement with the value 
determined from Equation 3. 

Metabolism of [l- l3 CJ Acetate— In principle, it is possible to 
explain the 1:2 ratio for the C-l/C-5 labeling of glutamate 
derived from [l- l3 C]acetate on the basis of a contribution of 
the unlabeled metabolites initially present in the cell, which 
would dilute the C-l label to a greater extent than C-5. 
However, this explanation proves inadequate since the con- 
centration of initial unlabeled metabolites is far less than 
would be necessary to produce the observed difference. Thus, 
the labeling difference must reflect a steady state phenomenon 
in which the C-l label is systematically diluted. This conclu- 
sion can be made more convincingly in the case of B. flauum, 
in which acetate can be used as the sole carbon source, 
eliminating the possible complication of an initial pool of 
unlabeled metabolites. 

A steady state analysis of the acetate labeling experiments 
indicates that the data can be explained on the basis of 
significant Krebs cycle activity. This reflects the fact that only 
half of the C-l carbons processed through the w pathway 
become labeled (Table II). The possibility of high Krebs cycle 
activity as an explanation for similar carbon- 14 acetate growth 
results was considered by Shiio and Tsunoda (24), but subse- 
quently rejected on the basis of the demonstrated inhibition 
of this cycle under glucose growth conditions. The alternative ' 
possibility of two acetate pools was proposed. In view of our 
demonstration of the existence of Krebs cycle activity in the 
glucose growth experiment, the possibility of increased activ- 
ity on acetate medium seems reasonable. In contrast, the 
possibility of two acetate pools for a prokaryotic cell seems 
less likely. We thus conclude that, as in the apparent time- 
dependent nature of the labeling, differences in carbon source 
may also have a significant effect on the relative contributions 
of the different pathways, possibly reflecting catabolite repres- 
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To determine the in vivo fluxes of the central metabolism 
we have developed a comprehensive approach exclu- 
sively based on the fundamental enzyme reactions 
known to be present, the fate of the carbon atoms of 
individual reactions, and the metabolite balance of the 
culture. No information on the energy balance is re- 
quired, nor information on enzyme activities, or the di- 
rectionalities of reactions. Our approach combines the 
power of 1 H-detected 13 C nuclear magnetic resonance 
spectroscopy to follow individual carbons with the sim- 
plicity of establishing carbon balances of bacterial cul- 
tures. We grew a lysine-producing strain of Corynebac- 
terium glutamicum to the metabolic and isotopic steady 
state with [1- 13 C]glucose and determined the fractional 
enrichments in 27 carbon atoms of 11 amino acids iso- 
lated from the cell. Since precursor metabolites of the 
central metabolism are incorporated in an exactly de- 
fined manner in the carbon skeleton of amino acids, the 
fractional enrichments in carbons of precursor metabo- 
lites (oxaloacetate, glyceraldehyde 3-phosphate, eryth- 
rose 4-phosphate, etc.) became directly accessible. A 
concise and generally applicable mathematical model 
was established using matrix calculus to express all me- 
tabolite mass and carbon labeling balances. An appro- 
priate all-purpose software for the iterative solution of 
the equations is supplied. Applying this comprehensive 
methodology to C. glutamicum, all major fluxes within 
the central metabolism were determined. The result is 
that the flux through the pentose phosphate pathway is 
66.4% (relative to the glucose input flux of 1.49 mmol/g 
dry weight h), that of entry into the tricarboxylic acid 
cycle 62.2%, and the contribution of the succinylase path- 
way of lysine synthesis 13.7%. Due to the large amount 
and high quality of measured data in vivo exchange re- 
actions could also be quantitated with particularly high 
exchange rates within the pentose phosphate pathway 
for the ribose 5-phosphate transketolase reaction. More- 
over, the total net flux of the anaplerotic reactions was 
quantitated as 38.0%. Most importantly, we found that in 
vivo one component within these anaplerotic reactions 
is a back flux from the carbon 4 units of the tricar- 
boxylic acid cycle to the carbon 3 units of glycolysis of 
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30.6%. © 1996 John Wiley & Sons, Inc. 
Key words: intracellular fluxes • metabolite balance • car- 
bon labeling balance • lysine • anaplerotic reactions • 
NMR spectroscopy 

INTRODUCTION 

Amino acid production has a long tradition in biotechnol- 
ogy. In particular, the bacterium Corynebacterium glu- 
tamicum as well as mutants derived from this species are of 
industrial interest for the production of L-lysine and L-glu- 
tamate. 24,29 Currently, enzymatic and genetic analysis of 
this species is mostly directed to a single reaction 6, 15,34,45 or 
to a very limited set of reactions. 7,39 However, the effective 
implementation of any targeted genetic engineering of such 
overproducers requires a basic understanding of the central 
metabolism as a whole. Therefore, it is imperative to quan- 
tify the actual fluxes in vivo through all the reactions of the 
central metabolism. Unfortunately, experimental access to 
these reactions in particular (e.g., citrate synthase) is diffi- 
cult to achieve in vivo. 

One useful approach which has very recently attracted 
great attention is metabolite balancing. 8 ' 32,53 Flux estima- 
tions are obtained using measurements of substrate con- 
sumption and product formation, while taking into account 
the known demands of precursor metabolites for cellular 
growth (oxaloacetate, pyruvate, a-ketoglutarate, etc.). This 
approach can be well suited for particular purposes, i.e., 
when only specific pathways need to be considered 53,60 or 
if the contribution of flux for cellular growth is weak, as 
with hybridoma cells 49,65 and entire organs. 21 However, in 
the more complex bacterial systems, where product forma- 
tion often occurs simultaneously with growth and the con- 
sequent use of the entire set of metabolic pathways is re- 
quired, the large number of branching points in the meta- 
bolic network usually prevents a detailed analysis 43 In such 
cases, the application of metabolite balancing- either re- 
quires sets of reactions to be lumped together 59 or forces the 
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use of data from in vitro enzyme assays to answer such an 
important question as, for instance, that on the in vivo use 
of the glyoxylate cycle. Moreover, additional assumptions 
on the energy balance of the cell have to be included, thus 
leading to further uncertainties. 23,31 

Therefore, additional experimental data are required to 
establish a reliable stationary flux analysis. Isotopic tracer 
experiments are a prominent candidate for such a purpose. 
This technique is well established for differentiating be- 
tween the relative use of two pathways, such as between the 
pentose phosphate pathway and glycolysis. 18,63 The use of 
I3 C-labeled substrate combined with nuclear magnetic res- 
onance (NMR) spectroscopy is particularly instrumental, 
since the isotopic labeling of individual carbons is directly 
accessible. As an example, the measured incorporation of 
label from [l- 13 C]glucose into L-glutamate derived from the 
precursor metabolite a-ketoglutarate was proof of an in vivo 
activity of the tricarboxylic acid (TCA) cycle in a Coryne- 
bacterium strain, 60 where previous in vitro studies had sug- 
gested that it was inactive. 50 In a similar approach, the 
quantification of glyoxylate cycle and TCA cycle usage was 
possible. 61 We used 13 C NMR to determine the fractional 
enrichments in carbons of L-alanine (derived from pyruvate) 
and L-lysine (derived from pyruvate and oxaloacetate) pro- 
duced by C. glutamicum. 52 This enabled us to quantitate the 
relative use of the two competing pathways of L-lysine syn- 
thesis present in this bacterium. The data indicated that 
about three-quarters of the total flux of L-lysine occurs in 
vivo through the succinylase variant, although the activities 
of the enzymes involved are about 10 times less than that of 
the competing dehydrogenase variant. This once more il- 
lustrates the power of the NMR technique and presents a 
case where including information from in vitro analyses to 
decide on in vivo fluxes clearly would be invalid. 

The aim of the present study was to provide basic quan- 
titative information on the in vivo fluxes of as many reac- 
tions in the central metabolism of C. glutamicum as possi- 
ble. We chose an L-lysine producer of C. glutamicum for 
flux analysis, which has high productivity in batch cul- 
ture, 45 is very well characterized with respect to its genetic 
situation, 7,33 and for which data on its flux characteristics 
already exist. 52 In addition, information is also available on 
other L-lysine-secreting strains of C. glutamicum. 151 *' 29 - 53 
Here we present and apply an experimental methodology 
to make maximum and simultaneous use of the balancing 
and NMR technique, which allows the fluxes to be deter- 
mined in an overall network of reactions. A generally ap- 
plicable, completely automated methodology for data han- 
dling and analysis was developed. In contrast to other stud- 
ies, we grew C. glutamicum to metabolic steady state as 
well as to isotopic steady state to increase the reliability of 
the metabolite balance and the isotopic balances as much as 
possible, thereby avoiding the inherent uncertainties of a 
dynamic study. Furthermore, we made full use of the fact 
that the carbon skeletons of the central metabolites are fixed 
at clearly predetermined positions in amino acids. Therefore 



the particularly valuable 13 C labeling of central metabolites 
is obtainable, to which access is otherwise difficult. 

To obtain very precise and accurate labeling information, 
we determined 13 C enrichments by ! H NMR spectroscopy. 
Due to the high information content of the labeling data, 
neither assumptions on the stoichiometry of the energy me- 
tabolism (e.g., K ATP , P/O ration H + /e ratio), which are a 
matter of controversial debate, 3 1,42 nor assumptions based on 
any in vitro enzyme activity have to be taken into account. 
The metabolic model uses only the well-established chem- 
istry of the carbon transitions in the enzyme reactions of the 
central metabolism. Our description takes full account of 
label scrambling reactions where symmetrical molecules 
occur and also of back fluxes. The formalization developed 
by us is of general validity and ready to be applied to other 
systems. It enables a detailed picture of the in vivo fluxes to 
be obtained under L-lysine-producing conditions in C. glu- 
tamicum. 

MATERIALS AND METHODS 
Organism and Cultivation 

The organism used was the L-lysine-producing Corynebac- 
terium glutamicum strain MH20-22B. 44,45 To stabilize the 
Leu - character of the strain, the plasmid pKK74I2 33 was 
introduced via conjugation into MH20-22B. Twelve kana- 
mycin-resistant clones were obtained due to vector se- 
quences integrated in the chromosome since pKK74I2 is 
nonreplicative in C. glutamicum. Because this plasmid car- 
ried an internal leuA fragment of 647 bp, the vector was 
integrated into the chromosomal leuA locus. The resulting 
strains were confirmed to be inactive in the leuA encoding 
isopropylmalate synthase, thus securing the stability of the 
Leu~ character and ensuring that no leucine biosynthetic 
intermediate accumulates. The resulting strain C. glu- 
tamicum MH20-22B leuA was pregrown in a medium con- 
taining 3.7% brain-heart infusion (Difco, Detroit, MI), 2% 
glucose, and 20 mg/L kanamycin. 

Continuous cultivation was started by inoculating 120 
mL of a late exponential growth phase culture into the cul- 
tivation vessel containing 205 mL medium with (amounts 
per liter): 2 g urea, 150 mg leucine, 20 mg kanamycin 
monosulfate, 30 mg protocatechuic acid, 0.2 mg biotin, 
75 mg Titriplex II, 0.02% polypropylene glycol, 10 g 
(NH 4 ) 2 S0 4 , 1.5 g NaH 2 P0 4 ■ H 2 0, 1 g KCI, 250 mg 
MgS0 4 . H O, 1 mg ZnS0 4 * 7H 2 0, 0.5 mg CuS0 4 , 0.2 mg 
NiCl 2 ■ 6H 2 0, and 0.12% HC1. The medium was supple- 
mented with 7.55 g/L of either unlabeled glucose (Sigma, 
St. Louis, MO) or [l- 13 C]glucose (99% enriched, Cam- 
bridge Isotope Laboratories, Cambridge, MA). 

Continuous culture with a dilution rate of 0. 10 h~ 1 and 
325 mL culture volume was carried out in a BioFlo C30 
reactor (New Brunswick Scientific Co., Edison, NJ) in a 
20% 0 2 and 80% N 2 synthetic air atmosphere at 30°C. The 
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rotation speed was 400 rpm and the aeration rate was 41.0 
L/h. The pH of the culture was maintained at 7.0 with 2 N 
NaOH by an automatic pH control system (Bioengineering 
AG, Wald, Switzerland). Peristaltic pumps were used for 
media supply and withdrawal. On-line observation of the 
C0 2 concentration in the effluent gas indicated metabolic 
steady-state conditions after 77 h of cultivation so that the 
medium containing unlabeled glucose was changed to a 
medium supplemented with [l- 13 C]glucose while maintain- 
ing the other components unchanged. After 30 h of subse- 
quent cultivation, the cells were harvested by centrifuga- 
tion, washed twice, lyophilized, and stored at -20°C. Un- 
less otherwise stated, all chemicals were obtained from 
Merck (Darmstadt, Germany) or Boehringer Mannheim 
(Germany). 

Analytical Methods 

Dry weights were obtained after centrifuging 5 mL of 8-mL 
samples (13000g, 10 min), washing the pellets twice with 
50 mM KH 2 P0 4 /K 2 HP0 4 buffer, pH 7.0, and drying for 48 
h at 100°C. The supernatants were analyzed for glucose and 
amino acids. Amino acids were analyzed as ortho- 
phthaldialdehyde derivatives by reversed -phase liquid chro- 
matography. 45 Glucose was assayed enzymatically. 3 The 
carbon dioxide content of the effluent air was measured 
with an infrared C0 2 analyzer (Hartmann und Braun, 
Frankfurt, Germany). The total carbon content of the 
freeze-dried biomass was determined using a CHNS 932 
elemental analyzer (Leco Co., St. Joseph, MI). The 13 C0 2 / 
12 C0 2 ratio in the purified gas phase 12 of the culture was 
estimated by analyzing the head space with a Sira II mass 
spectrometer (Fisions VG Isotech Instruments, Mid- 
dlewich, United Kingdom). 

Preparation of Amino Acids for 
NMR Spectroscopy 

Lyophilized cells prepared from 55-mL culture volume 
were resuspended in 160 mL of 6 jV HC1 and hydrolyzed at 
110°C in a sealed glass for 12 h. After the addition of 480 
mL of demineralized water the hydrolysate was evaporated 
to dryness and redissolved in 8 mL of 0.2 M triethylamine 
(pH 2.5). Amino acids were fractionated 52 by cation ex- 
change chromatography on an Ultrapac ll-p,m resin col- 
umn (Pharmacia Biotech GmbH, Freiburg, Germany) with 
triethylamine buffer (0.2 M, pH 3.2-10.5) used for elution. 
Freeze-dried powders of the amino acids were dissolved in 
700 p.L D 2 0 and passed through a 0.2-|xm Dyna Gard filter 
(Microgon Inc., CA) to remove insoluble impurities. 

NMR Spectroscopy 

All NMR experiments were performed using a Bruker 
AMX 400-WB spectrometer equipped with a 5-mm inverse 
probe head for 'H. Carbon- 13 enrichments were determined 



from normal and 13 C-decoupled 400-MHz proton NMR 
spectra by relating the area of the I3 C satellite signals in the 
difference spectrum to the area of the corresponding proton 
signals in the decoupled spectrum as described. 52 Broad- 
band l3 C decoupling was accomplished using the GARP-1 
composite pulse decoupling scheme. 48 

RESULTS 

Biochemical Flux Model 

The fluxes connecting the various precursor metabolites 
within the central metabolism are given in Figure 1. These 
basic pathways have been identified by a variety of bio- 
chemical analyses carried out with C. glutamicum. 24 This 
network of reactions provides the biochemical model for the 
flux analysis. The established pathways are that of glycol- 
ysis, the pentose phosphate pathway (PPP), TCA cycle, and 
glyoxylate cycle. The anaplerotic reactions are currently a 
matter of much debate. Phosphoenolpyruvate carboxylase, 
of which the existence had been shown by enzyme studies, 9 
was concluded to be of utmost importance in C. glutamicum 
from metabolite balance studies 53 but can be deleted with- 
out any consequences for product formation, 14-34 although 
an alternative reaction replenishing the TCA cycle has not 
yet been identified. 20 According to the various enzyme de- 
terminations, the principal catalytic activities carrying 
fluxes connecting phosphoenolpyruvate with oxaloacetate, 
pyruvate with oxaloacetate, and pyruvate with malate are 
present in C. glutamicum. 

For our purpose of in vivo flux analysis, not relying on in 
vitro activities, we regarded phosphoenolpyruvate plus py- 
ruvate as one pool, which seems justified because the de- 
termined 13 C enrichments were equal within experimental 
error (Table I). Oxaloacetate and malate were also treated as 
a single pool. The notation of all fluxes is given in Table II. 
The metabolite branch points included in the model are 
identified by rectangles in Figure I . Because, in our model, 
acetyl coenzyme A is always directly and exclusively de- 
rived from pyruvate while no reverse flux of acetyl coen- 
zyme A to pyruvate can occur, its 13 C labeling must exactly 
reflect that in pyruvate. Therefore, acetyl coenzyme A was 
not considered to be a metabolic branch point and was not 
incorporated as a separate metabolite pool in our model. 
Likewise, all other metabolites not representing branch 
points (e.g., 6-phosphogluconate, succinate, etc.) were not 
included as separate pools in the model. In no way does this 
approach affect the 13 C distribution in our network. Of 
course, for, e.g., experiments with acetate supply in the 
medium, one definitely would require a separate acetyl co- 
enzyme A pool to be included in the model. 

Within the biosynthetic reactions removing metabolites 
for synthesis of the macromolecular components, CO z gen- 
erating and consuming reactions occur. Therefore, in total 
11 anabolic reactions were incorporated into the model 
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Figure 1. Metabolic model and metabolite balance used for in vivo flux determination in the 
L-lysine-producing Corynebacterium glutamicum strain MH20-22B leuA. Numbers in oval black 
symbols state measured fluxes of substrate, products, and biomass precursor metabolites. The latter 
were inferred from the measured biomass concentration using the stoichiometric approach described 
in the text. All fluxes are expressed as percentage of the glucose uptake rate (1 .49 mmol/g dry weight 
h). Flux identifiers indicated (v,, v 2 , . . .) refer to the enzyme reactions listed in Table I. Abbrevi- 
ations: glc, glucose; lys, lysine; co2, C0 2 ; g6p, glucose 6-phosphate; fl6bp, fructose 1,6- 
bisphosphate; f6p, fructose 6-phosphate; gap, glyceraldehyde 3-phosphate; pyr, pyruvate; oaa, ox- 
aloacetate; fum, fumarate; icit, isocitrate; akg, a-ketoglutarate; ru5p, ribulose 5-phosphate; ri5p, 
ribose 5-phosphate, xu5p, xylulose 5-phosphate; s7p, sedoheptulose 7-phosphate; e4p, erythrose 
4-phosphate; ddh, diaminopimelate dehydrogenase pathway; sue, succinylase pathway. 



(Fig. 1). The anabolic contribution was determined from the 
observed biomass yield on glucose and stoichiometric data 
(see below). It is important to note that all reaction steps 
shown, as well as all the concomitant carbon dioxide fluxes, 



were included in the metabolite balances as well as in the 
13 C label distribution calculations. Thus, for example, the 
C0 2 balance included carbon dioxide released and/or fixed 
in the oxidative PPP, in the TCA cycle, in anaplerosis, as 
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Table I. Carbon- 1 3 label data for central metabolic intermediates derived from the l3 C enrichments 
of the purified amino acids by 'H NMR (Table IV). 



l3 C enrichment (%) by carbon position 



Metabolite 


C-l 


C-2 


C-3 


C-4 


C-5 C-6 


e4p 


2.0 ± 1.0 


3.6 ± 1.0 


2.0 ± 1.0 


16.7 ± 2 




gap 


2.9 ± 0.2 


2.6 ± 0.1 


26.7 ± 0.2 






pep 


n.d. 


2.5 ± 0.5 


26.5 ± 1.0 






pyr 


n.d. 


3.0 ± 1.0 


26.4 ± 0.5 






akg 


n.d. 


24.1 ± 0.3 


11.1 ± 0.5 


28.1 ± 0.6 


n.d. 


oaa 


n.d. 


7.7 ± 2.0 


20.9 ± 2.0 


16.8 ± 2.7 




lys 


n.d. 


6.8 ± 0.2 


21.9 ± 0.3 


18.9 ± 1.0 


22.2 ± 1.0 5.6 ± 0.3 


co 2 


23.0 ± 0.4 











Note: ,3 C fractional enrichment of C0 2 was determined by mass spectroscopy. Abbreviations (see 
Fig. 1): pep, phosphoenolpyruvate; n.d., not determined. 



well as in the anabolic reactions. The resulting net C0 2 
production of the culture was measured. 

Depending on thermodynamic considerations as found in 
standard biochemistry textbooks, in our model a number of 
reactions were formulated as being reversible (see below) 
while others were considered irreversible. This was done 
because the presence of back reactions may severely influ- 
ence 13 C labeling. 



Strain and Quantification of Extracellular Fluxes 

Corynebacterium glutamicum strain MH20-22B is an ex- 
cellent L-lysine producer, 45 which is characterized by a lysC 



mutation and thereby feedback-resistant aspartate kinase. 
This and an altered export carrier 6,45 result in a high flux of 
aspartate to L-lysine with this strain. In addition this strain 
is auxotrophic for L-leucine, 7 ' 33 but this character has been 
observed to be susceptible to reversion. Therefore, we in- 
troduced the defined leuA mutation into strain MH20-22B 
- by gene-directed mutagenesis as described in the Materials 
and Methods section. The resulting strain was cultivated in 
a chemostat at a dilution rate of 0.1 h _1 and found to be 
stable for at least 107 h. L-Lysine production was not re- 
duced with cultivation time, as is often observed with over- 
producers. 25 L-Lysine determinations were made, as well as 
dry weight determinations, and the C0 2 release was fol- 
lowed on-line. 



Table II. Enzyme reactions corresponding to the flux identifiers in the biochemical model of 



Figure 1. 


v i 


glucose 6- phosphate isomerase 


v 2 


phosphofructokinase 


v 3 


fructose 1 ,6-bisphosphate aldolase, triose phosphate isomerase 


v 4 


glyceraldehyde 3-phosphate dehydrogenase, bisphosphoglycerate kinase, phosphoglycerate 




mutase, enolase, pyruvate kinase, phosphoenolpyruvate:sugar phosphotransferase system 


v 5 


glucose 6-phosphate dehydrogenase, lactonase, 6-phosphogluconate dehydrogenase 


v 6 


ribose 5-phosphate isomerase 


v i 


ribulose 5-phosphate epimerase 


v 8 


transketolase (xu5p + e4P «-» gap + f6p) 


v 9 


transaldolase (gap + s7p e4p + f6p) 


v 10 


transketolase (xu5p + ri5p *■* s7p + gap) 


V l] 


pyruvate dehydrogenase complex, citrate synthase, aconitase 




isocitrate dehydrogenase 


v 13 


a-ketoglutarate dehydrogenase, succinate thiokinase, succinate dehydrogenase 




fumarase, malate dehydrogenase 




isocitrate lyase, malate synthase, fumarase, malate dehydrogenase 




phosphoenolpyruvate carboxylase, phosphoenolpyruvate carboxykinase, malic enzyme. 




oxaloacetate decarboxylase 


ddh 


transaminase A, aspartate kinase, aspartate semialdehyde dehydrogenase, 



dihydrodipicolinate synthase, dihydrodipicolinate reductase, diaminopimelate 
dehydrogenase, diaminopimelate decarboxylase (lysine production over the 
diaminopimelate dehydrogenase pathway) 
sue same as ddh, except diaminopimelate dehydrogenase replaced by the sequence of enzyme 
reactions tetrahydrodipicolinate succinylase, succinyl aminoketopimelate transaminase, 
succinyl diaminopimelate desuccinylase, diaminopimelate epimerase (lysine production 
via the succinylase pathway) 
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Because all glucose was consumed in the culture, the 
uptake rate for glucose could be directly calculated as 1 .49 
± 0.03 mmol/g dry weight h from the dry weight content of 
the broth (2.8 g dry weight/L), glucose concentration in the 
medium (7.55 g glucose/L), and the dilution rate (0. 1 h~ '). 
For simplicity, this value was normalized to 100% and all 
fluxes reported in this study are molar fluxes expressed as 
percentage values of this glucose uptake rate. The L-lysine 
concentration accumulating in the culture was 7.66 mAf. 
This corresponds to an L-lysine excretion flux of 18.3 ± 
0.4%. The C0 2 efflux was determined from the measured 
C0 2 release (4. 10 mmol C0 2 /g dry weight h) to be 275. 1 ± 
13.8%. From the biomass concentration in the fermentor 
(2.80 g dry weight/L) and thexarbon content of the biomass 
(0.408 g carbon/g dry weight) it was concluded that 35.3% 
of the carbon originating from glucose was used for ana- 
bolic purposes. Furthermore, it was calculated that of the 
total supplied carbon 45.8% was converted to C0 2 and 
18.3% was used for the synthesis of L-lysine. Therefore, 
after summation, the carbon balance was closed up to 



99.4%. Accordingly, only traces of other amino acids, keto 
acids, or organic acids were detected. Trehalose, often 
found with L-lysine producers, 59 could not be detected ei- 
ther. 

Fluxes for Precursor Use 

As determined, in the culture about one-third of the carbon 
is used for cellular growth. Therefore, the requirements of 
pyruvate and oxaloacetate and that of the other precursor 
metabolites which are removed at various steps in the cen- 
tral metabolism have to be quantitated. For the Escherichia 
coli cell a detailed stoichiometric approach is described 
where the cell is decomposed into 12 precursor metabo- 
lites. 30 These data have been used in a metabolite balance 
study with C. glutamicum ATCC 21253. 59 We directly 
quantitated the amino acid composition of C. glutamicum 
MH20-22B (Table III) to obtain the amounts of the 12 pre- 
cursor metabolites required for amino acid synthesis. These 
requirements are based on the standard biosynthetic path- 



Table III. Precursor and carbon dioxide requirements (u.mol/g dry weight) for biomass synthesis of C. glutamicum. 



Stoichiometry 



Precursor 


Amount g6p 


f6p 


ri5p 


e4p 


gap 


Pga 


pep 


pyr 


accoa 


oaa 


akg 


co2 


nadph 


Ala 


606 














1 










1 


Arg 


189 




















i 


1 


4 


Asx 


399 


















1 






1 


Cys 


87 










1 














5 


Glx 


806 




















i 




1 


Gly 


361 










1 














1 


His 


71 




1 




















I 


He 


202 














1 




1 




-1 


5 


Leu 


0 














2 


1 






-2 


2 


Lys 


202 














1 




1 




-1 


4 


Met 


146 


















1 






8 


Phe 


133 






1 






2 










-1 


2 


Pro 


170 




















l 




3 


Ser 


225 










1 














1 


Thr 


275 


















1 






3 


Trp 


54 




1 


1 






1 










-1 


2 


Tyr 


81 






1 






2 










-1 


2 


Val 


284 














2 








- 1 


2 


Dmp 


146 














1 




1 






4 


Protein 


0 


0 


125 


268 


0 


673 


482 


1,724 


0 


1,370 


1,165 


-767 


9,718 


RNA 






630 






368 








262 




630 


427 


DNA 






100 






50 








50 




100 


200 


Lipids 










129 


129 






2,116 






-2,116 


3,612 


LPS 


51 


16 


24 






24 


24 




329 






-329 


470 


Peptidoglycan 




55 










28 


83 


55 


28 


28 


-55 


193 


Glycogen 


154 
























0 


Cl-units 












49 














49 


Polyamines 






















59 




180 


Precursor 


g6p 


f6p 


ri5p 


e4p 


gap 


Pga 


pep 


pyr 


accoa 


oaa 


akg 


co2 


nadph 


Total 


205 


71 


879 


268 


129 


1,293 


534 


1,807 


2,500 


1,710 


1,252 


-2,537 


14,849 



Note: The amino acid composition was determined from C. glutamicum strain MH20-22B leuA grown in continuous culure as described. The data 
account for the very high intracellular glutamate pool of 160 mAf present in this strain. The data on other cellular constituents were taken from ref. 30. 
The experimentally determined specific glucose consumption was 14,949 u.mol/g dry weight. Abbreviations (see Fig. 1): pep, phosphoenolpyruvate; accoa, 
acetyl coenzyme A; nadph, (J-nicotinamide adenine dinucleotide phosphate (reduced form). 
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ways of amino acids which are identical to those in E. coli 
(except for L-lysine synthesis 45 ). The total amino acid con- 
tent represents the major part (64%) of the carbon content of 
C. glutamicum. Together with data on the other cellular 
polymers 30 which are also shown in Table III, a good fit to 
the experimentally determined total carbon content of the 
cells was found, confirming once more the established total 
carbon balance. 

The demands of precursor metabolites as determined us- 
ing the stoichiometric data in Table HI in combination with 
the measured biomass yield on glucose for C. glutamicum 
MH20-22B, as well as the directly measured fluxes, are 
given in the black ovals of Figure 1 . Inspection of the pro- 
posed anabolic pathways 30 shows that with a number of 
reaction steps involved C0 2 is released or fixed. Therefore, 
careful account was taken of the fluxes of carbon dioxide 
associated with those reactions in order to keep accurate 
metabolite balances and take accurate account of the con- 
comitant transfer of 13 C label. The stoichiometric parame- 
ters for these fluxes are also included in Table III. Thus, 
where appropriate, anabolic precursor metabolite fluxes 
were split up into two components; one with involvement of 
carbon dioxide release (or fixation) and one without such 
involvement. This is indicated in Figure 1. 

As an example, the total flux of pyruvate withdrawn from 
the central metabolism for anabolic purposes, which was 
calculated to be 41.0%, is composed of a flux of 23.0% 
involving the release of C0 2 and of a flux of 18.0% where 
no carbon dioxide is involved. These values were deter- 
mined as follows: The total stoichiometric demand for glyc- 
eraldehyde 3-phosphate, phosphoenolpyruvate, pyruvate, 
and acetyl coenzyme A (Table III), all condensed into py- 
ruvate, is 6,132 |xmol C3 unit/g biomass, which amounts to 
41% of the measured specific glucose uptake of 14^949 
junol/g dry weight. Of this total of 6,132 fxmol/g dry 
weight, 3,456 jxmol C3 unit/g dry weight, or 23% normal- 
ized on specific glucose consumption, is involved in decar- 
boxylating biosynthetic reactions, i.e., with release of C0 2 . 
This value of 3,456 u.mol C3 unit/g dry weight is obtained 
by summing the contributions given in Table III for the 
synthesis of L-isoleucine, L-lysine, L-phenylalanine, L-tryp- 
tophane, L-tyrosine, L-valine, lipids, lipopolysaccharids, 
and peptidoglycan. The remaining 18% represents reactions 
where the complete C3 unit is incorporated in cellular con- 
stituents. Likewise, the total flux of a-ketoglutarate di- 
verted from the central metabolism was 8.2%, of which a 
component of 1.2% involves C0 2 fixation for L-arginine 
biosynthesis (Table III). 

Quantification of Fractional Enrichments in 
Precursor Metabolites 

In the steady state of the culture identical to that from which 
the metabolite balance was made, the unlabeled glucose 
was replaced by [l 13 C]glucose and the cells were grown for 
a further 30 h, which is equivalent to three complete volume 
changes. Stability of the culture was ensured by monitoring 



the L-lysine concentration, cell density, and C0 2 . Applying 
first-order wash-out kinetics, the 95% isotopic equilibration 
calculated was accounted for in the data analysis. The cells 
were harvested and hydrolyzed, and the amino acids were 
isolated on a preparative scale by cation exchange chroma- 
tography. Since the amino acids are synthesized by prede- 
termined, exactly defined pathways from precursor metab- 
olites, as established for £. coli 30 and confirmed in all stud- 
ies with C. glutamicum 24 the fractional enrichments in 
individual carbons of the amino acids directly reflect the 
fractional enrichments in the precursors. (In fact, all deter- 
minations were in full accord with the standard biosynthetic 
pathways). For instance, the a-ketoglutarate skeleton is di- 
rectly accessible by analysis of L-glutamate (Table IV), or 
glyceraldehyde 3-phosphate by analysis of l- serine. 

Because in several cases one specific precursor metabo- 
lite is incorporated in several amino acids, redundant mea- 
surements were possible, further increasing the accuracy of 
the determinations. For example, pyruvate was accessible 
by analysis of L-alanine, L-valine, and L-isoleucine. We 
determined the fractional 13 C enrichments by ! H NMR 
spectroscopy from the protonated carbons of purified L-glu- 
tamate, L-aspartate, L-threonine, L-lysine, l- leucine, L-iso- 
leucine, L-alanine, L-valine, L-serine, L-glycine, and L-phe- 
nylalanine. Thirty-one out of the 52 carbons of these amino 
acids were measured. The results are shown in Table IV. 
From these data the 13 C enrichments of the precursor me- 
tabolites erythrose 4-phosphate, glyceraldehyde 3-phos- 
phate, pyruvate, oxaloacetate, and a-ketoglutarate were 
available. The condensed results are given in Table I. A 
typical recorded NMR spectrum is shown in Figure 2. Fur- 
ther important information was the fractional enrichment in 
the carboxyl group of glycine, which is identical to the C-l 
position of glyceraldehyde 3-phosphate. This was obtained 
with 'H NMR by quantitating the small fractional splitting 
of the H-2 signal unequivocally attributable to heteronuclear 
scalar coupling to 13 C-1 of glycine. As additional informa- 
tion, the 13 C fractional enrichment of C0 2 was also deter- 
mined. This was achieved by mass spectroscopy and the 
enrichment estimated as 23.0 ± 0.4%. 

Computer Modeling of the Metabolite Balance 
and 3 C Label Distribution 

Manual inspection of the data already provides useful in- 
formation. Thus, the label in C-2 of pyruvate cannot be 
explained by assuming the sole transfer of carbon step by 
step through glycolysis or the PPP. Instead a back flux 
within the PPP or a back flux from the TCA cycle to py- 
ruvate is necessarily required. Similarly, inspection of the 
label distribution in erythrose 4-phosphate conclusively 
shows that back fluxes occur in vivo within the PPP. How- 
ever, for accurate quantification by detailed simultaneous 
consideration of metabolite, as well as 13 C label balances 
throughout the complete metabolic network, a general 
mathematical formalism had to be developed. The mathe- 
matical principles and a simple application example are de- 
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Table IV. Carbon-13 label data of the purified amino acids as determined by l H NMR. 

,J C enrichment (%) by carbon position 



Amino acid 


C-l 


C-2 


C-3 


C-4 


C-5 




C-6 


Gly 


2.9 ± 0.2 


2.6 ± 0.1 












Ser 


n.d. 


1.5 ± 1.5 


26.7 ± 0.2 










Ala 


n.d. 


3.0 ± 1.0 


26.4 ± 0.5 








t 


Val 


n.d. 


n.d. 


n.d. 


26.3 ± 0.5 


26.3 ± 


0.5 




He 


n.d. 


7.9 ± 2.0 


n.d. 


n.d. 


16.5 ± 


0.8 


26.0 ± 0.8 


Leu' 


n.d. 


1.1 ± 0.3 


0.9 ± 0.7 


1.1 ± 0.6 


1.0 ± 


0.2 


1.0 ± 0.2 


Glu 


n.d. 


24.1 ± 0.3 


11. 1 ± 0.5 


28.1 ± 0.6 


n.d. 






Asp 


n.d. 


7.8 ± 2.0 


20.5 ± 2.0 


n.d. 








Thr 


n.d. 


7.4 ± 1.0 


21.3 ± 0.9 


16.8 ± 2.7 








Lys 


n.d. 


6.8 ± 0.2 


21.9 ± 0.3 


18.9 ± 1.0 


22.2 ± 


1.0 


5.6 ± 0.3 


Phe b - C 


n.d. 


2.5 ± 0.5 


26.5 ± 1.0 


n.d. 









Note: Measurement inaccuracy varied strongly due to differences in overlap between the spectra. 
Abbreviation: n.d., not determined. 

■Leucine with natural abundance, since it was supplied in the medium (see Methods section). 

The average enrichment of Phe C-5 and C-9 was 21.6 ± 1.0%, where one of these carbons is 
derived from phosphoenolpyruvate C-3 and therefore should be labeled according to Phe C-3. 

c Of the remaining three carbons Phe C-6, C-7, and C-8, which are derived from erythrose 
4-phosphate, 0.047 are 13 C. 



scribed in the Appendix. According to the biochemical net- 
work of C. glutamicum, the three fluxes corresponding to 
nonstationary metabolite pools (glucose, lysine, C0 2 ), as 
well as 11 depleting fluxes for biomass synthesis, were 



defined as well as 14 nonaccumulating metabolite pools 
(identified by squares in Fig. 1), which are connected by 16 
fluxes v ( (Table II). For instance, the PPP was not simply 
represented as a glycolysis bypass, but instead included the 






A. 




i 



ppm 



4.0 



— I — 
3.5 



— 1 — 
3.0 



— I — 

2.5 



2.0 



Figure 2. Illustrative 'H NMR spectrum of L-glutamate extracted from the isotopically enriched cell protein of C gluiamicum strain MH20-22B leuA. 
Lower trace, l3 C decoupled spectrum. Upper trace, difference between normal and 13 C decoupled spectra showing the 13 C satellites (2.6 times enlarged). 
Signal assignments: a, b, c: H-2, H-3, and H-4 of glutamate, respectively; d: CH 2 of residual triethylamine buffer used for elution. 



118 BIOTECHNOLOGY AND BIOENGINEERING, VOL 49, NO. 2, JANUARY 20, 1996 



metabolites glucose 6-phosphate, ribulose 5-phosphate, ri- 
bose 5-phosphate, xylulose 5-phosphate, sedoheptulose 
7-phosphate, erythrose 4-phosphate, glyceraldehyde 
3-phosphate and fructose 6-phosphate together with the re- 
actions of glucose 6-phosphate isomerase, ribulose 5-phos- 
phate epimerase, ribose 5-phosphate isomerase, transketo- 
lase, transaldolase, as well as the enzymes of the oxidative 
PPP branch. Except for the oxidative branch, all the PPP 
reactions were modeled as being reversible; i.e., they were 
allowed to carry a back flux simultaneously with a forward 
flux. This enabled the correct modeling of exchange fluxes 
as well as of their influence on 13 C label distribution in these 
reactions (see Appendix). In addition, the reactions v 3 , 
v 4 , and v i6 (Table II) were considered to be reversible, as 
was the reaction sequence of fumarase and malate dehydro- 
genase (v 14 ). The extent to which a reaction runs in a re- 
versible manner is represented by the exchange rate, which 
we define as the ratio of the flux component running simul- 
taneously backward and forward in a reaction to the net 
forward flux. 

Whereas in most studies only extreme situations have 
been analyzed so far, i.e., there is either no exchange or 
very large exchange rates prevail, 46 our analysis also mod- 
els intermediate exchanges, simply because the exact de- 
pendence on the exchange flux of the l3 C label of the in- 
termediates involved is adequately described by our math- 
ematical equation system. In addition, full label scrambling 
was assumed within the reactions converting ot-ketoglut- 
arate to fumarate and fumarate to oxaloacetate. In total, 82 
carbon atoms are included in the model and 26 correspond- 
ing measurement values were available. In addition, since 
we supplied the culture with C0 2 -free gas to make use of 
the respective fractional enrichment in evolved C0 2 , also 
the origin and fate of C0 2 in the biosynthetic pathways, was 
individually included in the model and formulated accord- 
ingly as described above. The assumptions used for mod- 
eling are given in Table V. 



Table V. Special features of the in vivo flux analysis for C. glutamicum 
and assumptions used for flux estimations. 

1. Only data for metabolite and carbon atom balances are used. Energy 

consideration are excluded. 

2. Specific enzyme activities and kinetic information on enzymes, as 

obtained from in vitro studies, are excluded. 

3. No directionality is preset for any of the fluxes. 

4. In specific cases fluxes are preset to be irreversible, but with 

arbitrary directionality, (These fluxes appear without hexagonal 
symbol in Fig. 3.) 

5. The conversion of citrate to a-ketoglutarate is accompanied by a 

conserved orientation of the citrate carbon skeleton (i.e., no 
scrambling of carbons). 

6. The orientation of carbons in reactions involving fumarate is 

randomized, i.e., the carbons C-l and C-4 as well as C-2 and C-3 
of fumarate are scrambled. 

7. The conversion of tetrahydrodipicolinate to diaminopimelate through 

the succinylase variant of L-lysine synthesis is accompanied by 
scrambling of the carbons C-l and C-7, C-2 and C-6, and C-3 
and C-5. 



Based on the experimental data and the metabolic net- 
work, the l3 C labeling in alt carbons of all branch point 
metabolites was related to the flux and exchange parameters 
of the metabolic model. This was coupled to an optimiza- 
tion program that, in addition to the net fluxes, also opti- 
mized the exchange rates such that altogether the predicted 
I3 C labelings fitted best to the experimentally observed val- 
ues. The optimal estimations for fluxes and exchange rates 
were obtained by minimizing the sum of squares of the 
deviations between model values and the measured data (for 
fluxes as well as t3 C labelings). These deviations were nor- 
malized to the experimental errors, i.e., a difference (for a 
given single carbon label or single flux value) of equal 
magnitude to the experimental error in the corresponding 
measurement value contributed one unit to the total residual 
sum of squares. 

Estimated Fluxes in Central Metabolism 

A series of flux estimations was made using different sets of 
starting values. Independent of starting values, reproducible 
results were obtained for the entire fluxes of the glycolysis, 
-PPP, TCA cycle, glyoxylate cycle, and split pathway of 
L-lysine synthesis. Therefore, the fluxes in these pathways 
of the central metabolism seemed to be well defined by the 
experimental data. For instance, a nonoxidative PPP flux 
for the supply of erythrose 4-phosphate and ribose 5-phos- 
phate is also possible, 54 ' 63 but evidently did not occur, al- 
though our model is well able to accommodate it. As can be 
seen from Table VI, the experimentally determined enrich- 
ments are in excellent agreement with the corresponding 
enrichments as estimated by the parameter fitting algorithm. 

In Figure 3 the values for the net fluxes are given in the 
black ovals. Thus, of the 100.0% glucose flux entering the 
cell, 32.2% is metabolized over glycolysis (vj) and 66.4% 
over the oxidative branch of the PPP (v 5 ). The remaining 
flux of 1.3% was used for cell material synthesis. The net 
flux from glyceraldehyde 3-phosphate to pyruvate is 
160.9%, which is higher than the glucose input flux, since 
of course 2 mol of C3 unit can basically be generated per 
mole of glucose. Thus, 39. 1% of flux, expressed in C3 unit 
flux, had been withdrawn from glycolysis and the PPP for 
biomass synthesis and C0 2 formation. The flux entering the 
TCA cycle (v n ) is 62.2%. The glyoxylate cycle flux is 
weak (0.2%) but is significantly and reproducibly greater 
than zero. This is proof that the current assumptions used 
for modeling on the in vivo use of the TCA cycle are rea- 
sonable. The flux of pyruvate and oxaloacetate toward ex- 
ported lysine was 4.9% through the dehydrogenase reaction 
and 13.7% through the succinylase reaction. 

The estimated values for the exchange rates are stated in 
hexagonal symbols in Figure 3. As stated before, exchange 
rates were normalized to the net flux, i.e., a value of 10 
means that the reaction back flux is 10 times larger than the 
net flux v f . It has been shown 46 that large exchange rates can 
only be estimated within large confidence bounds. How- 
ever, the existence and approximate magnitude can be well 
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Table VI. Estimated l3 C enrichments of central metabolic intermediates resulting from the best fit 
to the experimental data. 



I3 C enrichment (%) by carbon position 



Metabolite 


C-1 


C-2 


C-3 


C-4 


C-5 C-6 


e4p 


2.5 


2.0 


1.9 


15.3 




gap 


2.7 


2.6 


26.3 






pyr 


2.8 


2.7 


26.3 




f 


akg 


17.3 


22.6 


9.8 


26.3 


2.7 


oaa 


7.2 


9.8 


22.6 


17.3 




lys 


5.5 


7.1 


24.0 


17.3 


24.9 5.3 


co 2 


21.6 











Note: Abbreviations: see Table I. 



assessed. Such very high exchange rates of 8.8 in the glu- 
cose 6-phosphate isomerase reaction, 387.1 in the glycer- 
aldehyde 3-phosphate to pyruvate conversions, and 46.8 in 
the ribose 5-phosphate transketolase reaction were found. 
Intermediate high exchange rates of 0. 1-1 .7 for the ribulose 
5-phosphate epimerase, transaldolase, and fumarase/malate 
dehydrogenase reaction sequence were estimated, indicat- 
ing that exchange fluxes in these reactions were of the same 
order of magnitude as the net fluxes. 

The biological complexity of anaplerotic reactions and 
interconnections with decarboxylating reactions is demon- 
strated in various organisms. 20,22 As already stated in the 
biochemical model description (Fig. I), several enzyme re- 
actions carrying fluxes that connect phosphoenolpyruvate 
with oxaloacetate and pyruvate with oxaloacetate, or pyru- 
vate with malate, are present in C. glutamicum. 9 ,l 5 ,iq ' 3,4 ' 51 
However, since we wished to retain the limited set of as- 
sumptions as given in Table V, i.e., not speculate which 
anaplerotic enzyme plays the most important role in vivo, 
we regarded phosphoenolpyruvate and pyruvate as one net- 
work pool as well as malate and oxaloacetate. This ap- 
proach was also prompted by the measured identical frac- 
tional enrichments in the carbon atoms of phosphoenolpyru- 
vate and pyruvate (Table I). In all parameter estimates the 
estimates yielded a reproducible, positive value of 38.0% 
for the net anaplerotic flux of C3 units to C4 units. More- 
over, an exchange rate of 0.8 was consistently estimated, 
implying that the total forward in vivo flux of C4 unit gen- 
eration was 68.6% (1.8 times 38.0%) (while a back flux of 
30.6% (0.8 times 38.0%) simultaneously decarboxylated 
the C4 units oxaloacetate and/or malate]. This finding is in 
accord with the view that cycling occurs through the ana- 
plerotic reactions in mammals, 22,37 though to date no quan- 
titative information is available for bacteria. 

DISCUSSION 

Our approach to determine the fractional enrichments of 
individual carbons of precursor metabolites in the metabolic 
and isotopic steady state by amino acid analysis yields de- 
tailed information on the fluxes in the central metabolism. 
The method has crucial advantages over the methods con- 
ventionally applied for stationary flux determination using 



NMR or metabolite balancing. Zupke and Stephanopou- 
los 65 showed consistency between results obtained from 
metabolite balancing and information gained from 13 C 
NMR tracer data. While these authors conducted the two 
types of data analysis separately, in our approach the data 
can be analyzed simultaneously. 

In C. glutamicum MH20-22B we quantified the entry of 
the glucose 6-phosphate flux into the PPP as 66.4%, 
whereas 32.2% enters the glycolysis, at an absolute glucose 
consumption rate of 1.49 mmol/g dry weight h set as 
100.0%. Current data available on the relative use of both 
pathways in Corynebacterium species vary from 1 1 to 65%. 
In these cases the metabolic state differed greatly from our 
situation since washed cells 18 or strains producing histi- 
dine 17 or glutamate 60 were assayed. Moreover, cultivation- 
time-dependent variations of fluxes are documented for 
L-lysine 52 and L-glutamate synthesis 60 with Corynebacteri- 
um species, which we prevent in our continuous culture 
study. Data obtained from the early lysine production phase 
in a batch study based solely on metabolite balances may be 
most readily compared to our data with respect to the flux 
distribution and absolute flux rate at the glucose 6-phos- 
phate branch point. 59 However, the other fluxes reported in 
the latter study are entirely different from ours, probably 
due to the fact that in our continuous culture study glucose 
was completely consumed, whereas in the batch culture a 
large excess of glucose was present. 

In studies with hybridomas 27 ' 49 or bacteria 13,59 the re- 
duced forms of nicotinamide adenine dinucleotide (NAD/H) 
or NAD phosphate (NADP/H) balances were included since 
this information was required as an additional constraint on 
metabolism in order to be able to solve the flux network 
equation system. Since we do not need this information in 
our analysis, we can use our estimates to calculate the 
amount of energy and redox equivalents used and generated 
during metabolism. With respect to adenine triphosphate 
(ATP) and NADH this calculation is probably not useful for 
C. glutamicum, since for this organism the number of cou- 
pling sites in the respiratory chain is questionable. 23,51 Fur- 
thermore, in addition to the various possible anaplerotic 
reactions, other well-defined biochemical reactions have 
been identified that have the potential to operate in a futile 
cycle for ATP hydrolysis accompanied by proton translo- 
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pentose phosphate 




Figure 3. In vivo flux distribution in the central metabolism of the ly sine-producing C. 
glutamicum strain MH20-22B leuA resulting from the best fit to the experimentally deter- 
mined fluxes and l3 C labeling data. Numbers in oval symbols state the estimated net fluxes. 
Numbers in hexagonal symbols state the estimated exchange rates in reversible reactions 
expressed as multiples of the net flux. Al! net fluxes are expressed as a percentage of the 
glucose uptake rate (1.49 mmol/g dry weight h). 



cation. 59 However, for NADP/H the situation is different 
and our analysis allows an independent check on the bal- 
ance of its synthesis and requirement. 

It is documented that in C. glutamicum the PPP enzymes 
glucose 6-phosphate dehydrogenase 55 and 6-phosphoglu- 
conate dehydrogenase 56 are NADP linked. This is also the 
case for the isocitrate dehydogenase 10 of C. glutamicum. 
Thus, with the flux distribution shown in Figure 3 the total 
flux of NADPH synthesis is 194.8%. NADP is regenerated 
in the biosynthetic reactions with a stoichiometry (Table III) 
of 14,849 M-mol per 14,949 ^mol of glucose (i.e., 99.3% 
molar flux), as well as in lysine biosynthesis, with a flux of 



74.4% (Fig. 3, Table III). Thus, according to our calcula- 
tion, an excess of 194.8 - 99.3 - 74.4 = 21.1% molar 
flux relative to the glucose consumption rate is the result. 
This questions the use of a fixed NADP stoichiometry for 
flux determination. Instead, we have to conclude that C. 
glutamicum possibly regenerates NADP in vivo by an oxi- 
dase, or that alternative reactions are present whose opera- 
tion in vivo still has to be demonstrated. A carboxylating 
enzyme requiring NADPH (the malic enzyme) could in the- 
ory be a candidate. However, this seems extremely unlikely 
since the thermodynamic equilibrium favors the reverse re- 
action from oxaloacetate plus NADP to pyruvate plus C0 2 
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and NADPH. 54 Thus, at present the problem of the ob- 
served NADPH imbalance cannot be resolved. 

An important advantage over metabolic balancing is triv- 
ially and obviously that our method can handle bidirection- 
ality of fluxes where this is impossible with metabolic bal- 
ancing since it leads to singularities in the equation sys- 
tem. 58 We detected enrichment in C-4 of erythrose 
4-phosphate. Since the biochemistry dictates that this 
C-atom originates from C-5 of ribulose 5-phosphate, and 
ribulose 5-phosphate is not labeled via the oxidative part of 
the PPP, this is proof that an in vivo exchange of the glyc- 
eraldehyde 3-phosphate pool with that of ribulose 5-phos- 
phate occurs via the ribose 5-phosphate transketolase and 
epimerase reaction. Such exchange activities have been 
demonstrated in vitro in extracts of erythrocytes 4 and rat 
liver 11 by following the appearance of the 13 C label origi- 
nally present in an added sugar phosphate in new sugar 
phosphates of the PPP. Whereas in those studies exchange 
could only be investigated in isolated enzyme reactions in 
vitro, our approach allowed exchange reactions to be quan- 
titated as they are embedded in the complete network of the 
operative PPP in vivo. Since the PPP serves as a source of 
NADPH, ribose 5-phosphate as well as erythrose 4-phos- 
phate, all of which are required in different stoichiometric 
quantities for anabolism, restructuring of the excess gener- 
ated polyol-phosphates with glyceraldehyde 3-phosphate 
from the PPP and glycolysis is necessary to arrive again at 
the precursor metabolite phosphoenolpyruvate. 54 The ex- 
change rate carried by the ribose 5-phosphate transketolase 
reaction exceeds the net flux by about 30-fold. Intermediate 
exchange rates are present for the ribulose 5-phosphate epi- 
merase and the transaldolase, where the exchange fluxes are 
in the range of the net flux. 

From current data on the quantification of the exchange 
reactions of the separated enzymes based on in vitro mea- 
surements, 11 it is concluded that in various tissue prepara- 
tions the transketolase can exceed the net flux through PPP 
from 7- to 600-fold. Our finding of label mixing within the 
PPP also implies that using label in glyceraldehyde 3-phos- 
phate solely for the calculation of the PPP flux can lead to 
falsified net fluxes. The withdrawal of precursors for bio- 
synthesis is fully accounted for in our model, which is not 
the case in many NMR studies. The flux of ribulose 5-phos- 
phate and erythrose 4-phosphate withdrawal for anabolism 
also alters the fluxes and labeling patterns as already con- 
cluded by Walker et al. in the study on glutamate produc- 
tion with the coryneform Microbacterium ammoni- 
aphilum. 60 Disregarding the withdrawal for anabolism, we 
would only calculate a value of 58.7% entering the PPP. 
Further high exchange rates are determined in the glycolysis 
for the reactions involved in converting glyceraldehyde 
3-phosphate into pyruvate, which is in accord with full re- 
versibility of these reactions as quantified within erythro- 
cytes. 28 

In recent attempts to quantify the in vivo fluxes in C. 
glutamicum it had to be assumed that the glyoxylate cycle is 
inactive, since due to the restricted information on metab- 
olite balances and/or label incorporation, calculations of 



fluxes would otherwise be impossible. A quantification of 
the in vivo glyoxylate cycle flux for E. coli during growth 
on acetate was done by Walsh and Koshland. 61 Attention 
was drawn here to the fact that knowledge of C-l labeling 
in glutamate is essential to differentiate between the TCA 
cycle and glyoxylate cycle flux. These authors point out that 
flux information is also presenj in the amount of C0 2 re- 
lease, which correlates with TCA cycle activity. In our 
analysis the fractional enrichment in C-l of glutamate was 
not directly determined since it carries no protons. How- 
ever, this C-atom corresponds to the C-4 atom of oxaloac- 
etate (which became available by the quantitation of the 
enrichments in threonine und aspartate). This is incorpo- 
rated in our analysis as well as the quantitation of the C0 2 
labeling. We always found a small but reproducible nonzero 
value for the glyoxylate cycle flux. Apparently, the low 
activity of the isocitrate lyase determined in vitro (0.01- 
0.02 |xmol/mg protein min) 40 and of malate synthase 41 is 
sufficient to carry the flux entry in the glyoxylate cycle 
when C. glutamicum grows with glucose as substrate. 
Those in vitro activities are in the range of the ct-ketoglut- 
arate dehydrogenase activity (0.02-0.03 ^mol/mg protein 
min), which in our study carries a flux of 53.7%, or 0.013 
p,mol/mg dry weight min. Walker already concluded for M . 
ammoniaphilum that the glyoxylate cycle may be active in 
vivo. 60 

It would in principle be possible to replenish the TCA 
cycle via the glyoxylate cycle, but this is ruled out by the 
low glyoxylate cycle activity. Therefore, carboxytating ana- 
plerotic reactions are required to replenish the loss of TCA 
cycle metabolites for anabolic processes. The metabolite 
balance dictates that this anaplerotic net flux in C. glu- 
tamicum MH20-22B is 38.0%. This largely explains the 
labeling of a-ketoglutarate in C-2. In the anaplerotic reac- 
tions, it cannot be differentiated whether a carboxylation of 
phosphoenolpyruvate or of pyruvate occurs, since we re- 
garded these metabolites as one pool in our analysis, as we 
also did for malate and oxaloacetate. Which enzyme(s) 
carry the anaplerotic flux in vivo is currently a matter of 
much debate. 20,34 

It is questionable whether the phosphoenolpyruvate car- 
boxylase plays a major role in replenishing the TCA cycle, 
as originally concluded, 53 since a defined deletion mutant 
of this enzyme was not altered in L-lysine production, 14,34 
a situation in which there is a high demand for oxaloacetate. 
At this point, however, attention should be drawn to the fact 
that any flux analysis by deletion mutants is questionable 
since the flux is severely disturbed and the cell may react 
very flexibly. This is shown within the central metabolism 
of C. glutamicum, not only for phosphoenolpyruvate car- 
boxylase, but also for glutamate dehydrogenase 5 and pyru- 
vate kinase. 15 To quantify the in vivo fluxes of anaplerosis 
in the undisturbed situation, a detailed NMR study, such as 
an isotopomer analysis, would be particularly instrumental, 
as was used for the metabolism in the heart, where it was 
found that the carboxylation of pyruvate mainly contributes 
to the replenishment of the TCA cycle. 26 

It was pointed out as early as 1982 60,61 that the C-5 
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carbon of glutamate (which is directly derived from pyru- 
vate C-2) becomes 13 C labeled by a back flux from oxalo- 
acetate to acetyl coenzyme A. This back flux was ignored in 
subsequent studies of fluxes in C. glutamicum because C-5 
of glutamate was not included in the equations for tracer 
data analysis. Our study fully accounts for this important 
futile cycling and conclusively reveals that there is a back 
flux from the C-4 pool (oxaloacetate and malate) to the C-3 
pool (phosphoenolpyruvate and pyruvate) amounting to 
30.6%. According to the literature, 60 this back flux is re- 
flected in the presence of 13 C label in C-5 of a-ketoglut- 
arate, which is necessarily derived from C-2 of pyruvate 
since the sole synthesis of a-ketoglutarate results from py- 
ruvate (acetyl coenzyme A) plus oxaloacetate. Indeed, flux 
variation studies showed that the U C labels in pyruvate C-2 
and oxaloacetate C-2 reacted significantly to the exchange 
rate of the anaplerotic flux. At infinite exchange rates, these 
two labels would in fact become equal. Thus within C. 
glutamicum not only are high exchange rates present within 
glycolysis, or PPP (simply for thermodynamic reasons), but 
also within the complex sets of reactions exchanging carbon 
between phosphoenolpyruvate, pyruvate, oxaloacetate, and 
malate. It is speculated that such futile cycling serves to 
waste ATP. 37 It could also serve as a peculiar regulatory 
device, since in the presence of two regulated enzymes, one 
catalyzing the negative and one the positive flux, a regula- 
tory signal has a stronger effect on the net flux than if only 
the enzyme exclusively catalyzing the positive flux were 
present. 54 A flux from pyruvate to oxaloacetate was also 
quantified in rats, being as in our case of approximately the 
same size as that from oxaloacetate to phosphoenolpyru- 
vate. 22 

Our observed labeling data, especially the enrichment in 
ot-ketoglutarate C-2, could in principle also be partly ,ex- 
plained by the presence of channeling through the reactions 
catalyzed by succinate thiokinase, succinate dehydroge- 
nase, and fumarase (directed transfer of C-2, -3, -4, -5 of 
a-ketoglutarate to C-l, -2, -3, -4 of oxaloacetate), which 
seems to occur in eukaryotic cells. 35 However, the labeling 
data are in full accord with our current model based on the 
assumptions given in Table V, which includes full label 
scrambling through the TCA cycle reactions that are as- 
sumed to go without channeling. Based on this model, we 
also determined that only a small fraction of the oxaloace- 
tate formed by the anaplerotic reactions is subjected to label 
scrambling by back flux to fumarate and subsequent recon- 
version to oxaloacetate. If ever it should be confirmed that 
channeling is present in C. glutamicum, our data would be 
consistently interpreted by an increased exchange of the 
anaplerotically formed oxaloacetate with the symmetrical 
intermediate fumarate. The presence of back fluxes greatly 
influences the 13 C label distribution. For instance, it has 
been reported 61 that the flux through the isocitrate dehydro- 
genase derived from l3 C label data differed by more than 
10% depending on whether complete equilibration of the 
oxaloacetate and fumarate pools was assumed or no equil- 
ibration at all. This illustrates how, in most cases, simpli- 
fied models are used 46 that address the situation where ei- 



ther no back flux or complete equilibration is assumed. In 
contrast, our model, by virtue of its comprehensiveness and 
the large amount of experimental data, actually allows the 
in vivo exchange rates to be estimated from measured data. 

In summary, the approach presented here is a straight- 
forward method for deriving detailed information on meta- 
bolic fluxes in growing cultures. Many NMR studies con- 
ducted thus far used labeling models that describe only a 
limited part of the central metabolism 27, 60,61 (especially the 
TCA cycle and anaplerotic reactions) and, even so, did not 
fully take into account all possible carbon transitions, as 
stated by the authors themselves. 60,63 Consequently, these 
models could not be generalized 11,17,18,49 to other meta- 
bolic systems. 

In contrast, our equation system is fully comprehensive 
of all carbon transtitions and the mathematical formulation 
can be generally applied to other systems. Our approach 
offers the advantage of well-defined physiological condi- 
tions. We perform measurements under isotopic steady state 
of intracellular metabolites, whereas in batch fermentation 
conditions this is generally not the case, even after a 10-h 
period of incubation with the labeled substrate. 60 It should 
-prove useful to combine genetic analysis with the analysis 
of consequences on intracellular fluxes. This is of special 
interest for those mutants of the central metabolism where 
no balance disturbances occur. Furthermore, our approach 
provides the basis for comparing well-defined isogenic 
strains of C. glutamicum, which display great differences 
with respect to amino acid overproduction. 
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APPENDIX: MATHEMATICAL DETAILS OF 
FLUX DETERMINATION 

Intracellular flux determination under stationary conditions 
as applied in this contribution is based on a general model 
quantitatively explaining the relations between fluxes and 
13 C labels. This mathematical model, the automatic gener- 
ation of model equations, and the data analysis procedure 
are now briefly described to give an understanding of the 
theoretical background. 

Modeling Frameworks for Carbon Isotope 
Labeling Systems 

Carbon isotope labeling studies fall within the category of 
tracer experiments. General models and tools for tracer 
analysis were already developed in the 1970s; an excellent 
textbook is that by Anderson. 2 However, there are some 
new aspects in stationary carbon isotope labeling systems 
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that have not been explicitly taken into account before. In 
particular, the fluxes of all carbon atoms participating in the 
same reaction step are closely linked by an identical flux 
rate which considerably reduces the number of indepen- 
dently varying variables in the system. Furthermore, as will 
be shown below, reversible reaction steps can significantly 
influence the label patterns in metabolic intermediates. Al- 
though this has been described previously, 2 ' 46 metabolic 
exchange fluxes are considered in their complete generality 
for the first time in our approach. Moreover (as in refs. 64, 
65) a computer-supported framework for modeling, simu- 
lation, and data analysis is supplied. 

A Simple Example 

Before formulating a general (but rather abstract) mathe- 
matical model, the equations will be introduced for the case 
of the very simple example network shown in Figure 4. It 
was designed to demonstrate the effect of a bidirectional 
reaction step with a minimal number of carbon atoms. Al- 
though oversimplified, the example may be interpreted as a 
prototype of two branching pathways with different fates of 
carbon atoms. For example, in the central metabolism a 
comparable branching situation is encountered at the level 
of pyruvate. The bidirectional flux v 2 in Figure 4 then cor- 
responds to the anaplerotic flux between pyruvate and ox- 
aloacetate while the unidirectional flux v 3 represents the 
alternative transformation of pyruvate via acetyl coenzyme 
A through the citric acid cycle where carbon atoms change 
positions. The example motivates how in principle the in- 
formation on exchanging fluxes can be obtained from la- 
beling data. 

In Figure 4 the substrate uptake v, as well as the frac- 
tional labels in carbon atoms 5, and C 2 are assumed to be 
measurable. The intracellular flux v 2 takes place in both 
directions while v 3 is assumed to be irreversible. Finally, 
fluxes across the system borders (here v,, v 4 ) must always 
be assumed to be irreversible. 




Figure 4. Simple example network. Left: metabolite view. Right: carbon 
atom view. The extracellular flux v, and the labels B lt C 2 are assumed to 
be measurable. Input labels A t , A 2 are known. 



The system will be described in terms of fluxes v ( and 
labels Xj. Both flux directions of a bidirectional reaction step 
must be accounted for. In the example the forward and 
backward fluxes of v, are represented by the symbols v~*, 
vp. The irreversibility assumptions are then represented by 
= 0, = 0, - 0 while the corresponding net 
fluxes are vj** = vj* — v*~. 

The fractional labels at each carbon atom position of 
intracellular metabolites are denoted by B u B 2> C,, C 2 6 
[0,1]. Additionally the known input to the system is spec- 
ified by A lf A 2 . For example, a substrate A, which is 100 % 
labeled only at the first carbon atom position, is described 
by A! = 1, A 2 = 0. Clearly, the terminal nodes D t ,D 2 are 
not required for flux balancing. 

The metabolite net fluxes from and to an intracellular 
metabolite pool must add up to zero in the stationary state. 
This gives rise to the following linear metabolite flux bal- 
ance equations at stages B and C. 

B: v{* + v 2 *~ - v^ + 
C: + vr = + vj* u; 

from which follows 

v7 = v T and v T = v T - v r + v r = V - v T x - 

(2) 

This leaves v{~, v^* to be determined (remembering that v"p 
= v" et is measured directly). To this end the carbon isotope 
balance equations are established. As an example the flux 
vj*: £f J — > C 2 (interpreted as a molar carbon atom flux rather 
than a molar flux of metabolites as shown in the left and 
right part of Fig. 4, respectively) carries the amount of 
• B x labeled carbon atoms per unit time. This term con- 
tributes to the label balances of B x and C 2 , but with opposite 
signs. The complete set of equations for the pools B lf B 2 , 
C,, C 2 is now 

B x \ 0 = vfA, - v^B, + v^Cj - v^B x 

B 2 : 0 = v7A 2 - v?B 2 + v\~C 2 - v^B 2 

C,: 0 = v7B, - vrC, + vfB 2 - v?C? } 

C 2 : 0 = v?B 2 - v\-C 2 + v^B, - v?C 2 . 

From these equations it becomes clear that negative flux 
values do not produce a meaningful system description 
which has to be taken into account for flux determination. 

In the example the flux determination problem can be 
solved explicitly. We have finished when v^, have been 
computed from vj*, B lf C 2 . The label conservation equa- 
tions 

A x + A 2 - + B 2 = C, + C 2 

can be derived from Equations (2) and (3), Using this, the 
unknown fluxes are calculated to be 

V * A! + A 2 - B x - C 2 ' V{ 
Ci~ B x 

V2 = A I+ A 2 -2*/ (V2 +Vl) - 
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This shows how both rates of bidirectional reactions can in 
principle be computed from labeling data. With conven- 
tional measurement techniques nothing more than the quan- 
tity vf = = vf - vj~ + v 3* can be observed! 

General Flux Model 

A general model for arbitrary metabolic networks is now 
formulated. Based on this model, fractional labels can be 
computed when all fluxes are known. This enables 13 C 
isotope labeling experiments to be simulated. The model 
makes extensive use of matrix calculus to obtain a concise 
notation. All matrices are constant and represent structural 
qualities of the underlying network, 38 i.e., they do not de- 
pend on the state variables. Furthermore, they are more or 
less sparsely populated, which can be exploited for reduc- 
tion of memory requirements and speedup of numerical 
computations. 

To begin with, the state vectors have to be introduced. 
The flux vectors 

(of equal dimension) comprise both directions of all the 
molar fluxes in the reaction network. In the example of 
Figure 4 these fluxes are v~* = (vf, vf , vf, vf) r and v*~ 
= (vi~> v 2~. v 3~> V D T - From tnese vectors the composite 
overall flux vector (of twice the dimension) and the corre- 
sponding net flux vector is formed: 



v = 



- v 



This is how in our analysis full account can be taken of 
bidirectional fluxes. On the other hand, it will be explained 
in the section on exchange fluxes how available knowledge 
can be incorporated by introducing irreversibility or rapid 
equilibrium assumptions. 

Finally, the vector of intracellular fractional labels of 
each (enumerated) carbon atom in the intracellular metab- 
olites as well as the vector of input labels comprising the 
known fractional labels of all carbon atoms fed into the 
system are required. These two vectors are denoted by 



x and v; i.e., all terms of type ±x ; v t containing a certain 
flux v, and some involved pool variable x;. The correspond- 
ing coefficients ± 1 can therefore be collected in a square 
matrix P, and a rectangular matrix Pj n associated with v ( . 
Here, P, describes transitions between intracellular pools, 
while P} n describes transitions involving one extracellular 
partner. In the example [compare to Equation (3)] 



/-. • • 

-1 ■ 



P7 = 



l 



pr = 



-i 



-i 



M 



\ / 

Using matrix notation the carbon pool balance equations 
can now be written as 



2 v ' p < 



X + 



2 v, • P| n \ • x in - 0. (6) 



From these equations it follows immediately that 

x = X (v) = v, * pA- 1 /2 v < • p i' 



• l> ln \ * v~ M 



(7) 



is a function of v. This equation is the key to the numerical 
simulation of 13 C tracer experiments. 



In the example x = (5„ B 2 , C,, C 2 ) r and x in = (A,, A 2 ) T . 

The metabolite flux balances are formulated using the 
well-known stoichiometric matrix N and Equation (4), 
which leads to the linear equation system 



N • v net = N 



- N v^ = 0. 



(5) 



In the example equations (1) lead to 



N = 



1-1-1 

• 1 1 -1 



where points denote zeroes. The carbon isotope label bal- 
ance equations (3) have a bilinear structure with respect to 



Flux Estimation 

Clearly, flux estimation is an inverse problem associated 
with Equation (7); i.e., a flux vector v has to be found that 
reproduces all the measured fluxes and labels. Unlike the 
example, the flux determination problem cannot always be 
solved explicitly because the nonlinear equation system is 
too greatly entangled. In this situation a parameter Fitting 
approach is suitable, as is now described. 

The measurement matrices M" ct , M x are used to describe 
which net fluxes and intracellular labels are measured. They 
are both composed of unit vectors. In the example v" et , B lt 
C 2 are measured, and therefore we have 



MARX ET AL: QUANTIFICATION OF FLUXES IN CENTRAL METABOLISM 



125 




The complete measurement equations are then 

w = M; ct • v net + e v 

y = M x x + e x . (8) 

Here w, y denote the vectors of measured fluxes and labels; 
e v , e x are normally distributed measurement noise terms 
with expectation 0 and covariance matrices 2 V , X K . 

The familiar least-squares flux estimator is constructed 
from Equations (8) by minimizing the sum-of-squares func- 
tion 

k(v) = |[w - Mf • (v- - OIlL + |jy ~ M x • x(v)||| x 

(9) 

where H£|[| = £ T ■ 2P 1 • f; denotes the squared weighted 
norm corresponding to a covariance matrix 2. Putting equa- 
tions (5), (7), and (8) together we are left with a linearly 
constrained minimization problem for flux estimation: 

minimize k(v) 

with respect to N • v neI = 0 and v^O. 

(10) 

The resulting flux estimator is denoted by v. Numerical 
algorithms for solving this problem are described later. 

Exchange Fluxes 

In the example from Figure 4 only one intracellular flux was 
assumed to be bidirectional as expressed by the flux vari- 
ables v^" and vp. If v 3 would also split into two directions, 
the available amount of measured quantities would no 
longer be sufficient for flux determination. In general, the 
number of unknown flux variables must not exceed the 
number of measured values. 

To prevent the investigated labeling system from being 
underdetermined, reversibility assumptions on fluxes must 
be introduced. Three types of assumptions are frequently 
encountered in the literature: 

1. If a reaction step is known to be highly reversible (i.e. 
vp 4- vp > v" et ), we can eliminate some variables by 
using a rapid equilibrium technique. 46 

2. When reaction steps cross the system border (e.g., a 
substrate uptake step), they are assumed to be irrevers- 
ible. 

3. Reaction steps exposing a large AG 0 value in vitro may 
be assumed to be irreversible in vivo too, which means 
that either vp or vp is zero. 

A formal notation is needed to express these assumptions, 
i.e., to introduce further numerical constraints on the flux 
vector v in the minimization equation (10). To this end a 
quantity must be defined that describes the degree of re- 



versibility for some bidirectional reaction step. For exam- 
ple, the quantity v 7* + vp was introduced elsewhere 46 to 
study the rapid equilibrium behavior vp + vp -* oo. How- 
ever, this is not a useful formalism when irreversibility has 
to be expressed since in that case it depends on v" e \ For this 
reason we define the exchange flux of the ith reaction to be 

vf x = min(v^, vp), (11) 

which means that vf* is the amount of flux that is common 
to both reaction directions in relation to v" ct . This makes it 
easy to express the above-mentioned modeling assump- 
tions. 

v" = 0 means irreversibility of step i 

v" = oc means rapid equilibrium of step /. 

Clearly, any other value for v" ct can be fixed for simulation 
purposes. On the other hand, if nothing is kept fixed, both 
flux directions can be freely varied. The reversibility as- 
sumptions made for a given metabolic network must be 
carefully documented to ensure reproducibility of the re- 
sults. In our case they are given within Figure 3 where 
bidirectional fluxes are indicated by hexagonal symbols 
while all other fluxes are assumed to be irreversible. For 
better readability the numbers given in the hexagons are the 
normalized exchange fluxes |vf x /v" ct |. 

In previous approaches, exchange fluxes could not be 
taken into extensive consideration because not enough mea- 
surement information was available to quantify fluxes and 
exchanges simultaneously. As a consequence, previously 
obtained results may be biased due to oversimplifying as- 
sumptions concerning irreversibility or rapid equilibria that 
had to be made. The experimental technique presented in 
this contribution supplies a previously unattainable amount 
of (stationary) flux and label measurements that enables to 
quantitate both net fluxes and exchange rates with only a 
few assumptions. 

Numerical Algorithms 

Presentation of the rather technical details of incorporating 
additional constraints and their computational treatment ex- 
ceeds the scope of this article. Thus only some basic prin- 
ciples will be given here and more details will be published 
elsewhere. 

It is well known from general tracer kinetics 2 that the 
system matrix 2 v,. * P,. of Equation (6) is diagonally dom- 
inant and thus can be treated by iterative methods. 15 On the 
other hand, it can be proven that high exchange fluxes lead 
to ill-conditioned system matrices, so that QR decomposi- 
tion far outperformed an iterative solution for solving Equa- 
tion (7). 16 ' 64 

Having established a simulation algorithm, the flux esti- 
mation equation (10) can be solved by iterative optimization 
methods. 47 The linear constraints given by Equation (5) are 
treated by singular value decomposition. 36,58 It should be 
recognized that the inequality v ^ 0 has to be strictly rec- 
ognized and that exchange fluxes as defined by Equation 
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(11) are not everywhere different iable with respect to v. 
This makes flux estimation a delicate numerical problem. 
As a first practicable approach, a robust derivative- free 
method was used to solve the minimization problem at the 
cost of high computation times. Further work is devoted to 
the development of a modified Gauss-Newton type algo- 
rithm that enables the almost everywhere differentiable, lin- 
early constrained problem to be treated. 

The most important property of exchange fluxes with 
respect to the parameter fitting algorithm is that a fixed 
exchange flux does not prescribe a reaction direction; i.e., 
the sign of the corresponding net flux can still vary freely. 
For this reason, the algorithm will find the reaction direc- 
tion that best fits the data. 

Automatic Generation of Equations 

As shown in this Appendix all model equations required for 
describing the isotope labeling system can be built up from 
certain vectors and matrices. However, the dimensions can 
become quite large (about 100 for a reasonably complex 
metabolic network). A manual input of these structures 
would be time consuming and susceptible to errors. More- 
over, model variation studies would be quite time consum- 
ing if manual input were used. For these reasons a compiler 
program was implemented that generates the required struc- 
tures from a text file describing the biochemical network 
together with the measured data in a convenient way. 62 The 
following formal notation is used therein: 

1 . Each reaction is expressed using the familiar biochemi- 
cal sum notation. The atom transitions are specified af- 
terward. For example, the transaldolase reaction from 
the PPP is written as 

x9 : gap + s7p > e4p + f6p 

#ABC + #abcdefg > #defg + #abcABC . 



This means that the first carbon atom of gap (denoted A) is 
taken over to the fourth carbon atom of f6p and so on. 

2. It now becomes clear how integer stoichiometric coef- 
ficients other than 1 can be expressed. If a metabolite 
occurs twice, it has to be repeated. For example, the 
conversion of fructose 1,6-bisphosphate to glyceralde- 
hyde 3-phosphate is simply written as 

x3 : f16bp > gap + gap 
#ABCDEF > #CBA + #DEF . 

Similarly the scrambling reaction of fumarate can be ex- 
pressed by 

x14 : fum + fum > oaa + oaa 

#ABCD + #U VWX > #ABCD + #XWVU . 

3. Sometimes noninteger coefficients have been used to 



express biomass composition. Such equations make no 
sense when carbon atoms have to be traced through the 
network. However, they can be easily replaced by a set 
of equations that describe the incorporation of each pre- 
cursor metabolite into biomass. For example, incorpo- 
ration of glucose 6-phosphate is written as 

bl:g6p > g6pBiomass 
#ABCDEF > #ABCDEF . 

The corresponding "noninteger coefficient' 1 is then sup- 
plied separately as a flux measurement value computed 
from biomass composition and growth rate. This is in any 
case the more natural way to take biomass composition 
into account. 

The NMRFlux program translates this text input to the vec- 
tor and matrix notation described above by using the recur- 
sive descent method as a standard technique for compiler 
construction. 1 

A companion program takes the generated vectors and 
matrices as input and performs the numerical computations 
as described above. For performing simulation runs, all 
linear independent flux parameters can be interactively 
changed. Similarly, the exchange fluxes can be fixed to 
given values. This led to valuable insight concerning the 
sensitivity of estimated exchange fluxes. Finally, the itera- 
tive optimization procedure can be set up and all results are 
written to a protocol file for further use. 
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Two glycerol utilization mutants of Mycobacterium smegmatis that were unable to utilize most carbon sources 
except glucose were isolated. Supplementation of these media with small amounts of glucose restored growth 
in the mutants; these strains are therefore glucose autotrophs. The mutant phenotype is complemented by the 
gene encoding phosphoglucose isomerase (pgi), and direct measurement of enzyme activities in the mutants 
suggests that this gene product is absent in the auxotrophic strains. Mapping of the mutant allele by Southern 
analysis demonstrates the presence of a 1-kb deletion extending into the coding sequence of pgi. The possible 
roles of phosphoglucose isomerase in mycobacterial cell wall synthesis and metabolic regulation are discussed. 



The resurgence of diseases caused by mycobacteria has 
spurred a new interest in these acid-fast organisms. Although 
current research is largely directed towards virulence mecha- 
nisms in pathogenic strains of mycobacteria, basic biochemical 
information about these organisms is essential for understand- 
ing their interactions with their hosts and for developing new 
therapeutics. Knowledge of intermediary metabolism of the 
mycobacteria has been fragmentary, in large part due to the 
low growth rates, tendency to clump in culture, and lack of 
genetics systems of members of this genus. Mycobacteria ex- 
hibit unusual and complex patterns of carbohydrate utilization. 
In Mycobacterium tuberculosis, Mycobacterium phlei, and My- 
cobacterium smegmatis, glycerol utilization is blocked until all 
other carbon sources, including glutamate and asparagine, 
have been consumed (3, 23, 38). Nevertheless, glycerol is con- 
sidered the preferred carbon source under laboratory condi- 
tions. Growth of most species on this substrate results in 
greater mass per unit volume of growth medium (32). 

Glycerol and its metabolic derivatives are at the interface of 
catabolic (glycolytic and oxidative pathways) and anabolic (lip- 
id biosynthesis and gluconeogenesis) processes and therefore 
play important growth regulatory roles in cells grown on a 
variety of carbon sources. In addition, glycerol is an important 
precursor of key cell wall constituents. The presence of glyco- 
lipids such as the acyl glycerols and acylglucoses varies widely 
depending on the carbon source in the medium (4). A number 
of the enzymes of glycolysis, the hexose monophosphate shunt, 
and the Entner-Doudoroff pathways have been studied (re- 
viewed in reference 26). However, carbon metabolism in my- 
cobacteria has not been characterized at the genetic level and 
few genes encoding glycolytic enzymes have been cloned. 

We have initiated a genetic approach to the study of carbo- 
hydrate metabolism in M. smegmatis by isolating and charac- 
terizing mutants unable to utilize certain substrates as sole 
carbon sources. In this study, we describe two strains obtained 
in a screen for glycerol utilization mutants. Unexpectedly, 
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these mutants are unable to use any carbohydrate except glu- 
cose as a sole carbon source. Genetic complementation of the 
growth defect on glycerol was accomplished by expression of 
the gene encoding phosphoglucose isomerase (pgi), which we 
have cloned and sequenced. The mutants lack detectable levels 
of Pgi, a key enzyme in glycolysis that functions at the juncture 
of gluconeogenesis and catabolism. Finally, we demonstrate 
that pgi mutants in M. smegmatis are defective in the utilization 
of many carbohydrates as sole carbon sources because they are 
actually auxotrophic for glucose. 



MATERIALS AND METHODS 

Bacterial strains and plasmids. The bacterial strains and plasmids used are 
shown in Table 1. M. smegmatis mtrl55 (27) was used for all mycobacterial 
experiments. Escherichia coli DH5ot was used to propagate DNA for cloning, 
sequencing, and transforming M. smegmatis. mc 2 743 and mc 2 744 were derived in 
this study. The plasmids used were pYUB18 and a cosmid library contained 
therein (13), pMVI203 (28), and pDT252, created in this study, which is 
pMVI203 bearing a l,9-kb fragment of mcMSS DNA encoding the pgi gene. 

Media. M. smegmatis was grown as described previously (6). The complete 
medium used was Middlebrook 7H9 (liquid) or 7H10 (agar; Difco). Minimal 
media were composed of basal salts medium (0.1% KH 2 P0 4 , 0.25% NaHP0 4 , 
0.5% NH 4 CI, 0.2% KjSCX,) supplemented with either glucose (1.0%) or glycerol 
(0.5%). For auxotrophy studies, glucose was added at a concentration of 0.01%, 
a level insufficient to support growth in the absence of another carbon source. All 
liquid cultures were supplemented with 0.05% Tween 80 (Sigma Chemical Co.). 
Kanamycin sulfate (Sigma) was used at 50 ^.g/ml. 

Mutant isolation. M. smegmatis (mc 2 155) was mutagenized by treatment with 
A^-mcthylW-nitro-Af-nitrosoguanidinc (MNNG; Sigma) at 100 u.g/ml in basal 
salts medium. After treatment at 30°C for 2 h, the cells were pelleted and washed 
three times in basal salts medium with 0.05% Tween 80. The cells were resus- 
pended in Middlebrook 7H9 medium with glycerol at 5 g/liter and glucose at 10 
g/liter and incubated at 37°C for 3 h to allow expression of mutant gene products 
and then diluted onto minimal plates containing glucose and glycerol with 0.05% 
Tween 80 to facilitate replica plating. Colonies were replica plated onto minimal 
glycerol plates, and those that did not grow on glycerol were collected from the 
master plates. 

Enzyme assays. Cell extracts were prepared by growing cells in the designated 
minimal medium to an of 0.8 to 1.0. Cells were harvested by centrifugation 
at 3,000 X g for 10 min at 4°C and washed in 10 mM potassium phosphate buffer 
(1 mM EDTA and 1 mM dithiothrcitol). The cell pellet was resuspended in an 
equal volume of buffer, and the cells were disrupted by passage twice through a 
French press at 18,000 lb/in 2 . The extracts obtained were diatyzed overnight at 
4°C in the same phosphate buffer and stored at -70°C. 

Glucose phosphate isomerase activity was determined spectrophotometrically 
by measuring the decrease inAy^ caused by reduction of p-NADP. Reactions of 
a total volume of 3 ml contained 42 mM glycylglycine, 3.3 mM D-fructose-6- 
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TABLE 1. Strains and plasmids used 



Strain or plasmid 



Description 



Reference 



Strains 
E. coli K-12 DH5a 



M smegmatis 
rac 2 155 
mc 2 743 
mc 2 744 

Plasmids 
pYUB18 
pMVI203 
pDT252 



endAl hsdR17(xy~ m K + ) 
supE44 thi-1 recAl 
gyrA(NaV) relAl 
b(lacIZYA-arxF)U169 
deoR [4>80d/ac 
A(/acZ)M15] 



ept-l 

ept-1 M(pgi) 
ept-l pgi2 



12 



27 

This work 
This work 



24 
28 

This work 



phosphate, 0.67 mM p-NADP, 3.3.mM MgCl 2) and 5.0 U of glucose-6-phosphate 
(G6P) dehydrogenase (8). 

Triose phosphate isomerase activity was determined spectrophotometrically 
by measuring the decrease in A 3M caused by reduction of p-NADP. Reactions of 
a total volume of 3 ml contained 280 mM triethanolamine, 245 mM glyceralde- 
hyde-3-phosphate, 0.132 mM p-NADP, and 4 U of a-glycerolphosphate dehy- 
drogenase (19). 

a-Olycerol phosphate dehydrogenase activity was determined spectrophoto- 
metrically by measuring the decrease in A 340 caused by reduction of p-NADP. 
Reactions of a total volume of 3 ml contained 290 mM triethanolamine, 2.5 mM 
dihydrafyacetone phosphate, and 0.13 mM p-NADP (14). 

Complementation of mutants with wild-type genomic libraries. The carbohy- 
drate utilization mutant strains mc 2 743 and mc 2 744 were transformed with a 
cosmid library prepared from mc 2 155 (wild-type) DNA. This library was con- 
structed by ligating 32- to 40-kb fragments of genomic DNA into pYUB18 (24). 
Transformants were plated on glucose minimal medium to estimate the trans- 
formation efficiency and on glycerol minimal medium to select complementing 
cosmids. mc 2 743 was complemented for growth on glycerol by 12 cosmids (of 3 x 
Iff screened), and mc 2 744 was complemented by 7 cosmids (of 6 x 10 3 screened). 
Cosmid DNA was isolated from two representative Glp + colonies of strains mc 2 743 
and mc 2 744 and amplified in E. coli DH5a. These four cosmid DNAs were then 
used to transform mc 2 743 and mc 2 744. Each was fully capable of complementing 
mc 2 743 and mc 2 744 for growth on glycerol minimal plates. Two of the cosmids 
were then partially digested with Sau3A, and fragments in the range of 1,500 to 
5,000 bp were isolated from an agarose gel. The Saw 3 A fragments were ligated 
into the BamUl site of pMVI203. The resulting recombinant plasmids were 
pooled and used to transform mc 2 743 and mc 2 744. Nine plasmids were recov- 
ered. All nine complemented the glycerol-negative phenotype of mc 2 743; mc 2 744 
was complemented fully by four plasmids. The four plasmids which comple- 
mented the glycerol-negative phenotype of both mc 2 743 and mc 2 744 showed 
identical restriction maps; the smallest of these, pDT252, bearing a 1.5-kb insert, 
was chosen for sequencing. That five plasmids did not complement both mutant 
strains suggests that mc 2 744 may have an alternate or additional lesion, and this 
will be investigated further. 

DNA sequencing and deletion mapping. The Erase-a-Base system (Promega) 
was used to generate a series of deletions from one end of the 1.5-kb insert in 
pDT252. The plasmid was digested with Noil and Kpnl and treated with exonu- 
clease III for 0.5 to 10 min. Samples of the reaction were removed at 30-s 
intervals and separated on an agarose gel. DNA fragments differing by 200 bp in 
length were collected, treated with SI nuclease and Klenow large fragment to 
introduce blunt ends, and self-ligated to produce a series of plasmids for se- 
quencing. Sequencing was performed with an automated sequencer (Applied 
Biosystems). 

Southern hybridization was performed to determine whether sequences found 
in the wild-type pgi region were missing in the mutant chromosome of mc 2 743. 
Two micrograms of genomic DNA prepared from mc 2 155 and mc 2 !^ was 
digested with various restriction enzymes and analyzed on agarose gels. The 
Southern transfer procedure used was described previously (15). The probe used 
for Southern analysis was the plasmid pDT252 containing the entire pgi gene. 

Nucleotide sequence accession number. The l,92S-bp sequence bearing the 
phosphoglucose isomerase gene of M. smegmatis was deposited in GenBank 
under accession number U88433. 



RESULTS 

Isolation and complementation of glycerol utilization mu- 
tants. Glycerol-nonutilizing mutants of M. smegmatis were iso- 
lated as described in Materials and Methods. Of 40,000 colo- 
nies screened, two isolates, mc 2 743 and mc 2 744, were unable to 
use glycerol as a sole carbon source. mc 2 743 and mc 2 744 had 
colony morphologies and growth rates on glucose minimal 
plates similar to those of the wild-type parental strain, mc 2 155 
(data not shown). 

Strains mc 2 743 and mc 2 744 were transformed with a cosmid 
library of mc 2 155 (wild-type) DNA (see Materials and Meth- 
ods). Several cosmids which complemented the glycerol defi- 
ciency in both strains were isolated. Analysis of restriction 
enzyme digests of four of the cosmids revealed overlapping 
inserts of similar sizes from the cosmids. An insert of 1.9 kb 
which complemented the deficiencies in strains mc 2 743 and 
mc 2 744 was subcloned from one of the complementing cos- 
mids. The insert was cloned into pMVI203 to create pDT252. 
Initial restriction mapping of the clone and construction of 
internal deletions indicated that the entire 1.9 kb was required 
for complementation of growth on glycerol for both mutants 
(data not shown). 

Characterization of the mycobacterial gene encoding phos- 
phoglucose isomerase. The sequence of the region of wild-type 
M. smegmatis DNA complementing the two glycerol utilization 
mutants is shown in Fig. 1. The fragment contains a large single 
open reading frame (ORF) on one strand and two smaller 
ORFs on the opposite strand (data not shown). The single long 
ORF encodes a putative protein of 442 amino acids, beginning 
with GTG (Val) at position 420, with a predicted molecular 
mass of 47,649 Da. An alternative start codon is found at 
position 471, with a methionine codon (ATG), predicting a 
polypeptide product of 425 amino acids or 45,818 Da. The 
codon usage pattern is consistent with that of other mycobac- 
terial genes, namely, a G+C content of >65%, due largely to 
G or C in the third position (36). To determine the identity of 
the ORF contained in the complementing fragment, we trans- 
lated the 442-codon ORF and subjected the derived amino 
acid sequence to a homology search in GenBank by using 
FAST A. The predicted M. smegmatis sequence shows striking 
homology to that of the products of genes encoding the en- 
zyme phosphoglucose isomerase (Pgi). Multiple sequence 
alignment demonstrated extensive amino acid homology be- 
tween the M. smegmatis protein and the amino acid sequences 
of Pgi from Homo sapiens (37), Saccharomyces cerevisiae (31), 
Leishmania mexicana (22), Bacillus stearothermophilus (30), 
and E. coli (11) (data not shown). For example, human Pgi has 
49.4% identity over a 447-amino-acid overlap, S. cerevisiae has 
50.7% identity over a 436-amino-acid overlap, and E. coli has 
52.9% identity over a 448-amino-acid overlap. 

The mutant allele of mc 2 743 was characterized by designing 
six pairs of primers to amplify the mutant gene in six overlap- 
ping fragments. It was observed that the 3' end of the mc 2 743 
allele was resistant to amplification by PCR. Southern analysis 
identified a deletion of approximately 1 kb in the 3' end of the 
pgi gene. Chromosomal DNA from wild-type and mutant cells 
was restricted with three enzymes (Ncol, Pstl, and Sphl) rec- 
ognizing internal sites suggested by the sequence of the cloned 
fragment carrying the wild-type gene (data not shown). In all 
cases, there is an approximately 1-kb reduction in the size of 
one band of the mutant allele; a restriction map of the region 
is shown in Fig. 2, This result is consistent with the reversion 
frequency measured for this mutant, which is <10 -8 (data not 
shown). Precise determination of the deletion endpoints in the 
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GGATCCACAGTCCGGTGGGG ACGOATGAAATCAGGTCTTC GATGCOCATACATCCATCAT GTACACCCTCGACACGTTCA CCGCACTAAGGTCGACACAA 100 

TGAGCGCTGACATCACCGAA AC TCCC GC ATGGCAG GCTTT GTCGGACCACCACGCCGAGA TCGGCGACCGGCATCTGACA GAACTGTTTQCTOACGATCC 200 

CGCGCGGGGAACCGAGTTGG CGTTGACCGTCGGGGATCTC TACATCGACTACAGCAAGCA CCGCGTCACGCGGCGCACCC TGGACCTGCTCGTCGACCTG 300 

GCGCGCGCGGCCGQGTTGGA GGAGCGCCCCGATGCGATGT TCGCCGGCGAGCACATCAAC ACCTCCGAGGACCGCGCGGT GCTQCACACCGCGCTOCOCT 400 

CCCOCGGGACGCCAAGCTG G TGG TGCACGGTCAGGACGTG GTCGCCGATGTGCACGACGT GCTGGACCGCftTGGGCGATT TCACCGATCGGTTGCGCAGC 500 

V VDOQDV VADVHDV L D R H G D F TDRIRS 
GGAGAGTGG AC CGGCGCCAC CGGCGAACGC ATCACC AC GO TGGTCAACATCGGCATCGGC GGCTCGGATCTGGGTCCGGT GATGGTGTACGACGCGCTGC 600 
GEWTGAT GERITTV VNIGIG GSDLGPV H V Y D A L R 
GCC ATTAC GCCGACGC CGGC ATCTCGGCGCGCTTCGTGTC CAACGTCGACCCCGCCGATC TGOTCGCCAAACTCGACGGC CTGGAGCCGGCGAAOACGTT 700 

HTADAO ISARFV3 HVDPADL VAKLDG LEPAKTL 
GTTCATCGTCGCCTCCAAGA CGTTCTCGAC GCTGGAGACG CTGACCAATGCGACCGCOGC GC GGC GCTGGCTC ACCG ACQ CCCTGGGTGACGCCOCGGTG BOO 

PIVASXT rSTLET LTHATAA R R W L T D A LGOAAV 
GCCAAGCATTCGTCGCQGTG TCCACCAACAAGAAGCTGGT CGACGAAGTTCGGCATCAAC ACCGACAACATGTTCGGGTT CTQGGACTGGGTGGGCGGGC 900 
AKH3SRC P P T R S W 3 TKTGIN T D N N F G F WDWVGOR 
GCTACTCGGTGGACAGTGCG ATCGGGTTGAGCGTGATGGC CGTGATCGGCAAGGAACGCT TCGCCGAGTTCCTGGCCGGG TTCCACATCGTCGACGAGCA 1000 

YSVDSA IGLSVMA V I G K E R F AEFLAG FBIVDEB 
CTTCCGCACCGCCCCGCTGC ACCAGAACGCTCCGGCGCTG CTGGGGCTGATCGGGCTGTG GTACTCGAATTTCTTTGGCa CGC AATC CCGTQCGOTQCTG 1100 

FRTAPLH Q K A P A L LGLIGLW YSMFFGA Q8RAVL 
CCCTACTCCAACGACCTGTC ACGGTTC GCGGC GTACCTGC AGC AGTTGAC CATGGAATC C AACGGC AAGTCG GT ACGCGC CGACGGCACACCGGTGTCCA 1200 
PYSHOLS R F A A Y L Q Q 1 T M E S N G K S V R A OGTPVST 
CCGACACCGGTGAGATCTTC TGGOGCG AGC CG GGC ACC AA CGGCCAGCACGCCTTCTACC AGTTGCTGCACCAGGGCACC CGACTGGTGCCCGCOGACTT 1300 

DTGEIF WOEPGTN GQBAFYQ LLHQGT RLVPADF 
CATCGOQTTCTCCCAGCCCA CCGACGACCTGCCCACGGCC G ATGGC ACC GGCAGCATGC A CGAC CTGCTGATG AGC AACT TCTTCGCC C AG ACCC AQOTG 1400 

IQFSQPT DDLPTA DGTGSHB DLLHSNF FAQTQV 
CTGGC ATTCGG AAAG ACCG C CGACGCAATAGCTTCCGAOG GCACCCCCGCAGATGTGGTG CCGCACAAGGTGATGCCCGG CAACCGGCCCACCACCAOCA 1300 
LAFGKTA DAIABEG TPADV V PBKVMPG HRPTT8I 
TCCTCGCCACCAAGCTCACC CCGTCGGTCGTCGGGCAGTT G ATCG C CCTCT ACGAAC ATC AGGTGTTCACCGAAGGCGTG ATCTGGGGCATCGACTCGTT 1600 

LATRLT PSVVGQL IALYEBQ VFTEGV IWG10BF 
CGACCAGTGGGGCGTCGAAC TCGGCAAGACGCAGGCCAAG GCCCTGCTCCCGGTTCTGAC CGGC GACAAATCTC CTGCCG C GC AGTCGGATAC GTCCACC 1700 

DQWOVEL GKTQAK ALLPVLT G D K SPAA QSDTST 
GACGCCCTGOTGCGGCGCTA CCGGACCGAGCGGGGCCGGC CCGCTTAGGCTTCGGCGTCQ GCTTCTGCGTCCGAACGCCT GGTCCAGCCCAGGATCATCG 18O0 
DALVRRY RTERGRP A. 

CGGGCAGGCAGCTGCCGATC ATC AAGGCCGC C AAGGCC AT CAGCGCACCGGCATAGGCCA GGTCGCGGATGGTCTTGGCG TCGGCGCCGGTTCAGCTOCG 1900 
TCGAC AC G ATG ATCC AGCTG CAGAATTC 1928 

FIG. 1. Nucleotide and deduced amino acid sequences of the EcoRl-BamHl fragment containing the phosphoglucose isomerase gene of M. smegmatis. Two putative 
start codons are underlined. 



pgi gene of mc 2 743 and characterization of the mutant allele of 
iTKr744 are under way. 

Characterization of pgi mutants. The enzyme Pgi exerts 
considerable control at a pivotal point at the juncture of three 
metabolic pathways, i.e., the Embden-Meyerhoff-Parnass, the 
Entner-Doudoroff, and the phosphogluconate pathways (Fig, 
3) (see Discussion). We examined the growth of the parental 
strain (mc 2 155) and the two glycerol utilization mutants 
(mc 2 743 and mc 2 744) on nine carbon sources, which is sum- 
marized in Table 2, Each of the three strains carries the E. 
coli- mycobacterium shuttle vector, pMVI203, containing a 
kanamycin resistance gene and a multiple cloning site (28). 

The noncomplemented mutant strains mc 2 743 and mc 2 744 
grow only on glucose and poorly on galactose: no other car- 



bohydrate tested (glycerol, fructose, succinate, Casamino Ac- 
ids, gluconate, sorbitol, or ribose) was able to serve as a sole 
carbon source for growth. This suggests that the metabolic 
requirement is G6P or some derivative thereof, since of the 
nine carbon sources tested, only galactose typically enters the 
glycolytic pathway above G6P (Fig. 3). Pgi, which catalyzes the 
conversion of G6P to fructose-6-phosphate, is the only source 
of G6P during growth on substrates other than glucose or 
galactose. 

We tested the growth of the wild- type and mutant strains 
carrying the complementing plasmid, pDT252, on these carbon 
sources. Table 2 shows that the wild-type and mutant strains 
grow equally well on glycerol minimal media in the presence of 
pgi carried on a multicopy plasmid (pDT252). This indicates 




lkb 

FIG. 2. Restriction map of the pgi region of M. smegmatis. Pstl (P), Ncol (N), and Sph\ (S) were used to analyze chromosomal DNA of mc^SS and 111(^743 by 
Southern analysis (see Materials and Methods). The sequenced region is indicated by an unshaded box. A horizontal arrow indicates the direction of transcription of 
the pgi ORF (shaded area). A hatched bar indicates the extent of deletion in mutant strain mc 2 743. 
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FIG. 3. Intermediary metabolic pathways in bacteria known to be present in mycobacteria (26). The routes to cell wall components are indicated. 



that the pgi gene is sufficient to restore wild-type growth abil- 
ities to the two mutants. 

To corroborate these phenotypic observations, levels of 
three glycerol and glucose metabolic enzymes were measured 
in extracts prepared from wild-type cells, the pgi mutants, and 
all three strains carrying the cloned copy of the pgi gene. We 
prepared extracts from wild-type and mutant cells carrying 
either the cloning vector pMVI203 or the pgi plasmid pDT252 
growing in either glucose or glycerol. Table 3 shows the levels 
of Pgi in wild-type cells and in mc 2 743 and mc 2 744. The wild- 
type strain had high levels of Pgi in both glycerol and glucose 



minimal medium, indicating constitutive expression of the en- 
zyme. However, the mutant strains mc 2 743 and mc 2 744 showed 
at least a 1,000-fold reduction in Pgi activity when grown on 
glucose. In contrast, the presence of the cloned copy of the 
complementing gene restored Pgi activity to the mutant 
strains. 

As controls for the cell extracts and growth conditions, we 
measured the levels of a-glycerol phosphate dehydrogenase 
(GlpD) and triose phosphate isomerase (Tpi). The former 
should be regulated in response to a carbon source (i.e., re- 
pressed by glucose) (3), whereas the latter should be constitu- 



TABLE 2. Growth patterns of the wild type and glycerol mutants on various substrates" 



Growth pattern'' 

Carbon source (concn) 

mcrHSStpMVnOa) mc 2 743(pMVI203) mc 2 744(pMVI203) mc 2 155(pDT252) mc 2 743(pDT252) mc I 744{pDT252) 



Glucose 


+ 


+ 




+ 


+ 


+ 


Glycerol 


+ 






+ 


+ 


+ 


Glucose + glycerol 


+ 


+ 


+ 


+ 


+ 


+ 


Fructose 


+ 






+ 


+ 


+ 


Sorbitol 


+ 






+ 


+ 


+ 


Ribose 


+ 






+ 


+ 


+ 


Galactose 


+ 




± 


+ 


+ 


+ 


Gluconate 


+ 






+ 


+ 


+ 


Succinate 


+ 






+ 


+ 


+ 


Succinate + glycerol 


+ 






+ 


+ 


+ 


Casamino Acids 


+ 






+ 


+ 


+ 


Glucose (0.1%) 


-t- 








ND 


ND 


Glucose (0.01%) 








ND 


ND 


ND 


Glycerol + glucose (0.01%) 


+ 


+ 


+ 


ND 


ND 


ND 



" Cells were streaked for single colonies on minimal basal salts plates supplemented with the indicated carbon sources at 1 0 g/liter, except as indicated. 
* Symbols: +, large single colonies; ±, small single colonies; -, no single colonies. 
c ND, not done. 
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TABLE 3. Levels of phosphoglucose isomerase, triose phosphate 
isomerase, and a -glycerol phosphate dehydrogenase activities in 
wild-type and mutant strains with and without the pgji gene carried 
on plasmid pMV!203° 

Activity* 



I PTG-AG linker | 
t 

dTDP-Rhamnose 





Pei 




Toi 


GlpD 


Glycerol grown 










mc 2 155(pMVI203) 


59.6 ± 0.5 


16.3 


±0.2 


12.4 ± 0.4 


mc 2 743(pMVI203) 
mc 2 744(pMVI203) 


NG" 




NG 


NG 


NG 




NG 


NG 


mc 2 155(pDT252) 
mc 2 743(pDT252) 


120.9 ± 0.1 


9.8 


±0.03 


12.5 ± 0.6 


81.0 ± 0.8 


14.5 


±0.1 


12.4 ± 0.4 


mc 2 744(pDT252) 


78.2 ± 0.2 


7.9 


±0.05 


0.9 ± 0.1 


Glucose grown 










mc 2 155(pMVI203) 


51.8 ± 0.8 


10.8 


±0.2 


<0.6 


mc 2 743(pMVI203) 


<0.6 


8.2 


±0.2 


<0.6 


mc 2 744(pMVI203) 
mc 2 155(pDT252) 


<0.6 


8.5 


±0.3 


<0.6 


141.7 ± 0.3 


11.3 


±0.3 


<0.6 


mc 2 743(pDT252) 
mc 2 744(pDT252) 


95.2 ± 0.5 


11.0 


±0.9 


<0.6 


82.4 ± 0.3 


9.4 


± 0.1 


<0.6 



" Assays were performed three times on two extracts. 

h Activity is expressed in units per milligram of cell extract. 

c NG, no growth (mutants do not grow in glycerol minima! medium). 



tively expressed. Tpi is found at similar levels in extracts from 
all six strains (Table 3). As expected, GlpD levels were barely 
detectable in glucose-grown cells of all six strains, suggesting 
that the mutants bear no severe regulatory defects with respect 
to glycerol and glucose regulation of glycerol enzymes (3). 
However, mc 2 744(pDT252), which exhibits wild-type levels of 
Pgi activity, does show slightly reduced levels of GlpD activity 
when grown on glycerol (Table 3). The basis of this sole dif- 
ference between the two pgi mutants is under examination. 

Pgi mutants of M. smegmatis are glucose autotrophs. Since 
the pgi mutants of M. smegmatis grow only on glucose, we 
tested the possibility that they actually require glucose for 
growth. We supplied extremely low concentrations of glucose 
in media composed of other carbon sources. Glucose at a 
concentration of 0.01% is not sufficient to support the growth 
of wild-type M. smegmatis in the absence of another carbon 
source (Table 2). When this low concentration of glucose was 
added to glycerol (Table 2) or to the other carbon sources 
listed in Table 2, mc 2 743 and mc 2 744 grew at wild-type levels in 
every case (data not shown). This experiment demonstrates 
that low levels of glucose are required for growth of mycobac- 
terial pgi mutants on substrates other than glucose, indicating 
an absolute requirement in these mutants for glucose or G6P. 

DISCUSSION 

We have isolated mutants of M. smegmatis lacking a func- 
tional gene encoding the enzyme phosphoglucose isomerase 
(pgi). These mutants can grow only on glucose or galactose as 
a sole carbon source. However, any carbon source will support 
wild- type levels of growth if glucose is provided in trace 
amounts. We used the mutants to clone the gene encoding pgi 
and showed that its sequence is highly homologous to that of 
the pgi genes of other species. We measured the levels of Pgi 
in the wild type and mutants and showed that there is no 
detectable activity of Pgi in the mutants. Complementation^ 
the mutants with the cloned copy of pgi restored both intra- 
cellular levels of Pgi and the wild-type growth phenotype when 
tested on a large number of carbon sources. Finally, we dem- 
onstrated that the pgi mutants are glucose auxotrophs, since (i) 



i 



dTDP-glucose 



j GALACTAN | 

t 

UOP-galactose 

t 
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^ Glucose- I-phosphate m 

t 

Glucose-6-phosphate 

t 

Fructose-6-phosphate 

FIG. 4. Relevant cell wall synthetic pathways. Enzymes whose activity has 
been identified in mycobacteria are indicated (see Discussion). 



the mutants grow only on glucose as the sole carbon source and 
(ii) supplementation of carbon sources other than glucose with 
trace amounts of glucose restores wild-type growth. 

The absolute requirement for glucose by pgi mutants of M. 
smegmatis may be explained by examining the structure of the 
mycobacterial cell wall and its synthetic pathways. The major 
polysaccharide is arabinogalactan (AG), a large polymer posi- 
tioned by covalent linkage between the peptidoglycan (PTG) 
and mycolic acid (MA) (18). MA, constituting close to 50% by 
weight of the cell wall mass, is likely responsible for the noto- 
rious impermeability of the mycobacterial cell. The ability to 
synthesize any of the components of the PTG-AG-MA net- 
work, or cell wall core, is expected to be essential, and the 
enzymes involved in its synthesis are targets of current or 
future antimycobacterial drugs (17), 

A key feature in this covalent polymeric network is the linker 
disaccharide phosphate which links the PTG to the AG (16). 
The structure of the PTG-AG linker is L-rhamnose-(l-3)-D-yV- 
acetylglucosamine-(l-P-6)-D-muramic acid (16). As illustrated 
in Fig. 4, G6P, via glucose- 1 -phosphate, is the precursor for the 
L-rhamnose component of the linker molecule. The donor of 
this component, dTDP-rhamnose, is synthesized from glucose- 
1 -phosphate and TTP (18) in the same manner as it is in 
gram-negative bacteria, which use dTDP-rhamnose to insert 
rhamnose into lipo polysaccharide (LPS) O antigens (39). 
Thus, inability to synthesize G6P will inhibit cell wall core 
formation and lead to loss of viability of the cell. 

A second component of the mycobacterial cell wall requiring 
G6P for its synthesis is the galactan residue of AG. Again, by 
analogy with other polysaccharide synthetic systems (33, 34), 
G6P is converted to G1P and thence to UDP-glucose; an 
activity has been identified in mycobacteria which converts 
UDP-glucose to UDP-galactose for transfer into the AG com- 
plex (16). Inhibition of AG synthesis is lethal for the mycobac- 
terial cell and is most likely the mechanism of action of etham- 
butol (29). 

Pgi is a key enzyme in the pathway of cell wall biosynthesis 
in mycobacteria. Our observation that 0.01% glucose (10~ 5 
g/ml) supplies sufficient precursor for cell wall core synthesis is 
supported by calculations of the amount of carbon directed 
into L-rhamnose and galactan from G6P, assuming that there is 
1 rhamnose molecule for every 33 galactose molecules and that 
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the cell wall galactan is 1% of the dry weight of the cell 
(approximately 3 x 10~ 6 g/ml of saturated culture) (17a, 21). 

Phosphoglucose isomerase mutants have been well charac- 
terized in a number of other genetic systems. The salient phe- 
notypes of Bacillus subtilis (25) and E. colt (35) mutants lacking 
phosphoglucose isomerase are (i) very slow growth on glucose 
and galactose, (ii) normal growth on carbon sources entering 
glycolysis below the block, such as fructose and glycerol, and 
(iii) normal growth on gluconate, which is metabolized via 
either the Entner-Doudoroff pathway or the hexose mono- 
phosphate shunt. The slow growth on glucose is presumed to 
be the result of the inefficiency of growth on these latter path- 
ways, pgi mutants of E. colt can grow on gluconate as the sole 
carbon source because they do not require G6P or any of its 
derivatives for growth (8, 18); gluconate is metabolized via the 
pentose phosphate or Entner-Doudoroff pathways (Fig. 3). 
The glucose auxotrophy of the pgi mutants of M. smegmatis 
described here predicts that these mutants cannot grow on 
gluconate, since the enzymes leading from G6P to gluconate- 
6-phosphate are unidirectional (Fig. 3). Table 2 shows that this 
is indeed the case. 

In contrast, pgil mutants of S. cerevisiae cannot grow on 
glucose as the sole carbon source, yet they are glucose auxo- 
trophs (1). pgil null mutants of S. cerevisiae have been well 
characterized. These strains cannot grow on glucose, presum- 
ably because flux through the hexose monophosphate shunt is 
inadequate to support growth (5). Yet, pgil null cells also 
require glucose for growth and sporulation on other substrates 
(1 ). Indeed, the locus is also known as CDC30, since cdc30 and 
pgil cells arrest late in nuclear division at the nonpermissive 
temperature (7). A number of CDC genes encode intermediary 
metabolic genes. For example, CDC19 encodes a temperature- 
sensitive pyruvate kinase (20). The sensitivity of 5. cerevisiae 
pgil null strains to high concentrations of glucose may explain 
the cell cycle phenotype as the result of accumulation of G6P 
in arrested cells, since many phosphorylated intermediates are 
toxic when allowed to accumulate (1). 

In addition to its involvement in yeast cell cycle control, 
phosphoglucose isomerase has recently been implicated in 
growth phase gene regulation in E. coli. The pgi gene was 
identified in a screen for insertions which affect in trans the 
expression of a lacZ fusion to osmY, a gene regulated by the 
stationary-phase sigma factor RpoS (2). The authors of the E. 
coli study propose that UDP-glucose, a derivative of G6P, is an 
intracellular signal molecule in the control of expression of a 
number of ;paS-regulated genes. While the mycobacterial rpoS 
itself has not been cloned, a mycobacterial homolog of sigF (B. 
subtilis and Streptomyces coelicolor) and sigB (B. subtilis) has 
recently been cloned (13) and the regulatory targets of this and 
other putative stationary-phase and stress response regulatory 
molecules will soon be identified. It will be interesting to see 
whether pgi mutants of M. smegmatis are similarly altered in 
expression of rpoS- regulated genes, and these experiments are 
under way. 

That E. coli and S. typhitnurium pgi null mutants can grow, 
albeit slowly, on glucose suggests that there is no requirement 
for G6P and its derivatives in- these species (35). Thus, pgi 
mutants of 5. typhimurium grown on gluconate have incom- 
plete LPS, since G6P cannot be formed from gluconate (9). In 
pgi mutants of E. coli, complete polymers of LPS or glycogen 
are produced only in the presence of glucose, and this feature 
can be exploited in the study of LPS or glycogen synthesis (10). 
Similarly, the mutants described here could serve as useful 
labeling strains. M. smegmatis appears to have an absolute 
requirement for glucose or a derivative thereof: growth of the 
cells on succinate or glycerol with trace labels of high-specific- 



activity [ 14 C]glucose would direct radiolabel only into those 
lipid moieties derived from glucose. We have initiated the 
construction of pgi deletion mutants to investigate the roles of 
pgi in carbohydrate metabolism and cell wall synthesis in the 
pathogenic slow-growing species such as M. tuberculosis and 
Mycobacterium bovis. 
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